找回密码
 注册
查看: 425|回复: 0

R语言 tlmec包 UTIdata()函数中文帮助文档(中英文对照)

[复制链接]
发表于 2012-10-1 10:43:59 | 显示全部楼层 |阅读模式
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:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

手机版|小黑屋|生物统计家园 网站价格

GMT+8, 2025-6-18 21:05 , Processed in 0.020390 second(s), 16 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表