From: Zhongming Yang <zhongmingyang at yahoo.com>
To: "Francisco J. Zagmutt" <gerifalte28 at hotmail.com>,
r-help at stat.math.ethz.ch
Subject: RE: [R] how to draw xyplot figure like figure 4.18 of MASS (4th) ?
Date: Wed, 16 Mar 2005 09:01:07 -0800 (PST)
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