necesito que me manden lo mas pronto posible la macro que ejecuta la accion del
modelo lineal general 'glm', para una familia binomial, con funcion de enlace
la probit y la lo-log complementaria, bueno tengo la macro para lafuncion de
enlace logit es:
# macro logit
macro.logit <- function( y, m, X, b.0, mu, niter )
{
# Inicializaciones (BOPTIM = B0 o ) MU = y/m C5=g(y/m)=ETA
# y Observaciones de ocurrencias
# Considerar MU=PI, es decir Y/M seguido B(M,PI)/M
# Recordar V(y(NITER))=MUx(1-MU)/m
# X es la matriz de diseño
c1 <- y/m
mu <- c1
c5 <- log(c1/(1-c1)) # C5 = g(y/M) FUNCIÓN DE ENLACE
for (i in 1:niter)
{
c1 <- mu # C1 = MU
c2 <-((1/c1)*(1/(1-c1)))^2 # C2 = g'(C1)**2 (g'(MU)**2)
c3 <-c1*(1-c1)*c2/m # C3 = var(C1)*C2=(C1(1-C1)/M)*C2
=1/wkk) Depende
# de M y MU !!!
c4 <- 1/c3 # C4 = 1/C3 (=wkk) (Diagonal de W)
M6 <- as.matrix(diag(c4)) # M6 = W
c6 <- c5 # C6 = ETA (C5)
c7 <- (1/c1)*(1/(1-c1)) # C7 = g'(C1) (= g'(MU))
c8 <- y
c8 <- (c8/m)-c1 # C8 = y/m-MU
c9 <- c6+c7*c8 # C9 = C6 + C7*C8 (z=ETA+g'(MU)*(Y/M-MU))
M11 <- t(X)%*%M6 %*% as.vector(c9) # M11 = XTWz
M12 <- t(X)%*%M6 %*% X # M12 = XTWX
M13 <- solve( M12 ) # M13 = (XTWX)-1
# browser() # > n # para ejecutar siguiente línea de comandos
b.fi <- as.vector( M13 %*% M11 ) # BOPTIM = (XTWX)-1 XTWz
c5 <- as.vector(X %*% b.fi) # C5 = ETA = X*BOPTIM (=Xb)
mu <- exp(c5)/(1+exp(c5)) # MU = LINK-1(C5)
}
list(estimadors=b.fi,prediccions=mu*m)
}
b.0 <- as.vector(c(0,0))
mu <- as.vector(rep(0.1,6))
X <- as.matrix(data.frame(dobson$uns,dobson$log.x))
niter<-30
exe1 <- macro.logit( dobson$y, dobson$m, X, b.0, mu, 30 )
exe1
# LA INSTRUCCION DIRECTA EN R ES:
rexe1 <- glm(cbind(y,m-y)~log.x,data=dobson, family=binomial)
rexe1
bueno esto como un ejemplo de lo que quiero. Estare infinitamente agradecido.
__________________________________________________
[[alternative HTML version deleted]]
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
