2013-07-25 49 views
3

我正在使用R軟件包TraMineR來計算和分析狀態序列。 我正在嘗試使用seqrplot命令進行繪圖並更改x軸刻度標籤。 的代碼會是這樣,例如:使用seqrplot命令輸入新軸刻度標籤

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)) 

x軸消失,但與推算值厚厚的新中軸線的某處出現在情節的非常右下角。換句話說,它不會取代已刪除的軸。有誰知道我做錯了什麼? 謝謝

回答

4

的問題是,在默認情況下seqrplot(作爲seq?plot家族的其他情節功能)會自動與情節一起顯示的顏色圖例,並使用layout組織情節,並在單個圖形的傳說。如果您想要對軸進行操作,您應該通過設置withlegend = FALSE來禁用自動圖例,然後顯示我的手圖例。例如:

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) 
+0

現在,它的工作。謝謝,我會用其他'seq?plot'命令考慮它。 –