lunes, 6 de junio de 2016
domingo, 5 de junio de 2016
Verificación de normalidad univariado
### Verificación de normalidad - univariado
# Gráfica qq-plot
set.seed(555)
dat<-as.vector((rnorm(100,0,1)))
datos<-sort(datos)
head(datos)
n<-length(datos)
largo<-1:n
pi<-c((1:n)-0.5)/n
qn<-qnorm(pi)
tabla<-cbind(datos,pi,qn)
tabla
par(mfrow=c(1,2))
plot(qn,datos, main="qqplot-manual")
qqnorm(datos, main="qqplot-fórmula")
# Gráfica qq-plot
set.seed(555)
dat<-as.vector((rnorm(100,0,1)))
datos<-sort(datos)
head(datos)
n<-length(datos)
largo<-1:n
pi<-c((1:n)-0.5)/n
qn<-qnorm(pi)
tabla<-cbind(datos,pi,qn)
tabla
par(mfrow=c(1,2))
plot(qn,datos, main="qqplot-manual")
qqnorm(datos, main="qqplot-fórmula")
Suscribirse a:
Entradas (Atom)