AW: [R] Improving R Editing: New JEdit Edit Mode for R - R-help

Thu, Nov 21, 2002 3:16 AM #

Hi Daniel,

I have never tried out WinEdit but after reading your message I took a quick glance at their homepage. It would not be useful for JEdit to 
implement the provided interaction with R-GUI since JEdit has to work on several platforms. But creating a button or shortcut for 
sending a complete file  / line / selection of code to R (the command line version) via pipe and to get back the result into JEdit shouldn't 
really be a problem. JEdit has a powerful macro language providing support for such things.

Best regards, Tobias.


-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._

Andrew Beckerman

Fri, Nov 22, 2002 9:29 AM #

Dear Listers -

R1.61 on NT

We are trying to fit a polynomial usning nlme so that in ax^2+bx+c, c is 
fixed and a and b are random.

The data are
Frames = speed of flight (20-50)
Manlength = tail length (-20,0)
ID = individual (14 of them)

each individual has experienced each of 6 tail lengths in the range above 
and we are interested in the shape of the relationship between speed and 
manipulation.

using the following approach based on the Venables and Ripley (3rd ed) 
examples with the Sitka data set

nlme(Frames~a*Manlength^2+b*Manlength+c,
+ fixed=list(c~1),random=a+b~1|ID,
+ start=list(fixed=c(0.3,6,50)))

Error in nlme.formula(Frames ~ a * Manlength^2 + b * Manlength + c, fixed = 
list(c ~  :
         starting values for the fixed component are not the correct length

-realizing that something is wrong with how we have framed the start 
values... can anyone advise?

cheers
andrew

----------------------------------------------------------
Dr. Andrew Beckerman
Institute of Biological Science
University of Stirling
Stirling FK9 4LA
+44 (0)1786 then wk-467808 fx-464994
(APB is not responsible for anything below this)

The University of Stirling is a university established in Scotland by
charter at Stirling, FK9 4LA.  Privileged/Confidential Information may
be contained in this message.  If you are not the addressee indicated
in this message (or responsible for delivery of the message to such
person), you may not disclose, copy or deliver this message to anyone
and any action taken or omitted to be taken in reliance on it, is
prohibited and may be unlawful.  In such case, you should destroy this
message and kindly notify the sender by reply email.  Please advise
immediately if you or your employer do not consent to Internet email
for messages of this kind.  Opinions, conclusions and other
information in this message that do not relate to the official
business of the University of Stirling shall be understood as neither
given nor endorsed by it.

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._

Douglas Bates

Fri, Nov 22, 2002 11:17 AM #

Andrew Beckerman <a.p.beckerman at stir.ac.uk> writes:

First the good news - you don't need to use nlme to fit a polynomial
model.  A polynomial is linear in the coefficients so you can use lme,
which is much easier than nlme.

lme(Frames ~ Manlength + I(Manlength^2), random = Manlength)

The problem with your usage of nlme is that you gave 3 values for the
starting estimates of the fixed-effects parameters and there is only
one parameter for fixed-effects the way you wrote it.  You should have
used "fixed=list(a+b+c ~ 1)" (or, equivalently, "fixed=a+b+c~1") but
you really don't want to use nlme for a linear model.
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._

Douglas Bates

Fri, Nov 22, 2002 1:42 PM #

Douglas Bates <bates at stat.wisc.edu> writes:

<snip>

<snip>

Oops - I got that wrong.  The random argument should be
  random = ~ Manlength | ID
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._