AvgDailyGain(SASmixed)
AvgDailyGain()所属R语言包:SASmixed
Average daily weight gain of steers on different diets
的阉牛平均日增重在不同的饮食
译者:生物统计家园网 机器人LoveR
描述----------Description----------
The AvgDailyGain data frame has 32 rows and 6 columns.
AvgDailyGain数据框具有32行和6列。
格式----------Format----------
This data frame contains the following columns:
该数据框包含以下几列:
Id the animal number
ID的动物数量
Block an ordered factor indicating the barn in which the steer was housed.
表示的谷仓中,被安置的督导座有序的因素。
Treatment an ordered factor with levels 0 < 10 < 20 < 30 indicating the amount of medicated feed additive added to the base ration.
处理一个有序的因子与水平0<10<20<30表示药饲料添加剂的量添加到基底口粮。
adg a numeric vector of average daily weight gains over a period of 160 days.
一个为期160天的平均日增重ADG一个数值向量。
InitWt a numeric vector giving the initial weight of the animal
给予动物的初始重量InitWt一个数值向量
Trt the Treatment as a numeric variable
TRT Treatment作为数字变量
源----------Source----------
Littel, R. C., Milliken, G. A., Stroup, W. W., and Wolfinger, R. D. (1996), SAS System for Mixed Models, SAS Institute (Data Set 5.3).
LITTEL,RC,美利肯,GA,斯特普,WW,Wolfinger,RD(1996年),SAS系统的混合模型,SAS研究所(数据集5.3)。
实例----------Examples----------
str(AvgDailyGain)
if (require("lattice", quietly = TRUE, character = TRUE)) {
## plot of adg versus Treatment by Block[#图块助理总干事与治疗]
xyplot(adg ~ Treatment | Block, AvgDailyGain, type = c("g", "p", "r"),
xlab = "Treatment (amount of feed additive)",
ylab = "Average daily weight gain (lb.)", aspect = "xy",
index.cond = function(x, y) coef(lm(y ~ x))[1])
}
if (require("lme4", quietly = TRUE, character = TRUE)) {
options(contrasts = c(unordered = "contr.SAS", ordered = "contr.poly"))
## compare with output 5.1, p. 178[#与输出5.1,对比较。 178]
print(fm1Adg <- lmer(adg ~ InitWt * Treatment - 1 + (1 | Block),
AvgDailyGain))
print(anova(fm1Adg)) # checking significance of terms[检查意义的术语]
print(fm2Adg <- lmer(adg ~ InitWt + Treatment + (1 | Block),
AvgDailyGain))
print(anova(fm2Adg))
print(lmer(adg ~ InitWt + Treatment - 1 + (1 | Block), AvgDailyGain))
}
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|