search for: logistico

...datos.2 = data.frame( dosis.2 = c(1,10,10,10,50,50,50,100,100,100,1000,1000,1000), vivos.2 = c(20,18,18,18,7,6,6,1,1,2,0,0,0), muertos.2 = c(0,2,2,2,13,14,14,19,19,18,20,20,20) ) attach(datos.1) attach(datos.2) y1 = cbind(vivos.1,muertos.1) y2 = cbind(vivos.2,muertos.2) modelo.logistico.1 = glm(y1~log(dosis.1),binomial) modelo.logistico.2 = glm(y2~log(dosis.2),binomial) library(MASS) log.LD50.1 <- dose.p(modelo.logistico.1,p=0.5) log.LD50.2 <- dose.p(modelo.logistico.2,p=0.5) log.LD50.1 log.LD50.2 LD50.1 <- exp(log.LD50.1[1]) LD50.2 <- exp(log.LD50.2[1]) LD50....

...tstrap # ------------ B <- 1000 # numero de muestras bootstrap d.sample <- replicate(B, { i <- sample(NROW(y1), NROW(y1), TRUE) datos.1.i <- datos.1[i, ] datos.2.i <- datos.2[i, ] y1 = with(datos.1.i, cbind(vivos.1,muertos.1)) y2 = with(datos.2.i, cbind(vivos.2,muertos.2)) modelo.logistico.1 = glm(y1 ~ log(dosis.1),binomial, data = datos.1.i) modelo.logistico.2 = glm(y2 ~ log(dosis.2),binomial, data = datos.2.i) log.LD50.1 <- dose.p(modelo.logistico.1,p = 0.5) log.LD50.2 <- dose.p(modelo.logistico.2,p = 0.5) LD50.1 <- exp(log.LD50.1[1]) LD50.2 <- exp(log.LD50.2[1]) LD5...

...data as totalsinatipicos. I made selfstarting each nonlinear model. ###Library library(NRAIA) ###Data d<-totalsinatipicos f=unique(d$TRAT) dia=unique(d$DDT) REG=4:15 ###Models monomolecular=function(x,a,b,c){return(a*(1-exp(-b*(x-c))))} gompertz=function(x,a,b,c){return(a*exp(-exp(-b*(x-c))))} logistico=function(x,a,b,c){return(a/(1+exp(-b*(x-c))))} gauss=function(x,a,b,c){return(a*exp((-((x-c)*(x-c))/(2*(b*b)))))} richards=function(x,a,b,c){return(a*(1+(m-1)*exp(-b*(x-c)))**(1/(1-m)))} ###Variables #Altura ALTURA0=subset(d,d$TRAT==f[1],select=c(1,REG[1])) ALTURA_0=data.frame(dia,tapply(ALTURA0[,2...

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

...e 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...

.../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 en hotmail.com&gt...

...> 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 en hotmail.com> > > Hola, > > Estoy tratan...

.../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> &...

...ncia 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, progra...

...//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 compañeros de la lista. Tengo el siguiente set de datos: Repeticiones <- c(rep("RI", 14), rep("RII", 14), rep("RIII", 14)) Tiempo <- rep(c(0, 2, 4, 6, 8, 10, 12, 24, 36, 48, 60, 72, 96, 120), 3) Concentracion_celular <- c(0.4862, 0.5375, 0.4309, 0.4390, 0.4603, 0.4733, 0.3936, 0.9085, 0.5838, 0.5477, 0.6331, 0.8693, 1.0092, 0.6341, 0.5350,