Re: [R] A better way to do this

From: Dennis Murphy <djmuser_at_gmail.com>
Date: Sun, 22 May 2011 14:02:28 -0700

Hi:

Here are a few options for making this type of plot:

# ggplot2

library(ggplot2)
ggplot(dd, aes(x = Age.at.Sample.Collection, y = MCP1, group = Subject.ID, colour = Subject.ID,

               shape = Group)) +
      geom_point(size = 2.5) + geom_line(size = 0.6) + labs(x = 'Age')

ggplot(dd, aes(x = Age.at.Sample.Collection, y = MCP1, group = Subject.ID, colour = Subject.ID)) +

      geom_point(size = 2.5) + geom_line() + labs(x = 'Age') + facet_wrap( ~ Group, ncol = 1)

ggplot(dd, aes(x = Age.at.Sample.Collection, y = MCP1, group = Subject.ID, colour = Subject.ID)) +

      theme_bw() + geom_point(size = 2.5) + geom_line() +
      labs(x = 'Age') + facet_wrap( ~ Group, nrow = 1)

# lattice
subjectKey <- list(space = 'right',

                   title = 'Subject ID',
                   cex.title = 1.2,
                   text = list(levels(dd$Subject.ID)),
                   points = list(pch = rep(c(16, 17), each = 3), cex =
1, col = 1:6))

xyplot(MCP1 ~ Age.at.Sample.Collection | Group, data = dd, group = Subject.ID,

       type = c('p', 'a', 'g'), cex = 1, pch = rep(c(16, 17), each = 3),
       col = 1:6, xlab = 'Age',
       key = subjectKey)

xyplot(MCP1 ~ Age.at.Sample.Collection, data = dd, group = Subject.ID,
       cex = 1, pch = rep(c(16, 17), each = 3), col = 1:6, type =
c('p', 'a', 'g'), xlab = 'Age',
       key = subjectKey)

HTH,
Dennis

On Sun, May 22, 2011 at 9:46 AM, 1Rnwb <sbpurohit_at_gmail.com> wrote:
> here is the data set, yes i am doing spagghetti plots for each Subject for
> MCP1 with respect to the Age.at.Sample.Collection, with a final of all the
> controls and all the cases
>
>  Subject.ID sample   Group Age.at.Sample.Collection        MCP1
> 19      00173-0   3455 control                11.767282 212.4438625
> 104     00173-0   4992 control                 4.087611 222.6706897
> 148     00173-0   5046 control                 5.103353 257.5997389
> 192     00173-0   5205 control                 7.249828 191.8050229
> 223     00173-0   4991 control                 6.061601 171.2995728
> 235     00173-0   5295 control                 7.843942 240.2288223
> 275     00173-0   5366 control                 9.817932  64.9237708
> 319     00173-0   5473 control                10.704996 199.4368169
> 2102 00173-0   3140 control                15.260780 366.1121484
> 3174 00173-0   4950 control                14.012320 253.4390051
> 3883    00173-0   4103 control                12.974674 260.3490578
> 23      01039-0 901388 control                 3.058179 273.6916770
> 491     01039-0 901068 control                 1.284052 255.5981436
> 735     01039-0 900985 control                 1.968514 450.2616943
> 1199    01039-0 900733 control                 5.086926 161.5011670
> 3948    01039-0 901422 control                 4.041067 365.8560986
> 24      21038-0 901594 control                 8.418891 303.9160091
> 128     21038-0 900089 control                12.440793 250.4846256
> 381     21038-0 901157 control                 6.477754  35.3079462
> 486     21038-0 901739 control                 6.669404  77.1044173
> 582     21038-0 900313 control                10.464065 233.8540385
> 972     21038-0 900527 control                11.438740 194.8173368
> 1256    21038-0 900262 control                14.472279 175.7577419
> 1565    21038-0 901475 control                 2.015058 246.2465683
> 1641    21038-0 901216 control                 4.416153  99.7496760
> 1957    21038-0 901618 control                 6.669404 141.4459310
> 2291    21038-0 900563 control                 3.025325 278.6264989
> 3131    21038-0 900363 control                 5.470225 103.3293132
> 3704    21038-0 900872 control                 9.434633 153.7875157
> 3984    21038-0 900766 control                13.437371 269.4589484
> 2       00041-0   3810    case                17.368925 203.1979716
> 310     00041-0   4923    case                17.648186 168.8229783
> 444     00041-0   3842    case                16.375085 102.6136126
> 1082    00041-0   4078    case                16.813141  89.7548225
> 1459    00041-0   5005    case                14.817248 131.0459201
> 1464    00041-0   4816    case                16.082135 131.1120569
> 1470    00041-0   3495    case                15.323750 118.7034216
> 1677    00041-0   3991    case                 9.782340 169.1765831
> 1845    00041-0   4824    case                 8.520191 122.3348897
> 1928    00041-0   4599    case                16.651608   2.3051960
> 2142    00041-0   3146    case                15.164955 191.6823727
> 2184    00041-0   3220    case                15.526351 188.6667153
> 2283    00041-0   5239    case                18.269678 174.9620952
> 2542    00041-0   3536    case                10.910335 301.6238182
> 2883    00041-0   4787    case                13.527720 159.0956166
> 2921    00041-0   4058    case                17.845311 137.7700442
> 2932    00041-0   5015    case                 7.520876 101.6216590
> 3641    00041-0   4266    case                 5.497604 125.6831543
> 3683    00041-0   4360    case                 4.421629 187.2961108
> 3725    00041-0   4428    case                13.451060 157.4203222
> 3767    00041-0   4506    case                 6.507871 141.1115853
> 3889    00041-0   4455    case                12.147843 213.7489685
> 3       00709-0   5023    case                 2.004106 412.4962003
> 398     00709-0   3180    case                 5.234770   9.6108457
> 639     00709-0   3205    case                 4.974674 262.5951990
> 1058    00709-0   3788    case                 4.577686 281.6220685
> 1700    00709-0   3450    case                 4.156057 131.7354794
> 2050    00709-0   3183    case                 5.878165   0.3707844
> 2483    00709-0   4434    case                 2.162902 339.0410965
> 2789    00709-0   3849    case                 4.405201 405.9685980
> 3281    00709-0   3619    case                 5.541409 255.1460421
> 3842    00709-0   3943    case                 3.515400 624.2636701
> 4011    00709-0   5304    case                 3.312799   2.9206063
> 4       00174-0   4142    case                16.060232 240.4858235
> 13      00174-0   4244    case                 5.615331          NA
> 22      00174-0   4446    case                15.523613 223.3748681
> 57      00174-0   4292    case                 6.209445 168.6679707
> 100     00174-0   4314    case                 6.666666 147.6358300
> 183     00174-0   3689    case                16.224503   0.5689179
> 270     00174-0   4208    case                 3.468856 286.9223725
> 313     00174-0   4223    case                 4.427104 178.0005964
> 471     00174-0   3188    case                11.633127 150.2065531
> 554     00174-0   3212    case                12.377823 217.9799114
> 1304    00174-0   5193    case                14.694045 150.0169015
> 1463    00174-0   4443    case                13.626283 155.9339661
> 1589    00174-0   5294    case                13.990417 167.2262451
> 1884    00174-0   4349    case                15.271731 146.8486102
> 1975    00174-0   3857    case                10.611909  82.6978002
> 2084    00174-0   3174    case                10.924024 217.8106656
> 2196    00174-0   5126    case                16.892539 272.8141490
> 2748    00174-0   4607    case                16.481861 156.5503087
> 2963    00174-0   4861    case                11.967145 294.4989050
> 3294    00174-0   4859    case                12.783025 205.6801413
> 3566    00174-0   4380    case                 9.070499  84.6743348
> 3608    00174-0   4405    case                10.132785  92.1683089
> 3844    00174-0   4252    case                10.324435 305.3678248
> 3923    00174-0   4193    case                11.340177 336.7327758
>>
>
>
> --
> View this message in context: http://r.789695.n4.nabble.com/A-better-way-to-do-this-tp3536576p3542376.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> R-help_at_r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>



R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Sun 22 May 2011 - 21:06:11 GMT

This quarter's messages: by month, or sorted: [ by date ] [ by thread ] [ by subject ] [ by author ]

All messages

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Sun 22 May 2011 - 22:10:09 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive