我想看看蝙蝠叫聲的次數和幼崽飼養季節的時間之間是否存在關系。pup 變數具有三個類別：“前”、“中”和“后”。當我索要摘要時，它只包括幼崽生產前后的 p 值。我在下面創建了一個示例資料集。使用示例資料集，我只是得到一個錯誤......使用我的實際資料集，我得到了上面描述的輸出。
樣本資料集：
Calls<- c("55","60","180","160","110","50")
Pup<-c("Pre","Middle","Post","Post","Middle","Pre")
q<-data.frame(Calls, Pup)
q
q1<-lm(Calls~Pup, data=q)
summary(q1)
來自示例的輸出和錯誤訊息：
> Calls Pup
1 55 Pre
2 60 Middle
3 180 Post
4 160 Post
5 110 Middle
6 50 Pre
Error in as.character.factor(x) : malformed factor
In addition: Warning message:
In Ops.factor(r, 2) : ‘^’ not meaningful for factors
我的分析的實際輸入：
> pupint <- lm(Calls ~ Pup, data = park2)
summary(pupint)
這是我從實際資料集中得到的輸出：
Residuals:
Min 1Q Median 3Q Max
-66.40 -37.63 -26.02 -5.39 299.93
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 66.54 35.82 1.858 0.0734 .
PupPost -51.98 48.50 -1.072 0.2927
PupPre -26.47 39.86 -0.664 0.5118
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 80.1 on 29 degrees of freedom
Multiple R-squared: 0.03822, Adjusted R-squared: -0.02811
F-statistic: 0.5762 on 2 and 29 DF, p-value: 0.5683
總的來說，只是想知道為什么上面的輸出沒有顯示“中間”。抱歉，我的示例資料集沒有得到相同的結果，但該錯誤訊息可能有助于更好地理解問題。
uj5u.com熱心網友回復：
為了讓 R 正確理解虛擬變數，您必須Pup使用factor
> Pup <- factor(Pup)
> q<-data.frame(Calls, Pup)
> q1<-lm(Calls~Pup, data=q)
> summary(q1)
Call:
lm(formula = Calls ~ Pup, data = q)
Residuals:
1 2 3 4 5 6
2.5 -25.0 10.0 -10.0 25.0 -2.5
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 85.00 15.61 5.444 0.0122 *
PupPost 85.00 22.08 3.850 0.0309 *
PupPre -32.50 22.08 -1.472 0.2374
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 22.08 on 3 degrees of freedom
Multiple R-squared: 0.9097, Adjusted R-squared: 0.8494
F-statistic: 15.1 on 2 and 3 DF, p-value: 0.02716
如果您希望 R 顯示虛擬變數內的所有類別，則必須從回歸中洗掉截距，否則，您將處于變數 dummy trap 中。
summary(lm(Calls~Pup-1, data=q))
Call:
lm(formula = Calls ~ Pup - 1, data = q)
Residuals:
1 2 3 4 5 6
2.5 -25.0 10.0 -10.0 25.0 -2.5
Coefficients:
Estimate Std. Error t value Pr(>|t|)
PupMiddle 85.00 15.61 5.444 0.01217 *
PupPost 170.00 15.61 10.889 0.00166 **
PupPre 52.50 15.61 3.363 0.04365 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 22.08 on 3 degrees of freedom
Multiple R-squared: 0.9815, Adjusted R-squared: 0.9631
F-statistic: 53.17 on 3 and 3 DF, p-value: 0.004234
uj5u.com熱心網友回復：
如果您pup在回歸中包含一個分類變數，那么它會為該變數中的每個值都包含一個虛擬變數，默認情況下除外。pupmiddle如果你省略截距系數，你可以顯示一個系數，如下所示：
q1<-lm(Calls~Pup - 1, data=q)
轉載請註明出處，本文鏈接：https://www.uj5u.com/shujuku/358113.html
上一篇：R中多列的值計數