eval(predvars, data, env) 中的錯誤：找不到對象“下水道”

Question

我有一組包含物種名稱和編號 ( spp_data ) 的數據，我正在嘗試測試物種如何受到不同參數的影響，例如 pH、電導率以及下水道位置（上游/下游）（ env_data1 ）。 當我嘗試運行lm() ，出現以下錯誤：

lm1 <- lm(specnumber ~ Sewer + pH + Conductivity, data=spp_data,env_data1)

eval(predvars, data, env) 中的錯誤：找不到對象“下水道”

是不是因為Sewer列是非數字的？ 我還嘗試排除該列並運行lm()但它不起作用。

物種數據

summary(spp_data)
    Pisidium        G_pulex         C_pseudo        A_aquatic          V_pisc       
 Min.   :0.000   Min.   : 0.00   Min.   : 0.000   Min.   :0.0000   Min.   :0.00000  
 1st Qu.:0.000   1st Qu.: 3.00   1st Qu.: 0.000   1st Qu.:0.0000   1st Qu.:0.00000  
 Median :0.000   Median : 8.00   Median : 3.000   Median :0.0000   Median :0.00000  
 Mean   :1.429   Mean   :16.86   Mean   : 4.476   Mean   :0.5714   Mean   :0.04762  
 3rd Qu.:2.000   3rd Qu.:20.00   3rd Qu.:10.000   3rd Qu.:0.0000   3rd Qu.:0.00000  
 Max.   :7.000   Max.   :68.00   Max.   :16.000   Max.   :4.0000   Max.   :1.00000  

    Taeniopt          Rhyacoph         Hydropsy        Lepidost        Glossos     
 Min.   :0.00000   Min.   :0.0000   Min.   :0.000   Min.   :0.000   Min.   : 0.00  
 1st Qu.:0.00000   1st Qu.:0.0000   1st Qu.:0.000   1st Qu.:0.000   1st Qu.: 0.00  
 Median :0.00000   Median :0.0000   Median :0.000   Median :0.000   Median : 0.00  
 Mean   :0.09524   Mean   :0.2381   Mean   :1.286   Mean   :1.238   Mean   : 1.81  
 3rd Qu.:0.00000   3rd Qu.:0.0000   3rd Qu.:3.000   3rd Qu.:2.000   3rd Qu.: 1.00  
 Max.   :2.00000   Max.   :2.0000   Max.   :5.000   Max.   :7.000   Max.   :14.00  
    Agapetus         Hydroptil          Limneph         S_person           Tipula 
 Min.   : 0.0000   Min.   :0.00000   Min.   :0.000   Min.   :0.00000   Min.   :0  
 1st Qu.: 0.0000   1st Qu.:0.00000   1st Qu.:0.000   1st Qu.:0.00000   1st Qu.:0  
 Median : 0.0000   Median :0.00000   Median :0.000   Median :0.00000   Median :0  
 Mean   : 0.5714   Mean   :0.04762   Mean   :0.381   Mean   :0.09524   Mean   :0  
 3rd Qu.: 0.0000   3rd Qu.:0.00000   3rd Qu.:1.000   3rd Qu.:0.00000   3rd Qu.:0  
 Max.   :12.0000   Max.   :1.00000   Max.   :2.000   Max.   :2.00000   Max.   :0  
    Culicida         Ceratopo     Simuliid          Chrinomi         Chrnomus     
 Min.   :0.0000   Min.   : 0   Min.   : 0.0000   Min.   : 0.000   Min.   : 0.000  
 1st Qu.:0.0000   1st Qu.: 0   1st Qu.: 0.0000   1st Qu.: 0.000   1st Qu.: 1.000  
 Median :0.0000   Median : 1   Median : 0.0000   Median : 2.000   Median : 3.000  
 Mean   :0.5714   Mean   : 7   Mean   : 0.5238   Mean   : 7.286   Mean   : 6.095  
 3rd Qu.:0.0000   3rd Qu.: 8   3rd Qu.: 0.0000   3rd Qu.: 8.000   3rd Qu.: 6.000  
 Max.   :5.0000   Max.   :31   Max.   :10.0000   Max.   :67.000   Max.   :41.000  

環境數據

summary(env_data)
    Sample             Sewer                 pH        Conductivity      
 Length:21          Length:21          Min.   :7.780   Length:21         
 Class :character   Class :character   1st Qu.:7.850   Class :character  
 Mode  :character   Mode  :character   Median :8.100   Mode  :character  
                                       Mean   :8.044                     
                                       3rd Qu.:8.270                     
                                       Max.   :8.280                     
     Depth           %rock            %mud      %sand,,         
 Min.   : 7.00   Min.   :10.00   Min.   : 0   Length:21         
 1st Qu.: 8.00   1st Qu.:10.00   1st Qu.:20   Class :character  
 Median :11.00   Median :70.00   Median :30   Mode  :character  
 Mean   :17.14   Mean   :57.14   Mean   :40                     
 3rd Qu.:28.00   3rd Qu.:80.00   3rd Qu.:90                     
 Max.   :40.00   Max.   :90.00   Max.   :90

Answer 1

假設您的spp_data行與您的環境數據行匹配......我想如果你這樣做

lm1 <- lm(as.matrix(spp_data) ~ Sewer + pH + Conductivity, 
        data=env_data1)

您將獲得運行 44 個獨立線性模型的結果，每個物種一個。 （注意：有 44 個回歸和只有 21 個觀測值，您可能需要進行一些多重比較校正以避免誇大您的結論。）

有用於更復雜的多物種分析的 R 包，例如mvabund或gllvm ，但它們可能不適用於這種大小的數據集......