![]() The full dataset from the dredge can be found here and the full dataset here What does this mean? Should I be ignoring and removing the top model as that doesn't seem right and conducting the model averaging on the other models? Or is it OK to use this as the top model? The top model has a vastly different AIC to the others, but has -23 degrees of freedom GsGlob_dredge <- MuMIn::pdredge(gsGlob, mycluster) #required packages must be also loaded there #data must exported to the cluster - see 'details' ĬlusterExport(mycluster,"allVars_node_dat") mycluster = makeCluster(5, type = "SOCK") # also need snow installed S(code, bs = 're') + s(station, bs = 're'),įamily=nb(), data=allVars_node_dat, na.action = "na.fail", discrete = TRUE)Īnd I'm using pdredge from MuMIn so I can increase the speed of the dredge. Salinity2 + SST.anomaly2 + s(SST.variability2) + wind2 + gsGlob <- bam(gs~ species + season + sex + TL2 + year + s(ri, bs="ad") + I have been using the dredge function in the MuMIn package to conduct model averaging on my global GAM model (using bam from the mgcv package), with a priori selected explanatory variables and two random effects and a negative binomial distribution.
0 Comments
Leave a Reply. |
AuthorWrite something about yourself. No need to be fancy, just an overview. ArchivesCategories |