使用seqrplot命令输入新的轴刻度标签
我使用R包TraMineR来计算和分析状态序列。 我正在尝试使用使用seqrplot命令输入新的轴刻度标签,r,traminer,R,Traminer,我使用R包TraMineR来计算和分析状态序列。 我正在尝试使用seqrplot命令进行绘图,并更改x轴刻度标签。 例如,代码是: library("TraMineR") data("mvad") mvad.alphab <- c("employment", "FE", "HE", "joblessness", "school", "training") mvad.seq <- seqdef(mvad, 17:86, xtstep = 6, alphabet = mvad.alpha
seqrplot
命令进行绘图,并更改x轴刻度标签。
例如,代码是:
library("TraMineR")
data("mvad")
mvad.alphab <- c("employment", "FE", "HE", "joblessness", "school", "training")
mvad.seq <- seqdef(mvad, 17:86, xtstep = 6, alphabet = mvad.alphab)
mvad.dist<-seqdist(mvad.seq, method="DHD")
seqrplot(mvad.seq, dist.matrix = mvad.dist,
criterion = "density", nrep = 1, title = "End CS qualification",
border = NA, axes=FALSE)
axis(1, at = c(1, 6, 12, 18, 24, 30, 36, 42, 48, 54,60,66,70),
labels = c(1, 6, 12, 18, 24, 30, 36, 42, 48, 54,60,66,70))
opar <- par(mfrow=c(1,2))
seqrplot(mvad.seq, dist.matrix = mvad.dist,
criterion = "density", nrep = 1, title = "End CS qualification",
border = NA, axes=FALSE, withlegend=F)
axis(1, at = c(1, 6, 12, 18, 24, 30, 36, 42, 48, 54, 60, 66, 70),
labels = c(1, 6, 12, 18, 24, 30, 36, 42, 48, 54, 60, 66, 70))
seqlegend(mvad.seq)
par(opar)
库(“TraMineR”)
数据(“mvad”)
mvad.alphab问题在于默认情况下seqrplot
(与seq?plot
系列的其他绘图功能一样)会自动将颜色图例与绘图一起显示,并使用layout
在单个图形中组织绘图和图例。如果要对轴执行操作,应通过设置with legend=FALSE
禁用自动图例,然后显示图例my hand(我的手)。例如:
library("TraMineR")
data("mvad")
mvad.alphab <- c("employment", "FE", "HE", "joblessness", "school", "training")
mvad.seq <- seqdef(mvad, 17:86, xtstep = 6, alphabet = mvad.alphab)
mvad.dist<-seqdist(mvad.seq, method="DHD")
seqrplot(mvad.seq, dist.matrix = mvad.dist,
criterion = "density", nrep = 1, title = "End CS qualification",
border = NA, axes=FALSE)
axis(1, at = c(1, 6, 12, 18, 24, 30, 36, 42, 48, 54,60,66,70),
labels = c(1, 6, 12, 18, 24, 30, 36, 42, 48, 54,60,66,70))
opar <- par(mfrow=c(1,2))
seqrplot(mvad.seq, dist.matrix = mvad.dist,
criterion = "density", nrep = 1, title = "End CS qualification",
border = NA, axes=FALSE, withlegend=F)
axis(1, at = c(1, 6, 12, 18, 24, 30, 36, 42, 48, 54, 60, 66, 70),
labels = c(1, 6, 12, 18, 24, 30, 36, 42, 48, 54, 60, 66, 70))
seqlegend(mvad.seq)
par(opar)
opar现在可以工作了。谢谢,我还将使用其他seq?plot
命令将其考虑在内。