審查中的R線性模型的觀測子集

-3

在R，I有一個線性模型標題lmod452。從這個線性模型中，我想檢查來自第25-50號案例的觀察結果，以觀察哪一個是異化的。如何在選定數量的觀測值上使用繪圖函數？審查中的R線性模型的觀測子集

謝謝！

2017-07-14 boccultist

您可以使用[]到子集數據。所以對於假數據

x <- rnorm(100) 
y <- x + 0.5*rnorm(100) 
plot(x,y)

...你可以繪製觀察的一個子集...

plot(x[25:50], y[25:50])

...或使用不同的符標或顏色疊加他們...

plot(x,y) 
points(x[25:50], y[25:50], pch=16, col="red")

很多路要走，希望這可以讓你開始。

哦，既然你提到你的LM模型，而不是你的數據，這是值得注意的是，你可以從一個LM對象提取數據。只需撥打names(lmod452$model)即可獲取姓名。

lmod452 <- lm(y~x) 

# extracting stuff 
xxx <- lmod452$model$x 
yyy <- lmod452$model$y

來源

2017-07-14 22:56:09

不幸的是我無法按照你的代碼來使它工作。我嘗試： > X < - lmod452 $模型$ X >ý< - lmod452 $模型$ Y >積（X [25:50]，Y [25:50]）錯誤plot.window （...）：需要有限「XLIM」值此外：警告信息： 1：在分（X）：沒有非缺失參數分鐘;返回Inf 2：在max（x）中：沒有非缺失參數爲max;返回-Inf 3：在min（x）中：沒有非缺少參數min;返回Inf 4：在max（x）中：沒有非缺失參數爲max;返回-Inf 然後我試着定義xlim =和ylim =，它不適用於NULL或#s。我真的感謝您的幫助 – boccultist

嗯 - 我曾希望該數據將提取的這種方式。可能需要用str（lmod452）進行挖掘以查看數據的位置。最簡潔的解決方案是使用原始數據（如果可用）。 –

審查中的R線性模型的觀測子集

回答

相關問題