我有一组看起来像这样的数据:
> head(mydata)
name noise cond rho dpr se type
1 ah pow anim 0.5 1.777938 0.1867230 pla
2 ah pow anim 0.6 2.150651 0.1811027 pla
3 ah pow anim 0.7 1.672588 0.1775456 pla
4 ah pow anim 0.8 2.373414 0.1962687 pla
5 ah pow anim 0.9 2.806913 0.2299617 pla
6 ah pow mirr 0.5 1.707879 0.1645585 plm
我按类型划分子集(总共四个; pla,plm,mama,mamm) 我已经将线性模型应用于它们,如下所示:
rho<-c(.5,.6,.7,.8,.9) ## defines rho
x<-1/sqrt(1-rho)
x
x<-1/sqrt(1-pla$rho)
plafit<-lm(pla$dpr~x)
summary(plafit)
并给出了这个:
Residuals:
Min 1Q Median 3Q Max
-1.42392 -0.23589 0.01057 0.28924 2.06383
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.7819 0.2771 2.822 0.00653 **
x 0.5857 0.1297 4.517 3.14e-05 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 0.6267 on 58 degrees of freedom
Multiple R-squared: 0.2602, Adjusted R-squared: 0.2475
F-statistic: 20.4 on 1 and 58 DF, p-value: 3.135e-05
不同的类型是不同的条件,我想创建一个图表,在同一图表中有两个条件(pla和plm),y轴为dpr,x轴为rho。我很困惑如何将plm和pla放到同一个图表上我已经设法做到这一点看起来像这样:
http://i.stack.imgur.com/KlO8R.png
我是r的新手,很容易混淆,所以任何帮助都会很棒,抱歉这个超长问题!