title: Using the glmulti package in R for exhaustive search multiple regression for akaike weights
body: <p>I was wondering if someone could help me understand why I am getting an error message when I enter a script into R. For abit of background information I am looking into <strong><em>the effect 6 different variables (which I think is 63 combinations or models) (X)</em></strong> have on gross primary and net ecosystem production (Y) seperatly at different spatial scales for my environmental science honours project. I have decided to use exhaustive search multiple regression analysis with akaikes information criterion (AIC) to try and find a group of models for best fit. (and hierarchical partitioning to compare variance attributed to different X variables) I want to get the weights so I can rank which models "best meet" the criterion see if there is one or a group of them that outfit the rest and therefore be a more likely fit to the data. </p>
  
  <p>I recently posted a similar question on the hier.part package on Cross Validated received a great answer and was told to come here if I had any similar questions in the future.</p>
  
  <p>The package I am using for R is glmulti. <a href="http://cran.r-project.org/web/packages/glmulti/glmulti.pdf" rel="nofollow">which can be found here</a></p>
  
  <p>The script i am using is this</p>
  
  <pre><code>require(glmulti)
  GPPANDDRIVER&lt;-read.table("C:\\Databases at different scales for R\\River Rhine and Netherlands\\GPP and drivers rhineland (comma delimited).csv",header=T,sep=",")
  GPP&lt;-GPPANDDRIVER$GPP
  IND_VARS&lt;-subset(GPPANDDRIVER,select=-GPP)
  #  glmulti S4 generic 
  glmulti(y=GPP, xr=IND_VARS, data, exclude = c(), name = "glmulti.analysis", intercept = TRUE, marginality = FALSE, bunch=30, chunk = 1, chunks = 1,
  level = 2, minsize = 0, maxsize = -1, minK = 0, maxK = -1, method = "h", crit = "aic", confsetsize = 63, popsize = 40, mutrate = 10^-3, sexrate = 0.1, imm = 0.3, plotty = TRUE, report = TRUE, deltaM = 0.05, deltaB = 0.05, conseq = 5, fitfunction = "glm", resumefile = "id", includeobjects=TRUE,)
  </code></pre>
  
  <p>Here is the link for the .csv data for sites in the rhineland mentioned in the example, <a href="http://www.filedropper.com/gppanddriversrhinelandcommadelimited" rel="nofollow">http://www.filedropper.com/gppanddriversrhinelandcommadelimited</a></p>
  
  <p>I am extremely new to R so I presumed <strong><em>popsize</em></strong> means the number of replicates which is 40 for this scale so I used 40, I also assumed <strong><em>confsetsize</em></strong> meant number of possible models which I believe is 63 due to the 6 variables?</p>
  
  <p>If anyone could help it would be greatly appreciated</p>
  
  <p>Thanks for you patience and apologies for the basic question</p>
  
  <p>Richard </p>
  
  <p><em>edit</em> I just tried running the script this morning and it now crashes R.</p>