Использование пакета glmulti в R для множественной регрессии исчерпывающего поиска для весов akaike

Мне было интересно, может ли кто-нибудь помочь мне понять, почему я получаю сообщение об ошибке при вводе сценария в R.Для получения некоторой справочной информации я изучаювлияние 6 различных переменных (, которые, как я думаю, составляют 63 комбинации или модели)(X)на валовую первичную и чистую продукцию экосистемы (Y )отдельно в разных пространственных масштабах для моего проекта по науке об окружающей среде. Я решил использовать множественный регрессионный анализ исчерпывающего поиска с информационным критерием akaikes (AIC ), чтобы попытаться найти группу моделей для наилучшего соответствия. (и иерархическое разбиение для сравнения дисперсии, приписываемой разным X-переменным )Я хочу получить веса, чтобы ранжировать модели, которые «наиболее соответствуют» критерию, посмотреть, есть ли одна или группа из них, которые дополняют остальные и поэтому лучше подходит для данных.

Недавно я разместил аналогичный вопрос о пакете hier.part на сайте Cross Validated, получил отличный ответ и посоветовал обращаться сюда, если у меня возникнут подобные вопросы в будущем.

Для R я использую пакет glmulti. который можно найти здесь

Сценарий, который я использую, таков

require(glmulti)
GPPANDDRIVER<-read.table("C:\\Databases at different scales for R\\River Rhine and Netherlands\\GPP and drivers rhineland (comma delimited).csv",header=T,sep=",")
GPP<-GPPANDDRIVER$GPP
IND_VARS<-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,)

Вот ссылка на данные.csv для сайтов в Рейнской области, упомянутых в примере,http://www.filedropper.com/gppanddriversrhinelandcommadelimited

Я очень новичок в R, поэтому Я предположил, чтоpopsizeозначает количество повторений, которое равно 40 для этой шкалы, поэтому я использовал 40, я также предположил, чтоconfsetsizeозначает количество возможных моделей, которое, как я полагаю, равно 63 из-за 6 переменных. ?

Если кто-то может помочь, буду очень признателен

Спасибо за терпение и извинения за основной вопрос

Ричард

редактировать Я только что попробовал запустить скрипт сегодня утром, и теперь R.