[R-es] contorno noparametrico
Pablo Emilio Verde
PabloEmilio.Verde en uni-duesseldorf.de
Jue Jul 30 14:53:33 CEST 2009
Hola a todos,
En el ejemplo que sigue construyo una elipse de confianza del 95% basada en
una distribucion normal bivariada.
# comienzo del ejemplo
# parametros de la normal
muD <- 1
muS <- 2
sdD <- 2
sdS <- 4
rho <- 0.7
# contorno del 95%
cc95 <- sqrt(qchisq(0.95,2))
muSp <- muS + sdS * cc95 * cos(seq(0, 2*pi, 0.01))
muDp <- muD + sdD * cc95 * cos(seq(0, 2*pi, 0.01) + acos(rho))
plot(muDp, muSp, type="l")
# simulo 200 casos
library(MASS)
Sigma <- matrix(c(sdD^2, rho*sdD*sdS,
rho*sdD*sdS, sdS^2),2,2)
mu <- c(muD, muS)
y <- mvrnorm(n=200, mu, Sigma)
points(y)
# fin del ejemplo
Conoce alguien como construir la misma elipse pero con alguna de las
funciones de densidad no parametrica?
Pablo
Más información sobre la lista de distribución R-help-es