Dear Henrique Dallazuanna,

thank you very much, this really helped me.

by the way, do you thinks it is also possible to do the following loop in R?

thanks and best regards

Andreas




v1 <- rnorm(10)
v2 <- v1+rnorm(10)
v3 <- v2+rnorm(10)
v4 <- rnorm(10)
y <- ifelse(v1<0,1,0)
dat <- data.frame(v1=v1,v2=v2,v3=v3,v4=v4)

newdat<-dat[9:10,]
dat <- dat[1:8,]
newdat[2,1] <- NA
pred <- numeric(2)
vars <- c("v1","v3","v4")
y <- y[1:8]

for (i in 1:nrow(newdat))
{
   ind <- which(!is.na(newdat[i, ]) & colnames(newdat) %in% vars)
   vars.dat <- colnames(dat[ind])
   formu <- as.formula(paste("cbind(y, 1-y) ~ ", paste(vars.dat,
         collapse = "+")))
   glm_tmp <- glm(formu, family = binomial(link = logit), data = dat[, ind])
   pred[i] <- predict(glm_tmp, newdat = newdat[i, ind],
     type = "response")
}




On 14.05.2010 14:42, Henrique Dallazuanna wrote:
> Try this:
>
> mapply(function(x, y)x[match(x, y)], m2, m1)
>
> On Fri, May 14, 2010 at 9:23 AM, Andreas Wittmann 
> <andreas_wittmann@gmx.de <mailto:andreas_wittmann@gmx.de>> wrote:
>
>     Dear R-users,
>
>     after trying and searching a long time i have the following question.
>     is it possible to replace to following double loop by some apply
>     calls?
>
>     #######################################################################
>
>     m1 <- data.frame(v1=factor(letters[1:5]),
>                     v2=factor(letters[2:6]),
>                     v3=factor(letters[3:7]))
>
>
>     m2 <- data.frame(v1=factor(letters[3:7]),
>                     v2=factor(letters[1:5]),
>                     v3=factor(letters[4:8]))
>
>
>     ind <- matrix(logical(nrow(m2) * ncol(m2)),
>                  nrow = nrow(m2), byrow = TRUE)
>
>     for (j in 1:ncol(m2))
>     {
>      for (i in 1:nrow(m2))
>      {
>        ind[i, j] <- any(levels(m1[, j]) == m2[i, j])
>      }
>      ind[is.na <http://is.na>(ind[, j]), j] <- TRUE
>      m2[!ind[, j], j] <- NA
>     }
>
>     #######################################################################
>
>
>
>
>     thanks and best regards
>
>     Andreas
>     --
>     GRATIS für alle GMX-Mitglieder: Die maxdome Movie-FLAT!
>
>     ______________________________________________
>     R-help@r-project.org <mailto:R-help@r-project.org> 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.
>
>
>
>
> -- 
> Henrique Dallazuanna
> Curitiba-Paraná-Brasil
> 25° 25' 40" S 49° 16' 22" O


	[[alternative HTML version deleted]]

