[R] ayuda con macros

From: Dony Henry Clavel Quijada <donyhenryc_at_yahoo.com>
Date: Fri, 18 May 2007 15:18:22 -0500 (CDT)


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



R-help_at_stat.math.ethz.ch 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. Received on Fri 18 May 2007 - 20:30:54 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Fri 18 May 2007 - 21:31:48 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.