全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
5438 9
2013-04-11
悬赏 300 个论坛币 已解决
俺是R语言白丁。请帮忙解释下面一段R程序,比如:哪些变量,关系,每个命令和参数的含义,语句的含义和结果等等。越详细越好。。对解释最佳者,给予300论坛币略表心意。
另外,附上相关数据供参考。。
------------------------------------------------------------------------------
gdp <- read.csv("GDP P chain volume flat.csv",
na.strings="..", check.names=TRUE)
# Create shades of colors between blue and red for labels
shades <- colorRampPalette(c("blue", "red"))(nrow(gdp))
# Exploration of wholesale and retail trade - extract trend series
wholesale.d <- decompose(ts(gdp$Wholesale.Trade, frequency=4, start=c(1987,2)),
type="multiplicative")
wholesale <- wholesale.d$trend
retail.d <- decompose(ts(gdp$Retail.Trade, frequency=4, start=c(1987,2)),
type="multiplicative")
retail <- retail.d$trend
# Exploratory plot 1
plot(wholesale,retail, bty="l")
# Exploratory plot 2
plot(cbind(wholesale,retail), bty="l")
# Main plot
library(MASS)
eqscplot(diff(wholesale), diff(retail), bty="l", cex=.5, col=c("grey60", shades),
xlab="Growth since last quarter in wholesale contribution to GDP, $m",
ylab="Growth since last quarter in retail contribution to GDP, $m")
grid(col="grey90", lty=1)
# Draw diagonal line showing when the two variables are equal
abline(0,1)
# Draw line showing regression
abline(lm(diff(retail) ~ diff(wholesale), na.action="na.omit"), lty=2)



GDP P chain volume flat.rar

大小:7.66 KB

 马上下载

本附件包括:

  • GDP P chain volume flat.csv

最佳答案

求证1加1 查看完整内容

一句一句来吧,第一句read.csv,用来读取数据csv格式的,如果是其他格式的,那么点后面的内容也要变,如读txt就用read.table 第二句,是用来调色,colorRampPalette是一个调色函数 第三句,是时间序列的分解,l即趋势项、季节项、随机项,用滑动平均的方法,并把结果复制到wholesale.d(是个列表)中 第四句,读取列表中的trend(趋势项)部分 后面的语句都是用来作图的,解释一下参数吧。 第五句wholesale和retail分别是画图 ...
二维码

扫码加我 拉你入群

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

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

全部回复
2013-4-11 10:03:40
一句一句来吧,第一句read.csv,用来读取数据csv格式的,如果是其他格式的,那么点后面的内容也要变,如读txt就用read.table
第二句,是用来调色,colorRampPalette是一个调色函数
第三句,是时间序列的分解,l即趋势项、季节项、随机项,用滑动平均的方法,并把结果复制到wholesale.d(是个列表)中
第四句,读取列表中的trend(趋势项)部分
后面的语句都是用来作图的,解释一下参数吧。
第五句wholesale和retail分别是画图的x,y,bty是用于设置图形边框的样式,可选的参数有o, l, 7, c, u,n,这里l其实就是最正常的那种。lz可以自己尝试一下,印象会更深刻。
第六句就不用多说了吧,同理可得,cbind解释一下,就是矩阵的合并,cbind是按列合并,即横向合并,对应的rbind是按行合并。
第七句,加载MASS包,因为eqscplot函数在MASS包中,lz要确认已安装该包
第八句,eqscplot也是用来画图的,一般用来画散点图吧,里面的参数大多继承了plot的,cex指放大的倍数,col设置颜色,用到了之前所调的颜色,xlab,ylab用于给坐标轴命名。
第九句,grid用来添加网格,lty用于设置线条样式,等于1就是指画实线,其余选项:0 不画线,2 虚线,3 点线,4点划线,5 长划线,6 点长划线
最后两句,abline是一个低级绘图函数,用于画线,前提是之前必须调用plot等高级绘图函数。代码注释里也已经说了,用来画回归曲线的。
解释完毕,希望lz有用。
二维码

扫码加我 拉你入群

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

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

2013-4-11 16:29:33
求证1加1 发表于 2013-4-11 10:31
一句一句来吧,第一句read.csv,用来读取数据csv格式的,如果是其他格式的,那么点后面的内容也要变,如读t ...
非常感谢回复,很专业,佩服。
可俺还不是太明白,烦请进一步说明下:
1)详细些说明下“时间序列的分解”是如何做的?你说用“滑动平均”是不是说“moving average”?它是把变量wholesale.trade 做按季度的“滑动平均”了吗?时间序列里“gdp$Wholesale.Trade”是怎么回事(gdp 和wholesale.trade 分别是变量)?存到wholesale.d 里的trend 变量是怎样得到的?或者说,存到“wholesale”里的变量是什么?同样,对retail 也是这个问题。
2)"cbind(wholesale,retail)"把wholesale 和retail按列进行合并的含义,目的或表现是啥?
3)“diagonal line showing when the two variables are equal”是说对角线焦点两个变量相等吗?wholesale 和 retail 是x,y,的话,坐标轴上的标尺表示什么?交点说明什么?
4)最后的回归是什么对什么回归?“abline(lm(diff(retail) ~ diff(wholesale), na.action="na.omit"), lty=2)
”这句话能否详细解释下??
5)另外,将画出的图附上,如方便麻烦就图解释下。 比如 2008, 2002 等线段说明什么?

望不吝赐教。
二维码

扫码加我 拉你入群

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

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

2013-4-11 16:30:31
图片在这里。。

附件列表
11-04-2013 5-51-59 PM.jpg

原图尺寸 138.17 KB

11-04-2013 5-51-59 PM.jpg

二维码

扫码加我 拉你入群

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

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

2013-4-11 20:44:31
好的。
第一个问题,这个问题其实跟R无关,属于时间序列的分解,lz理解基本正确。
第二个问题,其实那个plot意义不大,注释里也说了,只是大概看下效果,而且后面调用了MASS包里高级绘图函数,所以并没有在你那张图上显示出来
第三个问题,那条命令只是为了画一条y=x的线而已,即那条实线
第四个问题,左边是因变量右边是自变量,画一下回归线,lz可以看看lm函数的help
最后一个问题,你这个数据是针对时间序列的,那么必然会有时间标度,所以那个就是代表时间
回答完毕,希望对lz有用。
二维码

扫码加我 拉你入群

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

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

2013-4-12 09:35:17
求证1加1 发表于 2013-4-11 20:44
好的。
第一个问题,这个问题其实跟R无关,属于时间序列的分解,lz理解基本正确。
第二个问题,其实那个p ...
再次感谢及时而又专业的回复。

不过还是问题5 不太明白。你是说X,Y轴的标度是时间吗?"-20" 表示什么?2008,2002等线段表示啥意思?虚线和实线对角线的交点表示啥?

望再次指点。


二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群