Re: [R] Inconsistent linear model calculations

2008-05-16 Thread e-letter
For the record, v6 represented the averages for columns v3-v5 because I could not quickly find the r instruction to calculate the average. Similiarly v7 represents square power of v2. At the time it was quicker than trawling through five or so user manuals for the formula syntax.

Re: [R] Inconsistent linear model calculations

2008-05-16 Thread Petr PIKAL
Hi Maybe you shall reconsider your decision to avoid manuals especially R-intro.pdf, where you can learn basic operations and data structures. [EMAIL PROTECTED] napsal dne 16.05.2008 09:52:22: For the record, v6 represented the averages for columns v3-v5 because I could not quickly find the

Re: [R] Inconsistent linear model calculations

2008-05-16 Thread Gavin Simpson
On Thu, 2008-05-15 at 14:50 +0100, e-letter wrote: 2008/5/15 Douglas Bates [EMAIL PROTECTED]: Did you happen to notice the part at the bottom of every message about provide commented, minimal, self-contained, reproducible code? Sorry, don't understand what that statement means in my

[R] Inconsistent linear model calculations

2008-05-15 Thread e-letter
Readers, Using version 251 I tried the following command: lm(y~a+b,data=datafile) Resulting in, inter alia: ... coefficients (intercept) a 1.2 3.4 Packages installed: acepack ace() and avas() for selecting regression transformations adlift

Re: [R] Inconsistent linear model calculations

2008-05-15 Thread Douglas Bates
Did you happen to notice the part at the bottom of every message about provide commented, minimal, self-contained, reproducible code? Considering that the result you quote from 251 has 2 coefficients and the result from 171 has 3 coefficients one might contemplate the possibility that you are

Re: [R] Inconsistent linear model calculations

2008-05-15 Thread Michael Kubovy
Hi, To get help with this problem, you will have to create an example that others can duplicate. That is why each message to the list says (at the bottom): PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained,

Re: [R] Inconsistent linear model calculations

2008-05-15 Thread e-letter
2008/5/15 Douglas Bates [EMAIL PROTECTED]: Did you happen to notice the part at the bottom of every message about provide commented, minimal, self-contained, reproducible code? Sorry, don't understand what that statement means in my context. Considering that the result you quote from 251 has

[R] Inconsistent linear model calculations

2008-05-15 Thread e-letter
Professor Kubovy As you instructed, below is the command terminal output: sessionInfo() R version 2.5.1 (2007-06-27) i586-mandriva-linux-gnu locale:

Re: [R] Inconsistent linear model calculations

2008-05-15 Thread Ben Bolker
e-letter inpost at gmail.com writes: 2008/5/15 Douglas Bates bates at stat.wisc.edu: Did you happen to notice the part at the bottom of every message about provide commented, minimal, self-contained, reproducible code? Sorry, don't understand what that statement means in my context.

Re: [R] Inconsistent linear model calculations

2008-05-15 Thread e-letter
On 15/05/2008, Michael Kubovy [EMAIL PROTECTED] wrote: Hi, Your reply doesn't help you with your lm() problem. You need to show the problem with lm() as well. Your original message was unclear on what is wrong. The problem I have with lm() is that when I enter this command, I obtain

Re: [R] Inconsistent linear model calculations

2008-05-15 Thread Kenn Konstabel
On Thu, May 15, 2008 at 3:05 PM, e-letter [EMAIL PROTECTED] wrote: Using version 251 I tried the following command: lm(y~a+b,data=datafile) Resulting in, inter alia: ... coefficients (intercept) a 1.2 3.4 [...] When using version 171 I entered the same command:

Re: [R] Inconsistent linear model calculations

2008-05-15 Thread Rolf Turner
On 16/05/2008, at 7:23 AM, e-letter wrote: Well, I sought a statistician who showed me how to do this task using mathematica. Despite advice otherwise, You are certainly getting bad advice! I am interested to learn how to perform this basic task using r because it's free and and

Re: [R] Inconsistent linear model calculations

2008-05-15 Thread e-letter
Below is direct copy from command terminals of both pcs (mandrake 92 with r 171; mandriva 2008 with r 251, respectively). R : Copyright 2003, The R Development Core Team Version 1.7.1 (2003-06-16) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under

Re: [R] Inconsistent linear model calculations

2008-05-15 Thread Charilaos Skiadas
On May 15, 2008, at 5:37 PM, e-letter wrote: Below is direct copy from command terminals of both pcs (mandrake 92 with r 171; mandriva 2008 with r 251, respectively). R : Copyright 2003, The R Development Core Team Version 1.7.1 (2003-06-16) R is free software and comes with ABSOLUTELY NO

Re: [R] Inconsistent linear model calculations

2008-05-15 Thread Rolf Turner
On 16/05/2008, at 9:44 AM, Charilaos Skiadas wrote: This is exactly why we asked you for a reproducible example and full code. That last entry in the V6 column is 69.65667 in this case, but 66.27667 in the other cases. So you clearly are working with two slightly different dodgy files,

Re: [R] Inconsistent linear model calculations

2008-05-15 Thread e-letter
Thank you to all you eagle eyes; amendment made accordingly and solved. Not sure how the difference occurred... [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] Inconsistent linear model calculations

2008-05-15 Thread Kingsford Jones
Also, it's worth pointing out the reason for the numerical instability of the parameter estimates: the predictors are nearly collinear. (dubious - read.table('clipboard')) V1 V2V3V4V5 V6 V7 1 1 300 39.87 39.85 39.90 39.87333 9 2 2 400 45.16 45.23 45.17 45.18667