Multilocation(SASmixed)
Multilocation()所属R语言包:SASmixed
A multilocation trial
一个多地点试验
译者:生物统计家园网 机器人LoveR
描述----------Description----------
The Multilocation data frame has 108 rows and 7 columns.
Multilocation数据框有108行7列。
格式----------Format----------
This data frame contains the following columns:
该数据框包含以下几列:
obs a numeric vector
OB的数字矢量
Location an ordered factor with levels B < D < E < I < G < A < C < F < H
位置有序因素与水平BDEIGAC<FH
Block a factor with levels 1 to 3
A座因子,水平13
Trt a factor with levels 1 to 4
TRT一个因素与水平1到4
Adj a numeric vector
调整后的数字矢量
Fe a numeric vector
铁的数字矢量
Grp an ordered factor with levels B/1 < B/2 < B/3 < D/1 < D/2 < D/3 < E/1 < E/2 < E/3 < I/1 < I/2 < I/3 < G/1 < G/2 < G/3 < A/1 < A/2 < A/3 < C/1 < C/2 < C/3 < F/1 < F/2 < F/3 < H/1 < H/2 < H/3
GRP一个有序的因素与水平B/1B/2B/3D/1D/2D/3E/1<E/2E/3I/1I/2I/3G/1<X >G/2G/3A/1A/2A/3C/1C/2 C/3F/1F/2F/3H/1H/2
源----------Source----------
Littel, R. C., Milliken, G. A., Stroup, W. W., and Wolfinger, R. D. (1996), SAS System for Mixed Models, SAS Institute (Data Set 2.8.1).
LITTEL,RC,美利肯,GA,斯特普,WW,Wolfinger,RD(1996年),SAS系统的混合模型,SAS研究所(数据集2.8.1)。
实例----------Examples----------
str(Multilocation)
if (require("lme4", quietly = TRUE, character = TRUE)) {
options(contrasts = c(unordered = "contr.SAS", ordered = "contr.poly"))
### Create a Block %in% Location factor[##创建一个程序段%%位置因素]
Multilocation$Grp <- with(Multilocation, Blockocation)
print(fm1Mult <- lmer(Adj ~ Location * Trt + (1|Grp), Multilocation))
print(anova(fm1Mult))
print(fm2Mult <- lmer(Adj ~ Location + Trt + (1|Grp), Multilocation), corr=FALSE)
print(fm3Mult <- lmer(Adj ~ Location + (1|Grp), Multilocation), corr=FALSE)
print(fm4Mult <- lmer(Adj ~ Trt + (1|Grp), Multilocation))
print(fm5Mult <- lmer(Adj ~ 1 + (1|Grp), Multilocation))
print(anova(fm2Mult))
print(anova(fm1Mult, fm2Mult, fm3Mult, fm4Mult, fm5Mult))
### Treating the location as a random effect[##治疗的位置作为随机效应]
print(fm1MultR <- lmer(Adj ~ Trt + (1|Location/Trt) + (1|Grp), Multilocation))
print(anova(fm1MultR))
print(fm2MultR <- lmer(Adj ~ Trt + (Trt - 1|Location) + (1|Block),
Multilocation))
print(anova(fm1MultR, fm2MultR))
## Not run: [#不运行:]
intervals(fm1MultR)
## End(Not run)[#(不执行)]
}
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|