Displaying 5 results from an estimated 5 matches for "treecov".
2009 Jan 26
1
glm StepAIC with all interactions and update to remove a term vs. glm specifying all but a few terms and stepAIC
...are 78 predictor terms
considered, the 12 primary terms and 66 interactions [n(n+1)/2]. I see this
with trace=T also. Here is the code used:
>glm1<-glm(formula = PRESENCE == "1" ~ SNOW + I(SNOW^2) + POP_DEN + ROAD_DE
+ ADJELEV + I(ADJELEV^2) + TRI + I(TRI^2) + EDGE + I(EDGE^2) + TREECOV +
I(TREECOV^2),family = binomial, data = wolv)
summary(glm1)
>library(MASS)
>stepglm2<-stepAIC(glm1,scope=list(upper=~.^2,lower=NULL),
trace=T,k=log(4828),direction="both")
> summary(stepglm2)
> extractAIC(stepglm2,k=log(4828))
This results in a 15 term m...
2009 Jul 13
2
Help me get this function to work...
...for (i in 1:751){
if (TotalCover.df[i,"Station"]>=ShrubCover.df[j,"Start"] && TotalCover.df[i,"Station"]<= ShrubCover.df[j,"Stop"])
TotalCover.df[i,"Shrub"]<- 1
}
}
TreeCover.df<-Tree.df[Tree.df$Quadrat==ID,]
for (j in 1:length(TreeCover.df[,"Quadrat"])){
for (i in 1:751){
if (TotalCover.df[i,"Station"]>=TreeCover.df[j,"Start"] && TotalCover.df[i,"Station"]<= TreeCover.df[j...
2009 Jul 13
1
Help get this simple function to work...
...for (i in 1:751){
if (TotalCover.df[i,"Station"]>=ShrubCover.df[j,"Start"]
&& TotalCover.df[i,"Station"]<= ShrubCover.df[j,"Stop"])
TotalCover.df[i,"Shrub"]<- 1
}
}
TreeCover.df<-Tree.df[Tree.df$Quadrat==ID,]
for (j in 1:length(TreeCover.df[,"Quadrat"])){
for (i in 1:751){
if (TotalCover.df[i,"Station"]>=TreeCover.df[j,"Start"]
&& TotalCover.df[i,"Station"]<= TreeCover.df[j...
2012 Jan 11
2
Vegan(ordistep) error: Error in if (aod[1, 5] <= Pin) { : missing value where TRUE/FALSE needed
...is(expr, envir, enclos)
2: eval.with.vis(ei, envir)
1:
source("~/Documents/Dropbox/thesis/CH3/Analysis/RDAPARTIALSexandAgeConnectandGEOGraphy2.R")
print(myrda1)
Call: rda(formula = mygenind@tab ~ pcnmTRE_25_100_CS25 + pcnmTRE_25_10_CS25
+ pcnmTRE_25_2_CS25 +
pcnmTRE_25_5_CS25 + mydata$TreeCov + mydata$Hab_Config +
pcnmEYR_EO_100_CS25 +
pcnmEYR_EO_5000_CS25 + pcnmEYR_TH_10_CS25 + pcnmEYR_TH_2_CS25 +
mydata$Site_No + mydata$Landscape
+ Condition(pcnmCS_NULL + mydata$LAT.x + mydata$LONG.x), na.action =
"na.omit")
Inertia Proportion Rank
Total 1.8110 1.0...
2012 Jan 10
0
Error message in vegan ordistep
...r, enclos)
2: eval.with.vis(ei, envir)
1: source("~/Documents/Dropbox/thesis/CH3/Analysis/RDAPARTIALSexandAgeConnectandGEOGraphy2.R")
>
> print(myrda1)
Call: rda(formula = mygenind at tab ~ pcnmTRE_25_100_CS25 + pcnmTRE_25_10_CS25 + pcnmTRE_25_2_CS25 +
pcnmTRE_25_5_CS25 + mydata$TreeCov + mydata$Hab_Config + pcnmEYR_EO_100_CS25 +
pcnmEYR_EO_5000_CS25 + pcnmEYR_TH_10_CS25 + pcnmEYR_TH_2_CS25 + mydata$Site_No + mydata$Landscape
+ Condition(pcnmCS_NULL + mydata$LAT.x + mydata$LONG.x), na.action = "na.omit")
Inertia Proportion Rank
Total 1.8110 1....