setwd("~/Dropbox/Clases_UPJ/Biología de poblaciones/Alumnos_BioP/2018_30/Clases/Práctica")
datos1<-read.table("ma2.csv", sep = ",", header=T) # carga los datos
datos1 # observar los datos, se desea saver si hay homogenedad de varianzas para la masa en función de la densidad (debe ser un factor, es decir, que agrupa)
## densidad masa
## 1 1566 0.1481
## 2 1566 0.1954
## 3 1562 0.2204
## 4 1562 0.1207
## 5 1542 0.2920
## 6 1542 0.3892
## 7 1540 0.3837
## 8 1540 0.2994
## 9 1538 0.2892
## 10 1538 0.1347
## 11 1536 0.2254
## 12 1536 0.7405
## 13 1518 0.4069
## 14 1518 0.1213
## 15 1516 0.2477
## 16 1516 0.3641
## 17 1514 0.4780
## 18 1514 0.3243
## 19 1512 0.7489
## 20 1512 0.8562
## 21 1494 0.5270
## 22 1494 0.4817
## 23 1492 0.3918
## 24 1492 0.3449
## 25 1490 0.3349
## 26 1490 0.6067
## 27 1488 0.5888
## 28 1488 1.0516
## 29 1470 0.4802
## 30 1470 1.7754
## 31 1468 1.4030
## 32 1468 0.9107
## 33 1466 0.5893
## 34 1466 1.2156
## 35 1464 1.3592
## 36 1464 1.8306
## 37 1446 1.6102
## 38 1446 1.6628
## 39 1446 1.3939
## 40 1446 1.8325
## 41 1444 1.7633
## 42 1444 1.2809
## 43 1444 1.7984
## 44 1444 2.0020
## 45 1442 1.4276
## 46 1442 1.9924
## 47 1442 1.8784
## 48 1442 1.2773
## 49 1440 1.2808
## 50 1440 1.7115
## 51 1440 1.9192
## 52 1440 1.6359
## 53 1422 2.8007
## 54 1422 2.7597
## 55 1422 2.2568
## 56 1419 1.9971
## 57 1419 1.4329
## 58 1419 2.8216
## 59 1419 2.4322
## 60 1416 2.8792
## 61 1416 1.7876
## 62 1416 1.8798
## 63 1416 2.0045
attach(datos1) # adjunta los datos para utilizar los nombres de variables directamente
densidad=as.factor(densidad) # el test de levene mide la homogeneidad de varianzas de una serie de observaciones a lo largo de un factor (que agrupa), entonces se debe convertir una de las variables en factor, en este caso la densidad ya que esta como una variable continua (con números no letras) para R
is.factor(densidad) # con esto se puede saber si la densidad se convirtió finalmente en factor
## [1] TRUE
## Una vez se tiene el factor puedo aplicar el test de levene
library(car)
## Warning: package 'car' was built under R version 3.3.2
leveneTest(masa,densidad)
## Levene's Test for Homogeneity of Variance (center = median)
## Df F value Pr(>F)
## group 24 1.6452 0.08293 .
## 38
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1