请教不同自变量水平下的条件分布作图

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

请教古扎拉蒂《计量经济学基础》第5版表2.1的作图如何用R实现

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：自变量 计量经济学基础 经济学基础 计量经济学 经济学基 自变量

图上的蓝字是stata的作图命令

x   y
1   80  55
2   80  60
3   80  65
4   80  70
5   80  75
6   80  65
7  100  65
8  100  70
9  100  74
10 100  80
11 100  85
12 100  88
13 100  77
14 120  79
15 120  84
16 120  90
17 120  94
18 120  98
19 120  89
20 140  80
21 140  93
22 140  95
23 140 103
24 140 108
25 140 113
26 140 115
27 140 101
在r中将数据如上形式输入并命名为data，使用命令plot(data)可作出散点图，但条件均值的连线不知如何做出

sunqinggang 发表于 2013-3-21 14:08
x   y
1   80  55
2   80  60

lm.fit=lm(y~x)
abline(lm.fit)

谢谢楼上
晚上又摸索了下，发现如下方法：
将数据输入改为
> data1
    x1  x2  x3  x4  x5  x6  x7  x8
1   80  80  80  80  80  80  80  80
2  100 100 100 100 100 100 100 100
3  120 120 120 120 120 120 120 120
4  140 140 140 140 140 140 140 140
5  160 160 160 160 160 160 160 160
6  180 180 180 180 180 180 180 180
7  200 200 200 200 200 200 200 200
8  220 220 220 220 220 220 220 220
9  240 240 240 240 240 240 240 240
10 260 260 260 260 260 260 260 260
> data2
    y1  y2  y3  y4  y5  y6  y7   y
1   55  60  65  70  75  NA  NA  65
2   65  70  74  80  85  88  NA  77
3   79  84  90  94  98  NA  NA  89
4   80  93  95 103 108 113 115 101
5  102 107 110 116 118 125  NA 113
6  110 115 120 130 135 140  NA 125
7  120 136 140 144 145  NA  NA 137
8  135 137 140 152 157 160 162 149
9  137 145 155 165 175 189  NA 161
10 150 152 175 178 180 185 191 173
使用命令
> matplot(data1,data2,type="p",pch=15:16)
> lines(data1$x8,data2$y)
可得近似书中的图片效果。

4楼的方法画下来没有条件均值点，今晚又想了个方法，还是有点笨拙
数据：
> da1
    V1  V2
1   80  80
2   80  80
3   80  80
4   80  80
5   80  80
6  100 100
7  100 100
8  100 100
9  100 100
10 100 100
11 100 100
12 120 120
13 120 120
14 120 120
15 120 120
16 120 120

> da2
    V1  V2
1   55  65
2   60  65
3   65  65
4   70  65
5   75  65
6   65  77
7   70  77
8   74  77
9   80  77
10  85  77
11  88  77
12  79  89
13  84  89
14  90  89
15  94  89
16  98  89
17  80 101
18  93 101
19  95 101
20 103 101
21 108 101
22 113 101
23 115 101
24 102 113
25 107 113
26 110 113
27 116 113
28 118 113
29 125 113
30 110 125
31 115 125
32 120 125
33 130 125
34 135 125
35 140 125
36 120 137
37 136 137
38 140 137
39 144 137
40 145 137
41 135 149
42 137 149
43 140 149
44 152 149
45 157 149
46 160 149
47 162 149
命令：
>  matplot(da1,da2,pch=16)
> lines(da1$V1,da2$V2)
得图

excel中列出原始数据（不包括各组均值）
x  y
······
复制后，进入r
da<-read.table("clipboard",header=T)
ada<-aggregate(da,by=list(da$x),FUN=mean)
ada1<-ada[,c(2:3)]
total<-rbind(da,ada1)
plot(total,pch=16)
lines(ada$x,ada$y)

同上的数据来源，
da<-read.table("clipboard",header=T)
ada<-aggregate(da,by=list(da$x),FUN=mean)
plot(da$x,da$y,pch=20)
lines(ada$x,ada$y)
points(ada$x,ada$y,pch=19,col="red")

使用ggplot可以更便捷
数据分为四列：x，y，x（ex），y的条件期望ey
p<-ggplot(da,aes(x,y))+geom_point(size=2,col="blue")
p+geom_point(aes(ex,ey),size=4,col="red")