domingo, 1 de abril de 2012

Un mono en ASCII y más ejemplos con R





Ejemplo de visionado y borrado de elementos.

> objects()
[1] "a" "b" "col" "et2" "etiquetas"
> rm(b)
> objects()
[1] "a" "col" "et2" "etiquetas"
> #le preguntamos qué es lo que hay, después borramos "b"

Ejemplo de indexación de matriz.

En la primera nos da la columna 2 de la matriz a. En la segunda instrucción nos da la fila 1 de la matriz a ( se puede pedir también un dato concreto conminando la fila y la columna).

> a
Aseo Amenaza
Sujeto A 0.451 0.16
Sujeto B 0.348 0.50
Sujeto C 0.125 0.20
Sujeto D 0.076 0.14
> a[,2]
Sujeto A Sujeto B Sujeto C Sujeto D
0.16 0.50 0.20 0.14
> a[1,]
Aseo Amenaza
0.451 0.160

Otro ejemplo para tratar los datos:

> agre<-matrix(scan("c:/r/agrefabre.txt")) Read 12 items 
> agre
[,1]
[1,] 2.05
[2,] 1.65
[3,] 0.68
[4,] 0.90
[5,] 15.38
[6,] 2.23
[7,] 1.66
[8,] 23.77
[9,] 0.61
[10,] 1.79
[11,] 7.44
[12,] 0.18


# Aquí le preguntamos los datos de 1 a 3.
> agre[1:3]
[1] 2.05 1.65 0.68
> a<-agre[1:3] 
> b<-agre[4:6] 
> c<-agre[7:9] 
> d<-agre[10:12] 
#Aquí hemos asignado nombre a unos vectores de datos que corresponden a 4 sujetos 
> t<-c(a,b,c,d) 
#Aquí el total 
> t
#Llamamos o invocamos a t y nos responde esto:
[1] 2.05 1.65 0.68 0.90 15.38 2.23 1.66 23.77 0.61 1.79 7.44 0.18


En R podemos realizar calculos de forma rápida, sin tener que hacer tantas filigranas como en el Excel (que parace más sencillo).


> range(t)
[1] 0.18 23.77
#El rango de todos lo suejtos
> max(range(t))-min(range(t))
[1] 23.59
#La amplitud del rango. Máximo del rango menos el mínimo.
> ra<-max(range(a))-min(range(a)) 
> ra
[1] 1.37
> rt<-max(range(t))-min(range(t)) 
> rt
#Les ponemos un aetiqueta para poder usar los datos más comodamente
[1] 23.59
> ra/rt
[1] 0.05807546
# Voila! La proporción del sujeto A respecto al rango del grupo.


Aquí un gráfico de 0 a 1.


> plot(0,0 , xlim=c(0, 1), ylim=c(0, 1))




Y material de la UB sobre R.
http://www.ub.edu/stat/docencia/EADB/Curso%20basico%20de%20R.htm#material

Un elegante "tutorial" sobre R.
http://www.statmethods.net/index.html

1 comentario:

Amaya dijo...

¡Muy bueno el mono! XDDD
Eso como portada de tesis también podría estar genial, jajajajaaaa!