Re: [R] Maximum likelihood estimation in R

2004-02-16 Thread DivineSAAM
Hello, Excellent, also the book: Pawitan, Yudi (2001). In all Likelihood: Statistical Modelling and Inference using Likelihood, Clarendon Press, Oxford. Is very good and the associated Web Site is full of MLE using R. Hope this also helps. /oal __ [

Re: [R] Maximum likelihood estimation in R

2004-02-15 Thread DivineSAAM
Hello, Use > x=rnorm(100, mean=3, sd=1) > library(MASS) >fitdistr(x, "normal") mean sd 2.9331 0.99673982 (0.09967398) (0.07048015) Hope this helps, Shrieb __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/ma

Re: [R] problems when compiling C code

2004-02-02 Thread DivineSAAM
Hello, Do you have S Programming by Venables and Ripley, Springer, 2000? There is an excellent discussion and examples there on compiling multifile codes. I think your problem is in the order of the compilation of the multiple files. Best, /oal __ [E

[R] MATLAB to R

2004-01-30 Thread DivineSAAM
Ladies and Gentlemen, In MATLAB, I can write: for J=1:M Y(J+1)=Y(J)+ h * feval(f,T(J),Y(J)); ... In R, I can write above as: for (J in 2:M) { y = y + h * f(t,y) ... } In MATLAB, I can write: for J=1:M k1 = feval(f,T(J),Y(J)); k2 = feval(f,T(J+1),Y(J)+ h * k1 How do I write k2 in R? k1 =

Re: [R] Fortran source code

2004-01-27 Thread DivineSAAM
Dear All, After spending 3 long days attempting to interface Fortran with R--having spent 1 week sifting through R-help and the horrific official documentation--I cannot emphasize in words the importance of consulting 1 and-only 1 reference: Venables, W.N., B.D. Ripley, S Programming. Springer,

RE: [R] Fitting compartmental model with nls and lsoda?

2004-01-23 Thread DivineSAAM
Dear Jesus, Yes there is a way and it is via Christoffer Torn\{o}e's package nlmeODE. I checked Chapter 4 of Bates and Watts. As usual, a little work involved, but doable and quite powerful. Thanks, olinares __ [EMAIL PROTECTED] mailing list https://

Re: [R] Fitting compartmental model with nls and lsoda?

2004-01-22 Thread DivineSAAM
Dear Colleagues, Our group is also working on implementing the use of R for pharmacokinetic compartmental analysis. Perhaps I have missed something, but > fit <- nls(noisy ~ lsoda(xstart, time, one.compartment.model, c(K1=0.5, k2=0.5)), +data=C1.lsoda, +start=list(K1=0.3,

Re: [R] ODE solvers in R (was:The Wrong Choice: Locked in by licenserestrict...

2003-05-30 Thread DivineSAAM
Dear Colleague: I am glad to hear from you. I was going to contact you and ask if you have notes on the steps involved on how you implemented LSODA. I will gladly take the project on. In addition to dassl, I want to implement ODESSA also. In other words, my project is a "suite" of ode solvers f

Re: [R] The Wrong Choice: Locked in by license restrictions

2003-05-28 Thread DivineSAAM
In a message dated 5/27/2003 11:18:33 PM Eastern Standard Time, [EMAIL PROTECTED] writes: > Can you comment on the benefits of odepack versus lsoda? The benefit of ODEPACK vs. LSODA is mainly that ODEPACK is a collection of solvers (A. C. Hindmarsh (1983) "ODEPACK: a systematized collection of

Re: [R] The Wrong Choice: Locked in by license restrictions

2003-05-27 Thread DivineSAAM
In a message dated 5/27/2003 7:11:00 PM Eastern Standard Time, [EMAIL PROTECTED] writes: > originalContent/0,289142,sid39_gci902076,00.html I run MATLAB v6.5 Release 13. In my view, the benefit of Matlab over R depends on your objectives. I am now using R exclusively, except for solving differe