Hi
r-help-boun...@r-project.org napsal dne 01.03.2010 16:36:27:
>
> Hi Petr,
>
> Thanks again for trying again with these data, I really appreciate it.
>
> Your script works perfectly, but the problem I'm having is how to store
the
> model results so after your script I would do:
>
> m1.R<-glm(cbind(res$BEH_T, res$BEH_F) ~ res$SITE + res$YEAR +
> res$PRECIP_MM_DAY + res$PUP_AGE_EST + res$MO_AIR_TEMP,
family="binomial")
> mod<-dredge(m1.R)
>
> where mod is a list not a vector.
>
> So your example has 10 iterations of the loop so there should therefore
be
> 10 different mod[1,] that I want to store and that is what I can't work
out
> how to do, for example I can do this:
>
> if (i>=1) print (mod[1,]) else print ("NO")}
>
> And I will get a print of each of the 10 model outputs that I want, but
I
> want to store these somewhere. I did try to adjust your value <- matrix
> section of the script but had no luck.
>
> I hope this is a little clearer?
Well, I am still a bit unsure what is your problem. You used slower way
but i believe you should use quicker one as you want to cycle it 100 times
Assuming your data are named data
dat.o<-data[order(data$ID2),]
len<-rle(dat.o$ID2)$lengths
shift.len<-c(0,cumsum(len))[-(length(len)+1)]
# declare object for storing your results (the structure depends on
structure of the result)
mod<- vector(10, mode="list")
for(i in 1:10) {
samp<-sapply(lapply(split(dat.o$ID2, dat.o$ID2), function(x) 1:length(x)),
sample, 1)
Sample.dat <- dat.o[shift.len+samp,]
# do any of your stuff here (put Sample.dat instead of res)
m1.R<-glm(cbind(res$BEH_T, res$BEH_F) ~ res$SITE + res$YEAR +
res$PRECIP_MM_DAY + res$PUP_AGE_EST + res$MO_AIR_TEMP, family="binomial")
mod[[i]]<-dredge(m1.R)
}
then you shall get a list with 10 components each one for one result.
You can extract each component acccordingly
mod[[1]]
Regards
Petr
>
> Thank you again for your help, I really appreciate it!
>
> Ross
>
>
> --
> View this message in context:
http://n4.nabble.com/Loop-overwrite-and-data-
> output-problems-tp1570593p1573703.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> 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.
______________________________________________
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.
