HR(SASmixed)
HR()所属R语言包:SASmixed
Heart rates of patients on different drug treatments
不同的药物治疗患者的心脏率
译者:生物统计家园网 机器人LoveR
描述----------Description----------
The HR data frame has 120 rows and 5 columns of the heart rates of patients under one of three possible drug treatments.
HR数据框有120行5列的患者的心脏率下三种可能的药物治疗。
格式----------Format----------
This data frame contains the following columns:
该数据框包含以下几列:
Patient an ordered factor indicating the patient.
患者一个有序的指示病人的因素。
Drug the drug treatment - a factor with levels a, b and p where p represents the placebo.
药物的药物治疗 - 一个因素与水平a,b和p其中p代表的安慰剂。
baseHR the patient's base heart rate
baseHR病人的基本心脏率
HR the observed heart rate at different times in the experiment
HR在不同的时间在实验中所观察到的心脏速率
Time the time of the observation
时的时间的观察
源----------Source----------
Littel, R. C., Milliken, G. A., Stroup, W. W., and Wolfinger, R. D. (1996), SAS System for Mixed Models, SAS Institute (Data Set 3.5).
LITTEL,RC,美利肯,GA,斯特普,WW,Wolfinger,RD(1996年),SAS系统的混合模型,SAS研究所(数据集3.5)。
实例----------Examples----------
str(HR)
if (require("lattice", quietly = TRUE, character = TRUE)) {
xyplot(HR ~ Time | Patient, HR, type = c("g", "p", "r"), aspect = "xy",
index.cond = function(x, y) coef(lm(y ~ x))[1],
ylab = "Heart rate (beats/min)")
}
if (require("lme4", quietly = TRUE, character = TRUE)) {
options(contrasts = c(unordered = "contr.SAS", ordered = "contr.poly"))
## linear trend in time[#线性趋势]
print(fm1HR <- lmer(HR ~ Time * Drug + baseHR + (Time|Patient), HR))
print(anova(fm1HR))
## Not run: [#不运行:]
fm2HR <- update(fm1HR, weights = varPower(0.5)) # use power-of-mean variance[使用电源的均方差]
summary(fm2HR)
intervals(fm2HR) # variance function does not seem significant[方差函数似乎并没有显著]
anova(fm1HR, fm2HR) # confirm with likelihood ratio[确认似然比]
## End(Not run)[#(不执行)]
print(fm3HR <- lmer(HR ~ Time + Drug + baseHR + (Time|Patient), HR))
print(anova(fm3HR))
## remove Drug term[#清除药物的术语]
print(fm4HR <- lmer(HR ~ Time + baseHR + (Time|Patient), HR))
print(anova(fm4HR))
}
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|