2017-07-14 74 views
-3

在R,I有一个线性模型标题lmod452。从这个线性模型中,我想检查来自第25-50号案例的观察结果,以观察哪一个是异化的。如何在选定数量的观测值上使用绘图函数?审查中的R线性模型的观测子集

谢谢!

回答

1

您可以使用[]到子集数据。所以对于假数据

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 
+0

不幸的是我无法按照你的代码来使它工作。 我尝试: > 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

+0

嗯 - 我曾希望该数据将提取的这种方式。可能需要用str(lmod452)进行挖掘以查看数据的位置。最简洁的解决方案是使用原始数据(如果可用)。 –