问题描述
我正在尝试使用dredge运行自动模型选择；通过使用全局模型中协变量的不同组合来生成一组候选模型。
我的模型如下：
occ_model <- occu(~ time_observations_started +
duration_minutes +
effort_distance_km +
number_observers +
protocol_type +
pland_04_deciduous_broadleaf +
pland_05_mixed_forest
~ pland_04_deciduous_broadleaf +
pland_05_mixed_forest +
pland_13_urban +
pland_08_woody_savanna +
pland_09_savanna +
pland_10_grassland,data = occ_um)
我收到此错误：
occ_dredge <- dredge(occ_model)
Fixed terms are "p(Int)" and "psi(Int)"
Error in h(simpleError(msg,call)) :
error in evaluating the argument 'x' in selecting a method for function 'diag': system is computationally singular: reciprocal condition number = 1.42463e-18
In addition: There were 50 or more warnings (use warnings() to see the first 50)
我似乎无法理解错误以及我出了什么问题，这可能是我正在使用的数据，而减小模型大小可能会使挖泥机更有效地工作？
我正在为模型使用的代码：
new("unmarkedFrameOccu",y = structure(c(0L,0L,1L,0L),.Dim = c(10L,10L),.Dimnames = list(NULL,c("y.1","y.2","y.3","y.4","y.5","y.6","y.7","y.8","y.9","y.10"))),obsCovs = structure(list(
time_observations_started = c(11.1833333333333,11.1833333333333,14.25,15,7.5,8.5,7.75,15.5,14.5,9.56666666666667,8.66666666666667,8.66666666666667),duration_minutes = c(100,100,125,30,120,70,20,20),effort_distance_km = c(0.644,0.644,0.805,1.609,2.414,1.127,3.219,0.161,0.161
),number_observers = c(1,1,2,2),protocol_type = structure(c(2L,2L,2L),.Label = c("Stationary","Traveling"),class = "factor"),pland_04_deciduous_broadleaf = c(0.266666666666667,0.266666666666667,0.2,0.241379310344828,0.258064516129032,0.451612903225806,0.451612903225806),pland_05_mixed_forest = c(0,0.517241379310345,0.1875,0.193548387096774,0.193548387096774)),class = "data.frame",row.names = c(NA,-100L)),siteCovs = structure(list(n_observations = c(10L,10L,latitude = c(36.3211505,35.9563584,35.8438348,35.848175,37.3624147,37.2937888,37.2723415,37.2738906,35.496358,37.2663858),longitude = c(-88.2169586,-88.9985061,-88.8894305,-88.768306,-76.7021227,-76.6597652,-76.7006765,-76.5842342,-89.2608529,-76.6323209),0.451612903225806
),0.193548387096774),pland_01_evergreen_needleleaf = c(0,0),pland_02_evergreen_broadleaf = c(0,pland_13_urban = c(0,0.387096774193548,0.03125,pland_12_cropland = c(0,0.9,0.032258064516129,pland_08_woody_savanna = c(0.7,0.4,0.733333333333333,0.483870967741935,0.125,0.290322580645161),pland_09_savanna = c(0.0333333333333333,0.366666666666667,0.0666666666666667,0.129032258064516,0.0645161290322581),pland_10_grassland = c(0,0.0967741935483871,0)),class = "data.frame"),mapInfo = NULL,obsToY = structure(c(1,1),10L
)))
解决方法
暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！
如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@）