Displaying 2 results from an estimated 2 matches for "mt_yop_new".
2010 Sep 28
3
calcular la variancia de gini por bootstrap
...tencion q la funcion q definia en el ejmeplo tenia 2
parametros d y w, este es el ejemplo del manual modificado para mi base:
ratio <- function(d, w) sum(d$EDAD * w)/sum(d$UNO * w)
Como ves tiene 2 parametros: d y w, entonces le puse eso mismo a mi funcion
grini: grini<-function(d,w){gini(d$MT_YOP_NEW, w=d$FEXP)}, donde el w q esta
adentro no es el w del parametro de la funcion sino el el parametro peso de
la funcion gini.
Me hizo el bootstrap pero me da siempre la misma estimacion para todas las
replicaciones de las muestras, por lo cual la variancia del estimador por
bootstrap da 0.
Entonces...
2010 Sep 29
0
Resumen de R-help-es, Vol 19, Envío 28
...t;*w" al final de cada parametro en la definicion
de la funcion "grini" q habia creado, cosa q deduje de mirar la funcion q le
habia robado al ejemplo del manual de boot
ratio <- function(d, w) sum(d$EDAD * w)/sum(d$UNO * w)
En tonces mi funcion paso de ser:
function(d,w){gini(d$MT_YOP_NEW, w=d$FEXP)
a ser:
grini<-function(d,w){gini(d$MT_YOP_NEW**w*, w=d$FEXP**w*)}
Como les decia este cambio hizo q el programa ande, (deje de darme variancia
cero) pero no se cual es su significado. Ahora el tema es q los resultados
son algo distintos del gini original, pero bueno sobre eso seguire...