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]]

