首页 > 其他 > 详细

重复测量方差分析教程——to 小卢

时间:2021-01-19 00:44:49      阅读:6      评论:0      收藏:0      [点我收藏+]

重复测量方差分析

读取代码

将代码repeat_anova.R放入E盘主目录

setwd("E:\\")
source("repeat_anova.R")

读取数据:

数据格式如下(第一列为分组,后面为各时间点顺序):

技术分享图片

若数据在此目录下技术分享图片, 则输入

setwd("E:\\我的坚果云\\work")
data = read.csv("data1.csv", header = T, sep=‘,‘)

#	> head(data)
#	   X T1 T2 T3 T4 T5 T6 T7
#	1 G1  4  3  4  4  4  4  3
#	2 G1  4  4  4  3  3  3  3
#	3 G1  4  3  4  4  4  4  3
#	4 G1  4  3  4  4  4  4  3
#	5 G1  3  3  4  4  3  4  3
#	6 G1  4  3  4  3  4  3  3

分析

输入数据信息,n.g为分组个数,n.t为时间点个数,N代表样本量

RES = repeat.anova(data,n.g = 5,n.t = 7,N = 150)
res = RES$res # res为主要结果
write.csv(res,"result1.csv")
RES$effect # effect是重复测量的分组效应,时间效应和交互效应
RES$note # 结果注释

重复测量方差分析教程——to 小卢

原文:https://www.cnblogs.com/biostat-yu/p/14295859.html

(0)
(0)
   
举报
评论 一句话评论(0
© 2014 bubuko.com 版权所有 鲁ICP备09046678号-4
打开技术之扣,分享程序人生!
             

鲁公网安备 37021202000002号