在同一图表上绘制具有两个条件的线性模型图?

时间:2015-04-11 12:08:52

标签: r regression linear

我有一组看起来像这样的数据:

  > 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的新手,很容易混淆,所以任何帮助都会很棒,抱歉这个超长问题!

0 个答案:

没有答案