Dear R-helpers,
I have run the code below which I expected to make an object called dd1,
but that object does not exist.
So, in summary, my problem is that my function is meant to make an object
(dd1), and it does indeed make that object (I know that the last line of
the function prints it out) but then, after the function has run, the
object has disappeared.
It's late on a Friday so I may be overlooking something obvious, but I'd
appreciate your help if you can see what I'm doing wrong.
Many thanks,
Mark Na
rm(list=ls())
#LOAD package, data, and keep clean copies of data
library(reshape2); library(MuMIn);library(AICcmodavg)
birdmatrix<-read.csv("birdmatrix.csv",stringsAsFactors=FALSE);
birdmatrixclean<-birdmatrix
habitatmatrix<-read.csv("habitatmatrix.csv",stringsAsFactors=FALSE);
habitatmatrixclean<-habitatmatrix
traits<-read.csv("traits.csv",stringsAsFactors=FALSE); traitsclean<-traits
birdhabitat<-merge(birdmatrix,habitatmatrix,by="SARBASINYEAR") #merge bird
and habitat data
data<-birdhabitat; data$SAR<-substr(data$SARBASINYEAR,1,3) #convenience;
create SAR variable
#make a FUNCTION to combine variables into models
dredgeit<-function(lm1){
dd1<-dredge(lm1,eval=FALSE) #create all possible combinations
dd1<-gsub("formula = ","",dd1,fixed=TRUE) #delete characters
dd1<-gsub(" + 1","",dd1,fixed=TRUE) #delete characters
dd1#inspect model formulae
}
lm1<-lm(data$BiRich.o~log(data$HaArea,10)+log(data$HaPeri,10)) #saturated
submodel
dredgeit(lm1)
[[alternative HTML version deleted]]
______________________________________________
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.
