pc4-partb

1
b) De las 1000 medias obtenidas promedio, desviación estándar Q5, Q25, Q50, Q75, Q90. >library(MASS) >edad=(Aids2$age) >MediaMuestra <- function(i){ >mean(sample(edad, size=i))} >i <- rep(500, 1000) >VMedias <- vapply(i, FUN= MediaMuestra, c(Media=0)) n=50 > mean(VMedias) [1] 37.39142 > sd(VMedias) [1] 1.377855 > quantile(VMedias,c(0.05,0.25,0.50,0.75,0.95)) 5% 25% 50% 75% 95% 35.199 36.440 37.400 38.380 39.640 n=150 mean(VMedias) [1] 37.40029 > sd(VMedias) [1] 0.7720011 > quantile(VMedias,c(0.05,0.25,0.50,0.75,0.95)) 5% 25% 50% 75% 95% 36.14000 36.86667 37.38667 37.89500 38.71367 n=500 > mean(VMedias) [1] 37.38992 > sd(VMedias) [1] 0.394892 > quantile(VMedias,c(0.05,0.25,0.50,0.75,0.95)) 5% 25% 50% 75% 95% 36.7159 37.1320 37.3830 37.6680 38.0460

Upload: katherinediazmatta

Post on 15-Sep-2015

215 views

Category:

Documents


0 download

DESCRIPTION

resistencia

TRANSCRIPT

b)De las 1000 medias obtenidas promedio, desviacin estndar Q5, Q25, Q50, Q75, Q90.

>library(MASS)>edad=(Aids2$age)>MediaMuestra mean(sample(edad, size=i))}>i VMedias mean(VMedias)[1] 37.39142> sd(VMedias)[1] 1.377855> quantile(VMedias,c(0.05,0.25,0.50,0.75,0.95)) 5% 25% 50% 75% 95% 35.199 36.440 37.400 38.380 39.640

n=150

mean(VMedias)[1] 37.40029> sd(VMedias)[1] 0.7720011> quantile(VMedias,c(0.05,0.25,0.50,0.75,0.95)) 5% 25% 50% 75% 95% 36.14000 36.86667 37.38667 37.89500 38.71367

n=500

> mean(VMedias)[1] 37.38992> sd(VMedias)[1] 0.394892> quantile(VMedias,c(0.05,0.25,0.50,0.75,0.95)) 5% 25% 50% 75% 95% 36.7159 37.1320 37.3830 37.6680 38.0460