))
text(60,36, paste(Yield at EONR=,yldeonr))
# I'd like to place the superscript on the r
text(60,32, paste(Pseudo r^2=,r_2))
# This does not seem to work ... I can get the subscript
# but now I cannot paste the r_2 object. Any ideas ?
text(60,28, expression(paste(Pseudo r ^2, = r_2)))
--
Jose
)
}
###
###
--
Jose A. Hernandez
Department of Soil, Water, and Climate
University of Minnesota
1991 Upper Buford Circle
St. Paul, MN 55108
Ph. (612) 625-0445, Fax. (612) 625-2208
__
[EMAIL PROTECTED
first. However, R still does not come so
naturally to me when I am coding, I have to look at my old codes,
manuals and google. If I ever end up teaching applied stats I'll teach
using R.
Sincerely,
Jose
... migrating to R ... and one of this days to linux
--
Jose A. Hernandez
Ph.D. Candidate
Please disregard my message I've seen the light now...
curve(predict(out, data.frame(x = x)), add = T)
Thanks,
Jose
--
Jose A. Hernandez
Department of Soil, Water, and Climate
University of Minnesota
1991 Upper Buford Circle
St. Paul, MN 55108
Ph. (612) 625-0445, Fax. (612) 625-2208
R-community,
Assuming a file with 3 columns: site, nrate, yield. In SAS I can easily
run the analysis BY SITE using the following code:
*--*;
proc nlin method=MARQUARDT;
by farm;
parms b0=100 b1=0.33 b2=-0.00111 x0=120;
model yield = (b0 + b1 * nrate
minor0.1
year 2004
month11
day 15
language R
--
Jose A. Hernandez
Ph.D. Candidate
Precision Agriculture Center
Department of Soil, Water, and Climate
University of Minnesota
1991 Upper Buford Circle
St. Paul, MN 55108
Ph. (612) 625-0445, Fax. (612) 625-2208
#
--
Jose A. Hernandez
Ph.D. Candidate
Precision Agriculture Center
Department of Soil, Water, and Climate
University of Minnesota
1991 Upper Buford Circle
St. Paul, MN 55108
Ph. (612) 625-0445, Fax. (612) 625-2208
__
R-help@stat.math.ethz.ch
(131.75,125.92,130.09,42.67,48.26,54.56,85.05,91.17,90.86)
d.1 - data.frame(event,chem)
attach(d.1)
s.err - function(x) sd(x)/(sqrt(length(x)))
print(mean_event - tapply(d.1$chem,factor(data.1$event),mean))
print(s.err_event - tapply(d.1$chem,factor(data.1$event),s.err))
print(total_chem - sum(mean_event))
s.err_total_chem -
--
Jose A. Hernandez
Ph.D
(qp.nls.fit, data.frame(nrate = x)), add = T)
text(100,175, paste(EONR=,eonr2))
text(100,173, paste(Yield at EONR=,yldeonr2))
--
Jose A. Hernandez
Ph.D. Candidate
Precision Agriculture Center
Department of Soil, Water, and Climate
University of Minnesota
1991 Upper Buford Circle
St. Paul, MN 55108
Ph. (612
5
6 5 6
--
Jose A. Hernandez
Ph.D. Candidate
Precision Agriculture Center
Department of Soil, Water, and Climate
University of Minnesota
1991 Upper Buford Circle
St. Paul, MN 55108
Ph. (612) 625-0445, Fax. (612) 625-2208
__
[EMAIL PROTECTED] mailing list
) cumsum(x$y)),
recursive=T))
--
Jose A. Hernandez
Ph.D. Candidate
Precision Agriculture Center
Department of Soil, Water, and Climate
University of Minnesota
1991 Upper Buford Circle
St. Paul, MN 55108
Ph. (612) 625-0445, Fax. (612) 625-2208
[[alternative HTML version deleted
11 matches
Mail list logo