Pessoal,

Estou com uma dúvida sobre a aplicação da função density, em um banco de dados de uma distribuição normal. Após fazer uma transformação encontrando o score (componente principal), ao aplicar a função density, eu esperava que os valores deveriam variar de 0 a 1. Plotei até o histograma da densidade. Por que a densidade não está no intervalo 0 a 1?

Gerei um código reproduzível, porém, minha variável que estou aplicando, não necessariamente é exatamente uma variável normal. 

set.seed(23)

X = cbind( 1:230, rnorm(230, 450, 20))

X[,2] <- (X[,2] - min(X[,2]))/(max(X[,2])-min(X[,2]))
#plot(X)
X[,1] <- X[,1]-mean(X[,1])
X[,2] <- X[,2]-mean(X[,2])
C <- cov(X)
E <- eigen(C)
score <- X %*% E$vectors
densidade = density(score[,2], n=100)  # os valores da densidade não estão dentro do intervalo [0,1]

a = hist(score[,2], freq = F)
lines(densidade)

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.