R help es - Jun 2013 - calculo de poder estadistico en glm

Hola,
 
Estoy tratando de calcular el poder estadistico de un GLM en R, programas
estadisticos como SPSS lo puede calcular pero no eh encotrado nada sobre R.
Alguien tiene alguna idea o sabe de literatura donde pueda encontrar ayuda?
 
Luego tengo otra pregunta mas bien teorica sobre la funcion offset utilizada en
negative poisson, por lo que eh entendido esta funcion posibilita trabajar con
fracciones a pesar de que la distribucion poisson solo sea para enteros
positivos. Mi duda es si esto es realmente asi, lo que eh leido a sido en
internet. Eh buscado en libros como Zuur 2007 y Hilbe 2007 sin encontrar nada.
 
modelA1F<-glm.nb(Tot.dameggs~STATN+YEAR+offset(log(FECUND)))
 
freq=(Tot.dameggs/FECUND)
 
Tot.dameggs= cantidad de huevos dañados por individuo
FECUND=cantidad total de huevos por individuo
STATN= estacion 
YEAR= año 
 
Muchas gracias,
Matias
 
 
 
 		 	   		  
	[[alternative HTML version deleted]]

Hola Matias,

Que tipo de GLM estas usando?  Cuantas variables explicativas?  Que tamaño
de efecto te interesa detectar?

Una forma de realizar el analisis de potencia es via simulacion.  Dale una
mirada a
http://r.789695.n4.nabble.com/a-priori-power-analysis-for-glm-family-poisson-td4659163.htmly
en especial a la referencia que alli aparece.  El paquete asypow
tambien
puede servirte.

En cuanto a tu 2a pregunta, en mi opinion, lo que necesitas es un modelo
logistico, no un modelo binomial negativo. En
http://ww2.coastal.edu/kingw/statistics/R-tutorials/logistic.html hay un
ejemplo parecido a lo que planteas.

Saludos,
Jorge.-


2013/6/11 Matias Ledesma <matutetote@hotmail.com>
> Hola,
>
> Estoy tratando de calcular el poder estadistico de un GLM en R, programas
> estadisticos como SPSS lo puede calcular pero no eh encotrado nada sobre R.
> Alguien tiene alguna idea o sabe de literatura donde pueda encontrar ayuda?
>
> Luego tengo otra pregunta mas bien teorica sobre la funcion offset
> utilizada en negative poisson, por lo que eh entendido esta funcion
> posibilita trabajar con fracciones a pesar de que la distribucion poisson
> solo sea para enteros positivos. Mi duda es si esto es realmente asi, lo
> que eh leido a sido en internet. Eh buscado en libros como Zuur 2007 y
> Hilbe 2007 sin encontrar nada.
>
> modelA1F<-glm.nb(Tot.dameggs~STATN+YEAR+offset(log(FECUND)))
>
> freq=(Tot.dameggs/FECUND)
>
> Tot.dameggs= cantidad de huevos dañados por individuo
> FECUND=cantidad total de huevos por individuo
> STATN= estacion
> YEAR= año
>
> Muchas gracias,
> Matias
>
>
>
>
>         [[alternative HTML version deleted]]
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>
	[[alternative HTML version deleted]]