[Traminer-users] seqplot

Weldon, Mat m.weldon at lancaster.ac.uk
Thu Aug 9 15:21:24 CEST 2012


Sorry, you're right. It doesn't produce an error but it produces a single plot for the entire dataset rather than the group you are interested in. I didn't spot this when I tried it.

Best wishes,

Mat

-----Original Message-----
From: traminer-users-bounces at lists.r-forge.r-project.org [mailto:traminer-users-bounces at lists.r-forge.r-project.org] On Behalf Of Jan Goebel
Sent: 09 August 2012 14:14
To: Users questions
Subject: Re: [Traminer-users] seqplot

I think the idea og Mat is right, but i guess that

seqdplot(mvad.seq[cluster4==4,])

makes more sense and not

 >   	seqdplot(mvad.seq, group=cluster4[cluster4==4])

which should result in an error.

Best,
Jan


On 08/09/2012 03:05 PM, Weldon, Mat wrote:
> Hi Maximilian,
>
> I'm answering based on experience with seqedplot, but it should be the same. If your group argument is just a vector of cluster assignments such as the result of a PAM clustering, for example:
> 	cluster4 <- pam(dd, diss=TRUE, k=4, cluster.only=TRUE) #where "dd" is 
> the dissimilarity matrix
>
> then to plot only one group you should be able to input something like:
>   	seqdplot(mvad.seq, group=cluster4[cluster4==4])
>
> Is that helpful?
>
> Kind regards,
>
> Mat
>
> -----Original Message-----
> From: traminer-users-bounces at lists.r-forge.r-project.org 
> [mailto:traminer-users-bounces at lists.r-forge.r-project.org] On Behalf 
> Of Maximilian Rothfuß
> Sent: 09 August 2012 12:12
> To: traminer-users at lists.r-forge.r-project.org
> Subject: [Traminer-users] seqplot
>
> Hello all together
>
> In my analysis I remarked that there is no possibility within the 'seqdplot´ - function for plotting just one level of the factor given by the `group` argument. Further I couldn´t find any solution for that problem, so I wrote a large function containing some loops, very complicated.
> It works but it's takes so much time...
> Of course, I know that it's possible to get the rigt group by seperating them like:
>
> data(mvad)
> mvad.seq <- seqdef(mvad, 17:86)
>
> ## Using the group argument
> seqdplot(mvad.seq, group=mvad$gcse5eq)
>
> ## Plotting groups separately
> seqdplot(mvad.seq[mvad$gcse5eq=="yes",])
> seqdplot(mvad.seq[mvad$gcse5eq=="no",])
>
> But if the grouping is a result from a cluster analysis, you cannot plot just one certin group.
>
> Did I overlook somthing within the code of seqdplot or TraMineR?
>
> Thank you in anticipation!
>
> Best regrads,
>
> Maximilian
> _______________________________________________
> Traminer-users mailing list
> Traminer-users at lists.r-forge.r-project.org
> https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/traminer-
> users _______________________________________________
> Traminer-users mailing list
> Traminer-users at lists.r-forge.r-project.org
> https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/traminer-
> users
>

--
-----------------------------------------
  Dr. Jan Goebel
  Head of the Division
  Data Operation and Research Data Center

  DIW Berlin
  Socio-Economic Panel Study (SOEP)
  Mohrenstr. 58
  D-10117 Berlin -- Germany --
  phone: +49 30 89789-377
-----------------------------------------

_______________________________________________
Traminer-users mailing list
Traminer-users at lists.r-forge.r-project.org
https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/traminer-users


More information about the Traminer-users mailing list