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")



0 comentarios:

Publicar un comentario