全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
4958 15
2011-11-13
kew.rar
大小:(4.05 KB)

 马上下载

本附件包括:

  • kew.dat



我真的无力了~!

一道R语言的题。。。各位大侠大哥美女牛人。。。能否帮我解答一下。。

关于Kruskal Wallis test。

这样子的。。
这个test statistic 有 X^{2} 分布, I - 1 的df, 至少3个groups,每个至少4个observation,
数据是关于每个月下雨量。

问题是
(a)用graph表示这个分布 in each month of the year,(命令month.abb可能会用的)
(b)不用命令kruskal.test,用其他方法,来test这个,来看是否有足够证据来说明每个月之间的下雨量有差、
必须有test statistic的值,degree of freedom,5%和1%的critical value,p值。
感谢!!

我已经快不行。。。。。求求你们了!!数据在附件。
二维码

扫码加我 拉你入群

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

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

全部回复
2011-11-13 18:23:46
给你顶一个 求高手
二维码

扫码加我 拉你入群

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

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

2011-11-13 19:36:26
我先来试试第一个问题

我对第一个问题的理解是,画出各年每个月分里面的下雨量的图。
首先,生成一个名为date的日期变量:
Sys.setlocale("LC_TIME", "C") #set the time zone;
kew$date <- as.Date(paste(1, kew$Month, kew$Year), format="%d %b %Y")
为了便于分析,我们把变量里面的月分提取出来,这样可以得到一个数字形式的月分:
kew$month <- substr(as.character(kew$date),6,7)
现在得到如下形式的数据:
Month Year       Rain             date       month
1      Jan 1900  74.400000 1900-01-01    01
2      Feb 1900  80.500000 1900-02-01    02
3      Mar 1900  23.600000 1900-03-01    03
4      Apr 1900  23.600000 1900-04-01    04
5      May 1900  25.100000 1900-05-01    05
6      Jun 1900  53.300000 1900-06-01    06
7      Jul 1900  31.700000 1900-07-01    07
8      Aug 1900  67.300000 1900-08-01    08
9      Sep 1900  26.400000 1900-09-01    09
10     Oct 1900  40.900000 1900-10-01    10
然后,统计每个月里面的平均下雨量(这里用到的是doBy Package里面的summaryBy函数):
library(doBy)
myFun <- Function(x) {c(ave=mean(x), sd=sd(x))}
a <- summaryBy(Rain~month,data=kew, FUN=myFun)
下面就a进行画图:
library(ggplot2)
p <- ggplot(a, aes(x=month,y=Rain.ave))
p+geom_point()
得到如下图所示的图:
1.png

把图修饰修饰得到如下的图:
1.png

再研究研究第二个问题。
二维码

扫码加我 拉你入群

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

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

2011-11-13 23:54:43
jacosis 发表于 2011-11-13 19:36
我先来试试第一个问题

我对第一个问题的理解是,画出各年每个月分里面的下雨量的图。
是。。。。

太感谢你了、/。/。/。。我等你接下来的答案!
二维码

扫码加我 拉你入群

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

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

2011-11-14 00:55:52
jacosis 发表于 2011-11-13 19:36
我先来试试第一个问题

我对第一个问题的理解是,画出各年每个月分里面的下雨量的图。
是不是需要先定义kew是什么。。因为

Error in paste(1, kew$Month, kew$Year) : object "kew" not found
二维码

扫码加我 拉你入群

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

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

2011-11-14 01:04:22
jacosis 发表于 2011-11-13 19:36
我先来试试第一个问题

我对第一个问题的理解是,画出各年每个月分里面的下雨量的图。
大神,您能把所有步骤都发出来吗。。。

恩,关于第二个问,可以用aov (ANOVA)来解答吗?
二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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