<div dir="ltr">Hi, <div><br></div><div>thanks to both, you are right : </div><div><br></div><div><div>> ncol(trajpro.seq)</div><div>[1] 41</div></div><div style>instead of 21. </div><div style><br></div><div style>In fact i have sequences that are misaligned like stairs year by year, so i do :</div>
<div style><br></div><div style><div>trajpro.seq <- seqdef(immimigpro, 1636:1676, states=labpro.court, labels=labpro.long, xstep=6, left="DEL", right="DEL")</div><div>summary(trajpro.seq)</div><div>
<div>> freq(seqlength(trajpro.seq), total=T)</div><div> n %</div><div>1 1 0.0</div><div>3 1 0.0</div><div>4 1 0.0</div><div>5 1 0.0</div><div>6 2 0.1</div><div>8 1 0.0</div>
<div>9 1 0.0</div><div>10 1 0.0</div><div>11 1 0.0</div><div>12 2 0.1</div><div>13 1 0.0</div><div>16 1 0.0</div><div>17 1 0.0</div><div>18 1 0.0</div><div>
19 2 0.1</div><div>20 7 0.3</div><div>21 2264 98.9</div><div>41 1 0.0</div><div>NA 0 0.0</div><div>Total 2290 100.0</div></div><div><br></div><div style>I guess these 2290-2264 individuals are "real missing values" and do not understand why one remains at 41. </div>
<div><br></div></div><div style><div>immimigpro <- subset(immimigpro, seqlength(trajpro.seq)==21)</div><div><div>freq(seqlength(trajpro.seq), total=T)</div></div><div><br></div></div><div style><div>> trajpro.seq <- seqdef(immimigpro, 1636:1676, states=labpro.court, labels=labpro.long, xstep=6, left="DEL", right="DEL")</div>
<div> [>] found missing values ('NA') in sequence data</div><div> [>] preparing 2264 sequences</div><div> [>] coding void elements with '%' and missing values with '*'</div><div> [>] state coding:</div>
<div> [alphabet] [label] [long label] </div><div> 1 1 Sa Salarié</div><div> 2 2 In Indépendant</div><div> 3 3 Ch Chômage</div><div> 4 4 Et Etudes</div>
<div> 5 5 Fo Foyer</div><div> 6 6 Aut Autre</div><div> 7 7 Var Variable</div><div> [>] 2264 sequences in the data set</div><div> [>] min/max sequence length: 21/21</div>
<div><br></div><div><div>summary(trajpro.seq)</div></div><div><br></div><div style>and then i obtain the graph with 41 periods even if all sequences are all 21-period long. </div><div style><br></div><div style>The solution with <span style="color:rgb(31,73,125);font-family:Calibri,sans-serif;font-size:15px">seqiplot(data.seq[,1:21]) works but it would be better to fix the problem above. Any idea ? </span></div>
<div style><span style="color:rgb(31,73,125);font-family:Calibri,sans-serif;font-size:15px"><br></span></div><div style><span style="color:rgb(31,73,125);font-family:Calibri,sans-serif;font-size:15px">Thanks a lot, best, </span></div>
<div style><span style="color:rgb(31,73,125);font-family:Calibri,sans-serif;font-size:15px"><br></span></div><div style><span style="color:rgb(31,73,125);font-family:Calibri,sans-serif;font-size:15px">Pierre </span></div>
</div><div style><br></div><div style><br></div><div style><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">
2014-03-24 16:14 GMT+01:00 Gilbert Ritschard <span dir="ltr"><<a href="mailto:Gilbert.Ritschard@unige.ch" target="_blank">Gilbert.Ritschard@unige.ch</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div lang="EN-US" link="blue" vlink="purple">
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">Hi Pierre,<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">It would be more useful to get help if you would provide the code you are using to generate your state sequence object.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">My guess is that you created your state sequence object by passing a table with 41 columns (the last 20 ones containing perhaps only NA’s). Assuming your sequence
data are in columns 1 to 21 of a data frame named data, a solution would be to use something like seqdef(data[,1:21]).<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">If your state sequence object is data.seq, you can use seqiplot(data.seq[,1:21]) to plot only the first 21 columns. Fixing the seqdef problem is preferable
however.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">Best.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">Gilbert<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><b><span style="font-size:11.0pt;font-family:"Calibri","sans-serif"">From:</span></b><span style="font-size:11.0pt;font-family:"Calibri","sans-serif"">
<a href="mailto:traminer-users-bounces@lists.r-forge.r-project.org" target="_blank">traminer-users-bounces@lists.r-forge.r-project.org</a> [<a href="mailto:traminer-users-bounces@lists.r-forge.r-project.org" target="_blank">mailto:traminer-users-bounces@lists.r-forge.r-project.org</a>]
<b>On Behalf Of </b>Pierre Blavier<br>
<b>Sent:</b> Sunday, March 23, 2014 18:02<br>
<b>To:</b> Users questions<br>
<b>Subject:</b> Re: [Traminer-users] plots with group option vs separate plots<u></u><u></u></span></p>
<p class="MsoNormal"><u></u> <u></u></p>
<div><div class="">
<p class="MsoNormal">Hi everybody, <u></u><u></u></p>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">here is a plot from the TramineR seqiplot command. It gives the first 10 sequences of my dataset, all my sequences are 21 period-long. Does anyone know how to reduce the length of the x-axis to have a length of the longest sequence, i.e.
21 periods ? It's probably simple but i could not fix it. I have the same problems for other plots <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Thanks, Best, <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Pierre <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div><div>
<p class="MsoNormal"><span style="border:solid windowtext 1.0pt;padding:0in"><img border="0" width="100" height="100" src="cid:~WRD000.jpg" alt="Image removed by sender."></span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div><div><div class="h5">
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><u></u> <u></u></p>
<div>
<p class="MsoNormal">2014-01-15 16:35 GMT+01:00 Gilbert Ritschard <<a href="mailto:Gilbert.Ritschard@unige.ch" target="_blank">Gilbert.Ritschard@unige.ch</a>>:<u></u><u></u></p>
<blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0in 0in 0in 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0in;margin-bottom:5.0pt">
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">Hi Jeremy,</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"> </span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">Thank you for your question. Could you provide a reproducible working example? Indeed you are right,
you should get the same plot in both cases. You do not provide enough information, however, to allow identifying the source of the problem.
</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"> </span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">Your question certainly is of interest for many TraMineR users and future users. I would therefore
suggest you post your question on StackOverflow (see </span><a href="http://mephisto.unige.ch/traminer/contrib.shtml" target="_blank"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif"">http://mephisto.unige.ch/traminer/contrib.shtml</span></a><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">)
using the “traminer” tag which is searchable, unlike this r-forge list.</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"> </span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">Best.</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">Gilbert</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"> </span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"> </span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"> </span><u></u><u></u></p>
<div style="border:none;border-left:solid blue 1.5pt;padding:0in 0in 0in 4.0pt">
<div>
<div style="border:none;border-top:solid #b5c4df 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">From:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">
</span><a href="mailto:traminer-users-bounces@lists.r-forge.r-project.org" target="_blank"><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">traminer-users-bounces@lists.r-forge.r-project.org</span></a><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">
[mailto:</span><a href="mailto:traminer-users-bounces@lists.r-forge.r-project.org" target="_blank"><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">traminer-users-bounces@lists.r-forge.r-project.org</span></a><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">]
<b>On Behalf Of </b>Jeremy Reynolds<br>
<b>Sent:</b> Monday, January 13, 2014 22:09<br>
<b>To:</b> </span><a href="mailto:traminer-users@lists.r-forge.r-project.org" target="_blank"><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">traminer-users@lists.r-forge.r-project.org</span></a><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif""><br>
<b>Subject:</b> [Traminer-users] plots with group option vs separate plots</span><u></u><u></u></p>
</div>
</div>
<p class="MsoNormal"> <u></u><u></u></p>
<div>
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt">Hello,<u></u><u></u></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt">I have been making sequence plots, and I seem to be getting very different results when I use the "group" option of the seqdplot or seqIplot command than when I draw separate plots for
each subgroup. <u></u><u></u></p>
</div>
<p class="MsoNormal">After creating a sequence object and performing optimal matching using PAM, I have chosen a 4 cluster solution. I then create a single plot that shows the distribution across states
in each of the 4 clusters like this:<u></u><u></u></p>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><br>
seqdplot(seq.hc, group = pam5vs$clustering$cluster4, border = NA, title="pam5vs")<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt">If I subset the data and make a separate plot for one of the 4 clusters as in the code below, the N matches the results above (the total N and the N across the states), but I get a very
different impression of how the cases are distributed across the states in the two graphs. Am I doing something wrong? I would be happy to provide more detail if needed.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">Thanks,<br>
<br>
Jeremy<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><br>
cluster4 <- subset(bhps, pam5vs$clustering$cluster4==(6875))<br>
seq.cluster4 <- seqdef(cluster4 [4:21], labels = c("M", "S", "F", "O", "U" ))<br clear="all">
<u></u><u></u></p>
</div>
<div>
<div>
<div>
<p class="MsoNormal">seqdplot(seq.cluster4, border = NA, title="pam5vs cluster 6875")<br>
<br>
-- <u></u><u></u></p>
<div>
<p class="MsoNormal">********************<br>
Dr. Jeremy Reynolds<br>
Associate Professor<br>
Undergraduate Coordinator<br>
Department of Sociology<br>
116 Baldwin Hall<br>
University of Georgia<br>
Athens, GA 30602-1611<br>
Phone: <a href="tel:%28706%29%20583-8072" target="_blank">(706) 583-8072</a><br>
Web: <a href="http://uga.edu/soc/people/faculty/reynolds_jeremy.php" target="_blank">
http://uga.edu/soc/people/faculty/reynolds_jeremy.php</a><br>
Fax: <a href="tel:%28706%29%20542-4320" target="_blank">(706) 542-4320</a><u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<p class="MsoNormal"><br>
_______________________________________________<br>
Traminer-users mailing list<br>
<a href="mailto:Traminer-users@lists.r-forge.r-project.org" target="_blank">Traminer-users@lists.r-forge.r-project.org</a><br>
<a href="https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/traminer-users" target="_blank">https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/traminer-users</a><u></u><u></u></p>
</blockquote>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div></div></div>
</div>
<br>_______________________________________________<br>
Traminer-users mailing list<br>
<a href="mailto:Traminer-users@lists.r-forge.r-project.org">Traminer-users@lists.r-forge.r-project.org</a><br>
<a href="https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/traminer-users" target="_blank">https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/traminer-users</a><br></blockquote></div><br></div>