[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