2014-09-05 45 views
3

我还是新的TraMineR;因此,我的问题可能对你们大多数人来说非常简单。 我正在使用我的数据处理一些序列图,并希望用调查权重和名义权重查看结果。我能够将数据导入到R中,并使用数据运行序列图,如,但无法使用权​​重运行它。我已经将我的语法包含在我如何创建状态序列对象(相当天真)和两个序列图示例中。我如何编写使用权重的语法?谢谢TraMineR使用权重

数据文件:DISDATAC; 调查权重:WTSURVY; 标称权重:NORMWEIGHTS

DISDATAC<-read.spss("/Users/tac0432/Desktop/DISDATAC.sav",to.data.frame=TRUE,use.value.labels=TRUE,) 
DISDATAC.shortlab <- c("2yr", "4yr", "Gov","Bus","UeL") 
DISDATAC.lab <- c("2year college or Other School System", " 4year college or Medical Institution", " Government","Business/Industry","Unemployed or Not in the labor force") 
DISDATAC.seq <- seqdef(DISDATAC, 09:12, states = DISDATAC.shortlab, labels = DISDATAC.lab) 
DISDATAC.seq[1:20, 1:4] ; print(DISDATAC.seq[1:20, ], format = "SPS"); alphabet(DISDATAC.seq) 
stlab(DISDATAC.seq) 

Sequence Plot Syntax examples 
seqIplot(DISDATAW.seq, border = NA, group = DISDATAW$NMRMEMG, sortv = "from.start") 
seqIplot(DISDATAW.seq, border = NA, group = DISDATAW$GENDER:DISDATAW$NMRMEMG, sortv = "from.start") 

回答

3

难是没有数据的样本一起玩肯定的,但它看起来像你需要在seqdef功能提供weights列。如果你这样做,默认情况下,绘图功能将考虑权重。

如果您按照seqIplot的帮助,可以进入讨论权重的帮助页面。下面是从帮助的关键部分为seqIplot

如果创建seqdata时已经指定的权重,加权 频率将由seqtab因为默认选项是 weighted=TRUE返回。有关可选参数的完整列表,请参见以下示例,seqtabplot.stslist.freq 手册页以及Mülleret al。,(2008)的关于序列频率图的描述。

所以我们需要看看函数seqdef,它创建了用于绘图的seqdata对象。对于seqdef帮助介绍weights说法:

weights可选数值包含矢量的权重,这是 通过绘制和统计功能适用时入账。

该帮助还有一个示例,说明如何在seqdef中指定权重。