请选择 进入手机版 | 继续访问电脑版
楼主: isodagreen
6639 1

求助:有关R软件多元总体的均值和方差检验程序 [推广有奖]

  • 0关注
  • 0粉丝

小学生

0%

还不是VIP/贵宾

-

威望
0
论坛币
2 个
通用积分
0
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
114 点
帖子
1
精华
0
在线时间
0 小时
注册时间
2008-5-13
最后登录
2022-6-16

isodagreen 发表于 2008-5-13 16:57:00 |显示全部楼层 |坛友微信交流群

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币

小弟想做多元总体的均值和方差检验,可是程序不知如何写!

期待各位大哥大姐的帮助!

谢谢先!

:)

二维码

扫码加我 拉你入群

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

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

关键词:有关R软件 方差检验 r软件 大哥大 如何写 方差 总体 程序 软件 均值

回帖推荐

adobephoebus 发表于2楼  查看完整内容

You could use summary.manova and mauchly.testfor help type ?summary.manova, ?mauchly.test in R.Examples from R help:   ## Example on producing plastic film from Krzanowski (1998, p. 381)     tear <- c(6.5, 6.2, 5.8, 6.5, 6.5, 6.9, 7.2, 6.9, 6.1, 6.3,               6.7, 6.6, 7.2, 7.1, 6.8, 7.1, 7.0 ...

本帖被以下文库推荐

adobephoebus 发表于 2008-5-13 22:06:00 |显示全部楼层 |坛友微信交流群

summary.manova and mauchly.test

You could use summary.manova and mauchly.test

for help type ?summary.manova, ?mauchly.test in R.

Examples from R help:

   ## Example on producing plastic film from Krzanowski (1998, p. 381)
     tear <- c(6.5, 6.2, 5.8, 6.5, 6.5, 6.9, 7.2, 6.9, 6.1, 6.3,
               6.7, 6.6, 7.2, 7.1, 6.8, 7.1, 7.0, 7.2, 7.5, 7.6)
     gloss <- c(9.5, 9.9, 9.6, 9.6, 9.2, 9.1, 10.0, 9.9, 9.5, 9.4,
                9.1, 9.3, 8.3, 8.4, 8.5, 9.2, 8.8, 9.7, 10.1, 9.2)
     opacity <- c(4.4, 6.4, 3.0, 4.1, 0.8, 5.7, 2.0, 3.9, 1.9, 5.7,
                  2.8, 4.1, 3.8, 1.6, 3.4, 8.4, 5.2, 6.9, 2.7, 1.9)
     Y <- cbind(tear, gloss, opacity)
     rate <- factor(gl(2,10), labels=c("Low", "High"))
     additive <- factor(gl(2, 5, len=20), labels=c("Low", "High"))

     fit <- manova(Y ~ rate * additive)
     summary.aov(fit)           # univariate ANOVA tables
     summary(fit, test="Wilks") # ANOVA table of Wilks' lambda


     example(SSD) # Brings in the mlmfit and reacttime objects

     ### traditional test of intrasubj. contrasts
     mauchly.test(mlmfit, X=~1)

     ### tests using intra-subject 3x2 design
     idata <- data.frame(deg=gl(3,1,6, labels=c(0,4,8)),
                         noise=gl(2,3,6, labels=c("A","P")))
     mauchly.test(mlmfit, X = ~ deg + noise, idata = idata)
     mauchly.test(mlmfit, M = ~ deg + noise, X = ~ noise, idata=idata)

已有 1 人评分论坛币 热心指数 收起 理由
crystal8832 + 10 + 1 观点有启发

总评分: 论坛币 + 10  热心指数 + 1   查看全部评分

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-3-29 21:37