[18979] in s-news-athena
[S] Need help on NLME procedure
daemon@ATHENA.MIT.EDU (Perevozskaya, Inna T.)
Thu Aug 19 15:42:27 1999
Message-Id: <199908191937.OAA13069@wubios.wustl.edu>
Date: Thu, 19 Aug 1999 15:33:08 -0400
From: "Perevozskaya, Inna T." <inna_pereyozskaya@merck.com>
To: "'s-news@wubios.wustl.edu'" <s-news@wubios.wustl.edu>
Mime-Version: 1.0
Content-Type: text/plain
Content-Transfer-Encoding: 7BIT
I am trying to run the following program( Splus 4.5 under Windows):
G1<-function(AGE,b0,b1,A)
{
tmp_A*exp(exp(b0)/b1*(exp(b1*AGE)-1));
G1_tmp;
G1;
}
G1(1,2,3,4);
names(last2);
last2.grouped<-groupedData(HTCM~AGE|ALLOC,data=last2, outer=~GENDER);
last2.fit<-nlme(model=HTCM~G1(AGE,b0,b1,A),
fixed=b0+b1+A~1,data=last2.grouped);
and get an error message "agrument model= not matched"
The data set last2 has 4 variables: ALLOC, GENDER, AGE, HTCM
and looks like(first 20 rows) :
5 1 9.76 129.60
5 1 10.32 132.90
5 1 10.93 134.50
5 1 11.86 138.60
5 1 12.77 142.80
5 1 13.86 147.90
5 1 14.79 156.50
5 1 15.94 165.00
5 1 17.00 167.20
5 1 18.05 168.30
6 1 9.28 143.90
6 1 10.11 148.90
6 1 10.68 151.70
6 1 11.52 156.80
6 1 12.62 164.10
6 1 13.57 172.20
6 1 14.72 179.60
6 1 15.72 182.60
6 1 16.74 185.40
6 1 17.37 183.30
...........
How can I fix it?
Thanks,
Inna Perevozskaya
-----------------------------------------------------------------------
This message was distributed by s-news@wubios.wustl.edu. To unsubscribe
send e-mail to s-news-request@wubios.wustl.edu with the BODY of the
message: unsubscribe s-news
