UTIdata(tlmec)
UTIdata()所属R语言包:tlmec
Data set for Unstructured Treatment Interruption Study
非结构化中断治疗研究中的数据
译者:生物统计家园网 机器人LoveR
描述----------Description----------
Data set from a study of Unstructured Treatment Interruption in HIV-infected adolescents in four institutions in the US. The main outcome is the HIV-1 RNA viral load, which is subject to censoring below the lower limit of detection of the assay (50 copies/mL). The censored observations are indicated by the variable RNAcens
从非结构化中断治疗,在感染艾滋病毒的青少年在美国的四家机构的研究数据集。会议的主要成果是HIV-1 RNA病毒载量,这是低于下限的检测测定(50拷贝/毫升)的审查。审查意见表示的变量RNAcens
用法----------Usage----------
data(UTIdata)
格式----------Format----------
A data frame with 146 observations on the following 5 variables.
有146个观察以下5个变量的数据框。
Patid patient ID
Patid的病人的ID
Days.after.TI days after treatment interruption
Days.after.TI天治疗后中断
Fup follow-up months
Fup后续个月
RNA viral load RNA
RNA病毒载量的RNA
RNAcens censoring indicator for viral load
RNAcens审查病毒载量指标
参考文献----------References----------
实例----------Examples----------
## Not run: [#不运行:]
## load data[#加载数据]
data(UTIdata)
## Sort the data by Patient and visit[#排序病人的数据,并参观]
o <- order(UTIdata$Patid, UTIdata$Fup)
UTIdata <- UTIdata[o,]
## Create censure vector[#创建谴责矢量]
cens = (UTIdata$RNAcens==1)+0
## Generate response vector [#生成的响应向量]
y = log10(UTIdata$RNA)
aa=y[cens==0]
## Create the design matrices[#创建设计矩阵]
x = cbind((UTIdata$Fup==0)+0, (UTIdata$Fup==1)+0, (UTIdata$Fup==3)+0, (UTIdata$Fup==6)+0, (UTIdata$Fup==9)+0, (UTIdata$Fup==12)+0, (UTIdata$Fup==18)+0, (UTIdata$Fup==24)+0)
z = matrix(rep(1, length(y)), ncol=1)
cluster = as.numeric(UTIdata$Patid)
## Create the nj vector[#创建新泽西州矢量。]
nj<-matrix(0,72,1)
for (j in 1:72) {
nj[j]=sum(cluster==j)
}
## Number of individuals [#个人]
m<-dim(nj)[1]
## Call the tlmec with Normal mixed-effects [#调用tlmec与普通混合效应]
out.N <- tlmec(cens,y,x,z,nj,family="Normal",criteria=TRUE)
## Call the tlmec with Student-t mixed-effects [#调用tlmec与学生-T混合效应]
out.T <- tlmec(cens,y,x,z,nj,nu=9,family="t",criteria=TRUE)
## End(Not run)[#(不执行)]
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|