立即注册 登录
公卫人 返回首页

ktorresjie的个人空间 https://www.epiman.cn/?505 [收藏] [复制] [RSS]

日志

如何用R实现细面图(Spaghetti plot)

热度 1已有 2721 次阅读2010-7-29 18:13 |个人分类:统计作图|

 

       如何用R实现细面图(Spaghetti Plot

 

     细面图是用于观察研究对象随时间整体变化趋势的一种方法。下面Applied Longitudinal Data Analysis: Modeling Change and Event Occurrence by Judith D. Singer and John B.中的数据为例,用R来实现。

> tolerance<-read.table("tolpp.csv",sep=",", header=T) #读入数据
> head(tolerance, n=10) #观察数据结构
 
 
 
     从数据结构上可以看到:每个对象在5个时间点上受到观察。下面我们以tolerance 为结局变量来作细面图。
> interaction.plot(tolerance$time,tolerance$id, tolerance$tolerance,xlab="time",ylab="Tolerance",col=c(1:10), legend=F) 
 
        显示如下:
                       Cited from:http://www.ats.ucla.edu/stat/R/faq/spagplot.htm

路过

鸡蛋

鲜花

握手

雷人

发表评论 评论 (1 个评论)

回复 ynsxx 2010-9-2 09:00
不错,据我所知,用stata和HLM也可以,相对来说,HLM较为方便,还可以只随机选择10%的样本量,另外可以在其中对性别这样的变量分别显示出来。感觉HLM的功能也不错。当然stata也可以,只是速度慢点。不过,话说回来,用R是一个好的选择,加油。

facelist

您需要登录后才可以评论 登录 | 立即注册

手机版|会员|至尊|接种|公卫人 ( 沪ICP备06060850号-3 )

GMT+8, 2024-4-30 15:22 , Processed in 0.029961 second(s), 7 queries , Gzip On, MemCached On.

Powered by Discuz! X3.4

© 2001-2023 Discuz! Team.

返回顶部