Dear Nick,
I attach the .pvs file from which I note that the
predicted splines start at line 26 of the file.
So the Splus command needs to be
asr.pvsplot("c:\\docs\\dracup_m\\96wh14\\test.pvs",26,5,12)
^^
Arthur
> Date: Thu, 29 Oct 1998 09:11:31 +0800
> X-Sender: ngalwey@cyllene.uwa.edu.au
> Mime-Version: 1.0
> To: asreml@ram.chiswick.anprod.csiro.au
> From: "N.W. Galwey" <ngalwey@cyllene.uwa.edu.au>
> Subject: Problem with prediction
>
> I am trying to follow the example on p. 88 of the ASREML manual(version of 2
> October 1998).
>
> I have created a file with five levels for group and 12 unique values in X,
> presented in test.dat below. I have analysed this with the command file
> presented in test.as below. I have then produced predicted values using
> the file presented in test.pin below. I have then tried to produce a plot
> of splines in Splus using the command
>
> asr.pvsplot("c:\\docs\\dracup_m\\96wh14\\test.pvs",1,5,12)
>
> This produces some output - i.e. it finds the file okay - but the lines
> zigzag all over the place. What am I doing wrong?
>
> Nick Galwey
> ____________________
>
> test.dat
>
> group X y
> 1 12 5.11
> 1 13 5.45
> 1 15 6.08
> 1 16 6.32
> 1 18 6.97
> 1 19 7.29
> 1 21 7.87
> 1 24 8.78
> 1 25 9.08
> 1 26 9.40
> 1 27 9.61
> 1 29 10.24
> 2 12 10.24
> 2 13 10.85
> 2 15 12.03
> 2 16 12.66
> 2 18 13.90
> 2 19 14.50
> 2 21 15.61
> 2 24 17.44
> 2 25 18.07
> 2 26 18.65
> 2 27 19.29
> 2 29 20.41
> 3 12 15.39
> 3 13 16.26
> 3 15 18.09
> 3 16 18.95
> 3 18 20.76
> 3 19 21.62
> 3 21 23.47
> 3 24 26.19
> 3 25 27.04
> 3 26 27.92
> 3 27 28.89
> 3 29 30.68
> 4 12 20.42
> 4 13 21.70
> 4 15 24.09
> 4 16 25.21
> 4 18 27.68
> 4 19 28.86
> 4 21 31.30
> 4 24 34.86
> 4 25 36.01
> 4 26 37.21
> 4 27 38.48
> 4 29 40.89
> 5 12 25.53
> 5 13 27.09
> 5 15 30.09
> 5 16 31.60
> 5 18 34.56
> 5 19 36.00
> 5 21 39.02
> 5 24 43.60
> 5 25 45.05
> 5 26 46.56
> 5 27 48.01
> 5 29 51.07
>
> test.as
>
> Example - Manual page 88
> group 5
> X
> y
>
> C:\DOCS\Dracup_M\96WH14\test.dat !skip 1
>
> y ~ mu group group.X !r spl(X) group.spl(X)
>
>
> test.pin
>
> # A is the vector of the 12 knot points
> # B is the vector of the average predicted spline
> C = C+B # Add average curve to Group 1 specific curve
> D = D+B # Add average curve to Group 2 specific curve
> E = E+B # Add average curve to Group 3 specific curve
> F = F+B # Add average curve to Group 4 specific curve
> G = G+B # Add average curve to Group 5 specific curve
> GROUP B6 [ 1:5 ] A6 [6:10 ] 11 # Group means at knot point 6
> GROUP1 C 1 A 6 11 # Spline for group 1
> GROUP2 D 2 A 7 11 # Spline for group 2
> GROUP3 E 3 A 8 11 # Spline for group 3
> GROUP4 F 4 A 9 11 # Spline for group 4
> GROUP5 G 5 A 10 11 # Spline for group 5
> _____________________________________________________________________
> N.W. Galwey,
> Faculty of Agriculture,
> University of Western Australia,
> Nedlands, WA 6709, Australia.
>
> Tel.: +61 9 380 1959 (direct line)
> +61 9 380 2554 (switchboard)
> Fax: +61 9 380 1108
>
<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>
Arthur Gilmour PhD email: Arthur.Gilmour@agric.nsw.gov.au
Senior Research Scientist (Biometrics) fax: <61> 2 6391 3899
NSW Agriculture <61> 2 6391 3922
Orange Agricultural Institute telephone work: <61> 2 6391 3815
Forest Rd, ORANGE, 2800, AUSTRALIA home: <61> 2 6362 0046
ASREML is currently free by anonymous ftp from pub/aar on ftp.res.bbsrc.ac.uk
or point your web browser at ftp://ftp.res.bbsrc.ac.uk/pub/aar/
To join the asreml discussion list, send the message
subscribe
to asreml-request@chiswick.anprod.CSIRO.au
The address for messages to the list is asreml@chiswick.anprod.CSIRO.au
<> <> <> <> <> <> <>
"Ye are not your own for ye are bought with a price: therefore glorify God
in your body, and in your spirit, which are God's" 1 Corinthians 6:20
<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>