[R] how to draw xyplot figure like figure 4.18 of MASS (4th) ?

Thu Mar 17 15:10:00 CET 2005

```I think you have restated your question a couple of days ago here
http://tolstoy.newcastle.edu.au/R/help/05/03/14095.html
I will reply to that thread as it is simpler and easier.

The only new thing this thread adds is "where is the data a5 used to
create Figure 4.18 in MASS (4th Ed)" ?

On Wed, 2005-03-16 at 09:01 -0800, Zhongming Yang wrote:
> Dear Francisco and All:
>
> Following is the source code to create figure 4.18 from MASS (4th).
>
> library(MASS)
> library(lattice)
> if(F) { # no data supplied
> xyplot(ratio ~ scant | subject, data = A5,
>       xlab = "scan interval (years)",
>       ylab = "ventricle/brain volume normalized to 1 at start",
>       subscripts = T, ID = A5\$ID,
>       strip = function(factor, ...)
>          strip.default(..., factor.levels = labs, style = 1),
>       layout = c(8, 5, 1),
>       skip = c(rep(FALSE, 37), rep(TRUE, 1), rep(FALSE, 1)),
>       panel = function(x, y, subscripts, ID) {
>           panel.xyplot(x, y, type = "b", cex = 0.5)
>           which <- unique(ID[subscripts])
>           panel.xyplot(c(0, 1.5), pr3[names(pr3) == which],
>                        type = "l", lty = 3)
>           if(which == 303 || which == 341) points(1.4, 1.3)
>       })
> }
>
> But there is no data set available for that, I can't figure out many stuff. So I provide my data set, hope you can help me.
>
> In my data set, there are 3 group patients, and 5 patients in each group. The mds are the repeat measurements on month. I want draw a xyplot with 3 lines, each line for a group.
>
>
> Many thanks
>
> Zhongming Yang
>
>
>
>
> "Francisco J. Zagmutt" <gerifalte28 at hotmail.com> wrote:
> Dear Zhongming,
>
> By asking for the figure in the book you are restricting you question to
> only the people that has the 4th edition. I would love to help you but
> unfortunatelly I have the 3rd edition of MASS and there is no figure 4.18
> since chapter 4 is "Programming in S". Please give us an idea of what you
> want to do and we might be able to help you without looking at the book.
>
> Cheers
>
> Francisco
>
> >From: Zhongming Yang
> >To: r-help at stat.math.ethz.ch
> >Subject: [R] how to draw xyplot figure like figure 4.18 of MASS (4th) ?
> >Date: Wed, 16 Mar 2005 08:07:21 -0800 (PST)
> >
> >Dear All:
> >
> >Could you please tell me how I can draw figure formatted like figure 4.18
> >of MASS (4th) with the attached data set?
> >
> >Thanks
> >
> >Zhongming Yang
> >
> >
> >
> >---------------------------------
> >
> > md idno month group
> > -0.090 521 0.000000 NN
> > -1.330 521 12.460274 NN
> > -0.570 521 14.991781 NN
> > -2.130 559 0.000000 NN
> > -0.920 559 3.978082 NN
> > -1.730 559 6.016438 NN
> > -2.390 559 9.665753 NN
> > -1.300 559 12.460274 NN
> > -2.170 559 15.682192 NN
> > -2.620 559 17.950685 NN
> > -1.830 559 21.664406 NN
> > -3.090 559 25.533258 NN
> > -2.760 559 30.811947 NN
> > -1.240 559 33.830137 NN
> > -0.650 559 36.394521 NN
> > -1.200 559 38.893151 NN
> > -1.320 559 42.147945 NN
> > -2.000 559 45.600000 NN
> > -2.150 559 48.328767 NN
> > -1.430 559 52.569863 NN
> > -2.650 559 56.120548 NN
> > -0.870 559 60.032877 NN
> > -1.920 559 63.550685 NN
> > -2.320 559 65.786301 NN
> > -2.640 559 68.547945 NN
> > -2.030 559 71.533258 NN
> > -1.900 559 75.500472 NN
> > -1.420 559 78.254570 NN
> > -1.330 559 81.863014 NN
> > -2.120 559 84.657534 NN
> > -2.400 559 86.926027 NN
> > -2.410 559 90.147945 NN
> > -3.370 559 93.369863 NN
> > -0.720 559 96.624658 NN
> > -2.075 598 0.000000 NN
> > -2.260 598 3.484932 NN
> > -2.150 598 5.917808 NN
> > -1.420 598 9.828879 NN
> > -1.490 598 12.353470 NN
> > -3.040 598 16.025601 NN
> > -1.960 598 18.845273 NN
> > -2.780 598 22.224658 NN
> > -3.310 598 24.591781 NN
> > -4.420 598 27.517808 NN
> > -1.890 598 30.312329 NN
> > -3.000 598 33.304110 NN
> > -2.260 598 36.065753 NN
> > -0.730 598 40.010959 NN
> > -2.710 598 43.232877 NN
> > -2.430 598 45.994521 NN
> > -1.050 598 49.183562 NN
> > -2.190 598 52.438356 NN
> > -0.580 598 54.969863 NN
> > -2.050 598 57.501011 NN
> > -1.480 598 60.681338 NN
> > -0.960 598 63.664945 NN
> > -2.090 598 66.419043 NN
> > -3.180 598 69.435616 NN
> > -1.230 598 72.394521 NN
> > -2.170 598 74.761644 NN
> > -2.400 598 78.378082 NN
> > -3.630 598 81.698630 NN
> > -2.100 598 84.920548 NN
> > -3.940 598 87.156164 NN
> > -3.780 598 90.378082 NN
> > -5.690 598 95.967123 NN
> > -13.525 622 0.000000 NN
> > -12.670 622 3.680934 NN
> > -13.500 622 7.123557 NN
> > -13.700 622 9.910442 NN
> > -11.260 622 12.926836 NN
> > -12.100 622 15.846575 NN
> > -13.400 622 18.410959 NN
> > -12.330 622 21.435616 NN
> > -12.940 622 24.197260 NN
> > -12.540 622 27.156164 NN
> > -12.550 622 29.917808 NN
> > -11.900 622 33.830137 NN
> > -11.640 622 36.394521 NN
> > -12.100 622 38.893151 NN
> > -11.010 622 41.457534 NN
> > -12.050 622 46.684932 NN
> > -12.650 622 50.136986 NN
> > -11.620 622 54.303885 NN
> > -12.150 622 58.303885 NN
> > -12.990 622 60.959623 NN
> > -12.590 622 63.912329 NN
> > -13.570 622 66.936986 NN
> > -12.480 622 69.336986 NN
> > -12.800 622 72.263014 NN
> > -12.640 622 75.452055 NN
> > -13.830 622 81.073973 NN
> > -13.250 622 84.295890 NN
> > -13.700 622 90.213699 NN
> > -13.250 622 96.197260 NN
> > -2.190 801 0.000000 NN
> > -1.250 801 2.663014 NN
> > -1.550 801 5.589041 NN
> > -1.570 801 8.575612 NN
> > -0.880 801 11.755940 NN
> > -0.230 801 15.428071 NN
> > -0.770 801 17.985448 NN
> > -1.050 801 20.975342 NN
> > -2.260 801 24.394521 NN
> > -2.760 801 27.715068 NN
> > -1.550 801 29.950685 NN
> > -2.330 801 34.093151 NN
> > -1.390 801 36.394521 NN
> > -2.090 801 39.353425 NN
> > -3.050 801 41.917808 NN
> > -1.930 801 44.646575 NN
> > -4.040 801 48.328767 NN
> > -3.260 801 51.780822 NN
> > -2.810 801 54.115068 NN
> > -0.660 801 57.395284 NN
> > -1.800 801 62.903481 NN
> > -1.540 801 65.755940 NN
> > -0.880 801 69.567123 NN
> > -2.000 801 73.117808 NN
> > -1.610 801 76.142466 NN
> > -4.360 801 79.463014 NN
> > -3.200 801 80.712329 NN
> > -3.080 801 83.013699 NN
> > -2.350 801 88.306849 NN
> > -3.170 801 89.917808 NN
> > -3.570 801 95.342466 NN
> > -3.300 525 0.000000 TP
> > -2.500 525 2.991781 TP
> > -5.480 525 5.917808 TP
> > -1.140 525 8.712329 TP
> > -2.260 525 11.967123 TP
> > -0.670 525 14.926027 TP
> > -0.630 525 17.950685 TP
> > -0.600 525 21.336986 TP
> > -1.870 525 24.427397 TP
> > -1.790 525 26.860274 TP
> > -3.640 525 30.771345 TP
> > -2.860 525 33.328722 TP
> > -2.140 525 37.197575 TP
> > -1.560 525 39.787739 TP
> > -0.360 525 42.542466 TP
> > -1.900 525 45.304110 TP
> > -1.190 525 48.032877 TP
> > -3.090 525 51.024658 TP
> > -2.620 525 53.786301 TP
> > -1.360 525 58.126027 TP
> > -1.130 525 60.953425 TP
> > -1.390 525 64.142466 TP
> > -0.270 525 66.904110 TP
> > -1.640 525 69.665753 TP
> > -1.390 525 72.887671 TP
> > -2.550 525 75.221918 TP
> > -2.540 525 78.180822 TP
> > -1.640 525 81.623804 TP
> > -2.410 525 84.574624 TP
> > -2.530 525 88.017247 TP
> > -1.980 525 91.463014 TP
> > -3.060 525 93.895890 TP
> > -1.920 525 96.131507 TP
> > -1.650 525 99.320548 TP
> > -2.625 527 0.000000 TP
> > -2.370 527 2.991781 TP
> > 0.050 527 6.180822 TP
> > -0.270 527 9.041096 TP
> > 1.060 527 11.736986 TP
> > -3.030 527 15.221918 TP
> > -1.000 527 18.180822 TP
> > -1.270 527 21.895890 TP
> > -1.910 527 24.361644 TP
> > -2.350 527 27.649315 TP
> > -3.720 527 30.142825 TP
> > -3.310 527 33.355940 TP
> > -2.280 527 35.847743 TP
> > -1.860 527 39.290366 TP
> > -5.050 527 42.345205 TP
> > -3.160 527 45.106849 TP
> > -2.270 527 47.802740 TP
> > -1.850 527 50.630137 TP
> > -4.010 527 53.621918 TP
> > -4.520 527 56.383562 TP
> > -3.820 527 60.000000 TP
> > -6.440 527 65.293151 TP
> > -4.600 527 68.350685 TP
> > -5.280 527 71.802740 TP
> > -3.920 527 75.320548 TP
> > -5.080 527 78.175612 TP
> > -8.550 527 80.831350 TP
> > -4.900 527 84.601841 TP
> > -11.620 527 87.585448 TP
> > -6.090 527 91.956164 TP
> > -7.120 527 95.408219 TP
> > -6.365 570 0.000000 TP
> > -6.100 570 2.761644 TP
> > -4.950 570 5.983562 TP
> > -5.530 570 9.205479 TP
> > -3.640 570 11.967123 TP
> > -3.330 570 15.879452 TP
> > -4.950 570 18.310442 TP
> > -5.250 570 21.392410 TP
> > -5.540 570 24.376016 TP
> > -4.300 570 27.359623 TP
> > -4.530 570 29.884932 TP
> > -5.970 570 32.646575 TP
> > -4.950 570 37.019178 TP
> > -5.090 570 39.320548 TP
> > -4.000 570 42.312329 TP
> > -6.160 570 45.073973 TP
> > -6.490 570 49.216438 TP
> > -4.890 570 52.668493 TP
> > -6.620 570 56.350685 TP
> > -7.010 570 61.183562 TP
> > -7.290 570 63.978082 TP
> > -6.440 570 66.736672 TP
> > -7.950 570 69.917000 TP
> > -7.570 570 73.130115 TP
> > -8.150 570 76.113721 TP
> > -8.090 570 78.805479 TP
> > -9.670 570 82.553425 TP
> > -9.480 570 86.169863 TP
> > -9.120 570 90.147945 TP
> > -11.330 570 97.019178 TP
> > -3.750 615 0.000000 TP
> > -3.510 615 3.287671 TP
> > -4.160 615 5.848641 TP
> > -5.170 615 9.651920 TP
> > -4.600 615 12.438805 TP
> > -5.940 615 15.619133 TP
> > -6.270 615 18.608219 TP
> > -4.690 615 21.632877 TP
> > -7.820 615 24.361644 TP
> > -5.940 615 27.189041 TP
> > -5.680 615 30.180822 TP
> > -6.750 615 32.712329 TP
> > -5.130 615 36.854795 TP
> > -6.630 615 39.386301 TP
> > -7.450 615 41.917808 TP
> > -7.290 615 45.369863 TP
> > -9.590 615 48.361644 TP
> > -9.700 615 51.353425 TP
> > -9.570 615 54.799461 TP
> > -9.840 615 57.783068 TP
> > -9.800 615 60.766674 TP
> > -11.000 615 63.291264 TP
> > -11.360 615 66.805479 TP
> > -10.860 615 69.731507 TP
> > -10.610 615 72.493151 TP
> > -12.180 615 78.706849 TP
> > -12.570 615 84.690411 TP
> > -13.280 615 90.213699 TP
> > -15.720 615 96.197260 TP
> > 0.375 818 0.000000 TP
> > 1.640 818 2.491803 TP
> > -5.190 818 5.967213 TP
> > -3.080 818 9.674018 TP
> > -4.020 818 12.402785 TP
> > -6.380 818 15.854839 TP
> > -4.570 818 17.926072 TP
> > -5.380 818 21.608264 TP
> > -7.210 818 25.060319 TP
> > -8.560 818 27.854839 TP
> > -7.300 818 30.156209 TP
> > -7.590 818 33.312374 TP
> > -7.050 818 36.600045 TP
> > -7.350 818 38.046620 TP
> > -9.040 818 41.893196 TP
> > -8.340 818 45.573770 TP
> > -6.890 818 48.131148 TP
> > -7.190 818 54.098361 TP
> > -7.690 818 57.312374 TP
> > -9.340 818 60.468538 TP
> > -9.580 818 62.572648 TP
> > -10.180 818 65.827442 TP
> > -9.000 818 69.443881 TP
> > -10.760 818 71.547990 TP
> > -10.580 818 74.112374 TP
> > -10.500 818 78.221963 TP
> > -10.460 818 83.252100 TP
> > -9.980 818 89.071278 TP
> > -4.540 535 0.000000 SN
> > -6.050 535 3.221918 SN
> > -4.040 535 5.983562 SN
> > -2.940 535 8.843836 SN
> > -4.360 535 12.197260 SN
> > -3.760 535 14.958904 SN
> > -2.780 535 18.641096 SN
> > -6.140 535 21.402740 SN
> > -6.070 535 24.164384 SN
> > -6.530 535 27.450528 SN
> > -7.320 535 31.057085 SN
> > -6.270 535 33.778397 SN
> > -6.480 535 36.991511 SN
> > -7.440 535 39.747945 SN
> > -9.400 535 42.509589 SN
> > -6.640 535 45.534247 SN
> > -6.170 535 48.098630 SN
> > -5.090 535 51.090411 SN
> > -4.660 535 53.786301 SN
> > -5.420 535 58.158904 SN
> > -8.420 535 60.920548 SN
> > -4.830 535 63.682192 SN
> > -5.210 535 66.673973 SN
> > -5.030 535 69.961644 SN
> > -4.720 535 72.460274 SN
> > -5.570 535 75.253806 SN
> > -4.240 535 79.122659 SN
> > -5.620 535 82.302987 SN
> > -5.180 535 85.286593 SN
> > -4.560 535 88.076712 SN
> > -3.730 535 91.265753 SN
> > -4.490 535 95.473973 SN
> > -5.010 535 99.221918 SN
> > -3.645 541 0.000000 SN
> > -2.130 541 3.024658 SN
> > -2.240 541 6.016438 SN
> > -4.490 541 9.073973 SN
> > -3.310 541 12.000000 SN
> > -3.320 541 16.339726 SN
> > -4.130 541 18.246575 SN
> > -3.640 541 21.271233 SN
> > -3.510 541 24.624658 SN
> > -5.260 541 27.612666 SN
> > -3.340 541 30.825780 SN
> > -3.990 541 33.678239 SN
> > -5.390 541 37.022502 SN
> > -9.960 541 39.616438 SN
> > -3.740 541 43.002740 SN
> > -4.550 541 45.600000 SN
> > -3.370 541 48.295890 SN
> > -4.810 541 50.893151 SN
> > -3.130 541 53.654795 SN
> > -3.770 541 57.632877 SN
> > -3.820 541 60.558904 SN
> > -5.160 541 63.090411 SN
> > -3.530 541 66.016438 SN
> > -4.690 567 0.000000 SN
> > -4.030 567 2.991781 SN
> > -2.080 567 6.049315 SN
> > -2.430 567 9.698630 SN
> > -3.500 567 12.723288 SN
> > -2.690 567 15.189041 SN
> > -3.260 567 18.476712 SN
> > -4.890 567 21.460768 SN
> > -3.940 567 24.673883 SN
> > -4.150 567 28.313227 SN
> > -2.570 567 31.068493 SN
> > -3.290 567 33.402740 SN
> > -2.650 567 36.394521 SN
> > -1.220 567 39.353425 SN
> > -2.230 567 45.764384 SN
> > -2.780 567 48.526027 SN
> > -3.230 567 51.780822 SN
> > -3.470 567 54.542466 SN
> > -2.440 567 57.041096 SN
> > -4.170 567 61.249315 SN
> > -3.440 567 63.715068 SN
> > -1.880 567 66.279452 SN
> > -2.150 567 69.231260 SN
> > -2.150 567 72.673883 SN
> > -7.010 567 89.720548 SN
> > -4.640 567 94.520548 SN
> > -5.600 567 96.723288 SN
> > -2.375 572 0.000000 SN
> > -2.920 572 3.452055 SN
> > -1.310 572 6.246575 SN
> > -1.890 572 10.158904 SN
> > -2.450 572 12.558904 SN
> > -2.590 572 15.189041 SN
> > -2.040 572 18.308826 SN
> > -1.570 572 21.685875 SN
> > -1.930 572 24.472760 SN
> > -1.300 572 27.817022 SN
> > -2.820 572 30.641096 SN
> > -1.240 572 33.600000 SN
> > -0.720 572 36.394521 SN
> > -0.440 572 39.156164 SN
> > -0.750 572 41.917808 SN
> > -1.310 572 45.567123 SN
> > -1.060 572 48.756164 SN
> > -0.670 572 51.517808 SN
> > -1.110 572 54.279452 SN
> > -0.870 572 57.501370 SN
> > -2.400 572 60.493151 SN
> > -1.790 572 63.484932 SN
> > -1.150 572 66.308826 SN
> > -0.550 572 68.800629 SN
> > -1.830 572 72.112104 SN
> > -1.830 572 75.030137 SN
> > -1.640 572 78.772603 SN
> > -2.340 572 81.468493 SN
> > -2.490 572 86.136986 SN
> > -1.840 572 90.147945 SN
> > -1.760 572 93.369863 SN
> > -0.970 572 96.657534 SN
> > -0.960 583 0.000000 SN
> > -1.150 583 3.189041 SN
> > -2.060 583 7.364384 SN
> > -1.710 583 9.402740 SN
> > -0.870 583 12.427397 SN
> > -2.100 583 15.185179 SN
> > -1.290 583 18.660588 SN
> > -1.350 583 21.185179 SN
> > -0.410 583 24.824523 SN
> > -2.750 583 27.550685 SN
> > -1.860 583 30.345205 SN
> > -1.750 583 33.336986 SN
> > 0.360 583 36.098630 SN
> > -1.200 583 38.893151 SN
> > -0.900 583 42.575342 SN
> > -1.350 583 45.501370 SN
> > -0.660 583 48.526027 SN
> > -0.950 583 51.254795 SN
> > -1.690 583 54.279452 SN
> > -0.410 583 57.073973 SN
> > -1.090 583 60.295890 SN
> > -0.590 583 63.021244 SN
> > -0.800 583 66.004851 SN
> > -0.700 583 69.250752 SN
> > -0.990 583 72.627801 SN
> > -0.240 583 75.879452 SN
> > -0.180 583 78.673973 SN
> > -1.000 583 81.402740 SN
> > -1.360 583 84.164384 SN
> > -1.240 583 86.958904 SN
> > -1.500 583 90.246575 SN
> > -1.410 583 93.468493 SN
> > -4.260 583 96.230137 SN
> >______________________________________________
> >R-help at stat.math.ethz.ch mailing list
> >https://stat.ethz.ch/mailman/listinfo/r-help
> >http://www.R-project.org/posting-guide.html
>
>
>
>
>
> ---------------------------------
>
>
> 	[[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help