Re: [R] Summary information by groups programming assitance

2008-12-22 Thread Ranney, Steven
lect Lake, psd as Category, max(Length), max(vol) from DF group by Lake, psd") There are many other good solutions too using various packages which have already been mentioned on this thread. On Mon, Dec 22, 2008 at 4:51 PM, Ranney, Steven wrote: > All - > > I have data that looks li

[R] Summary information by groups programming assitance

2008-12-22 Thread Ranney, Steven
All - I have data that looks like psd Species Lake Length WeightSt.weightWr Wr.1 vol 432 substock SMB Clear150 41.00 0.01 95.12438 95.10118 0.0105 433 substock SMB Clear152 39.00 0.01 86.72916 86.70692 0.0105 434 substock

[R] Contextstack overlow

2008-12-18 Thread Ranney, Steven
All - I have a number of rows that I am assigning length classes to via l.class<-with(wae, ifelse((Length>=120)&(Length<130),"125", ifelse((Length>=130)&(Length<140),"135", ifelse((Length>=140)&(Length<150),"145", ifelse((Length>=150)&(Length<160),"155", ifelse((Length>=160)&(Length<17

[R] nls.control()

2008-09-03 Thread Ranney, Steven
>> >> and I'm trying to fit a simple Von Bertalanffy growth curve with program: >> >> VonB=nls(TL~Linf*(1-exp(-k*(age-t0))), data=box5.4, >> start=list(Linf=1000, k=0.1, t0=0.1), trace=TRUE) >>... >> Everything works as it should right up until the confint(VonB) statement. >>When I ask >> f

[R] nls.control()

2008-09-02 Thread Ranney, Steven
All - I have data: TL age 388 4 418 4 438 4 428 5 539 10 432 4 444 7 421 4 438 4 419 4 463 6 423 4 ... [truncated] and I'm trying to fit a simple Von Bertalanffy growth curve with program: #Creates a Von Bertalanffy growth model VonB=nls(TL~

[R] lm() and dffits

2008-08-29 Thread Ranney, Steven
All - My question is a bit involved, so bear with me. I have some data that looks like: LakeLL LW 81 2.176091259 1.342422681 81 2.176091259 1.414973348 81 2.176091259 1.447158031 81 2.181843588 1.414973348 81 2.181843588 1.447158031 81 2.

[R] Newbie programming help

2008-08-22 Thread Ranney, Steven
All - Not sure if this is a real programming question, but here goes: I have data that looks like LakeLength Weight 1 158 45 1 179 70 1 200 125 1 202 150 1 206 145 1 209 165 1 210 140 1 215 175 1 216

[R] Confidence intervals with nls()

2008-08-01 Thread Ranney, Steven
I have data that looks like O.lengthO.age 176 1 179 1 182 1 ... 493 5 494 5 514 5 606 5 462 6 491 6 537 6 553 6 432 7 522 7 625 8 661 8 687 10 704 10 615 12 (truncated) with a simple VonB growth model from within nls(): plot

Re: [R] lines() warning message

2008-07-03 Thread Ranney, Steven
7 PM To: Ranney, Steven Cc: r-help@r-project.org Subject: Re: [R] lines() warning message Ranney, Steven wrote: > I have data that looks like > > Year,Recruit,Spawner,Mtempcv > 1958,4532,775,0.24125 > 1959,22996,2310,0.16319 > 1960,628,2990,0.46056 > 1961,879,14

[R] lines() warning message

2008-07-03 Thread Ranney, Steven
I have data that looks like Year,Recruit,Spawner,Mtempcv 1958,4532,775,0.24125 1959,22996,2310,0.16319 1960,628,2990,0.46056 1961,879,1400,0.33028 1962,14747,1130,0.22618 1963,13205,790,0.20596 1964,31793,1195,0.19229 1965,10621,981,0.20363 1966,22271,870,0.3452 1967,8736,1104,0.27511 1968,8761,88

[R] lm() question

2008-07-02 Thread Ranney, Steven
I have data that looks like YC Age Num 82 11 2 83 10 0 84 9 8 85 8 21 86 7 49 87 6 18 88 5 79 89 4 28 90 3 273 91 2 175 with a program mod1=lm(log(Num+1)~YC, data=box44) plot(lo

Re: [R] Quartile regression question

2008-06-13 Thread Ranney, Steven
Thanks for your help. Worked great. SR Steven H. Ranney Graduate Research Assistant (Ph.D) USGS Montana Cooperative Fishery Research Unit Montana State University PO Box 173460 Bozeman, MT 59717-3460 phone: (406) 994-6643 fax: (406) 994-7479 [[alternative HTML version deleted]] _

[R] Quartile regression question

2008-06-13 Thread Ranney, Steven
I have data that looks like lake,loglength,logweight 1,2.369215857,1.929418926 1,2.426511261,2.230448921 1,2.434568904,2.298853076 1,2.437750563,2.298853076 1,2.442479769,2.230448921 1,2.445604203,2.356025857 ... 102,2.722633923,3.310268367 102,2.781755375,3.502153893 102,2.836324116,3.683407299 1

Re: [R] Predicting from an nls model

2008-06-12 Thread Ranney, Steven
phone: (406) 994-6643 fax: (406) 994-7479 -Original Message- From: Prof Brian Ripley [mailto:[EMAIL PROTECTED] Sent: Thu 6/12/2008 7:00 AM To: Ranney, Steven Cc: r-help@r-project.org Subject: Re: [R] Predicting from an nls model To predict from Weight~alpha*Length^beta you need to

[R] Predicting from an nls model

2008-06-12 Thread Ranney, Steven
I keep running up against the same error when I try to plot a line from a nls model. The data is fisheries length/weight data. Code follows: require(graphics) pow = nls(Weight~alpha*Length^beta, data=wae, start=list(alpha=0.001, beta=3.0), trace=TRUE) pre

Re: [R] newbie nls question

2008-06-11 Thread Ranney, Steven
: Michael Manning [mailto:[EMAIL PROTECTED] Sent: Tue 6/10/2008 5:08 PM To: Ranney, Steven; r-help@r-project.org Subject: Re: [R] newbie nls question Hi On 11/06/2008 at 10:56 a.m., "Ranney, Steven" <[EMAIL PROTECTED]> wrote: > I'm tyring to fit a relatively simple nls mode

[R] newbie nls question

2008-06-10 Thread Ranney, Steven
I'm tyring to fit a relatively simple nls model to some data, but keep coming up against the same error (code follows): Oto=nls(Otolith ~ Linf*(1-exp(-k(AGE-to))), data = ages, start = list(Linf=1000, k=0.1, to=0.1), trace = TRUE) The error message I keep getting is "Error in eva