[R] Bug or my misunderstanding?

Bert Gunter gunter.berton at gene.com
Mon Jan 27 18:53:08 CET 2014


It is! I apologize for the noise.

Cheers,
Bert

Bert Gunter
Genentech Nonclinical Biostatistics
(650) 467-7374

"Data is not information. Information is not knowledge. And knowledge
is certainly not wisdom."
H. Gilbert Welch




On Mon, Jan 27, 2014 at 9:48 AM, Brian Ripley <ripley at stats.ox.ac.uk> wrote:
> I think it is already a bug report.  But I am away from base and cannot check easily.
>
>> On 27 Jan 2014, at 18:15, Bert Gunter <gunter.berton at gene.com> wrote:
>>
>> Folks:
>>
>> Before I waste someone's time with a stupid bug report, could I get
>> feedback as to whether the following really appears to be a (minor)
>> bug?
>>
>> Summary: get_all_vars does not seem to handle multiple responses correctly.
>>
>> Example:
>>
>>> y <- matrix(runif(12),nc=3)
>>> x <- 1:4
>>> lmfit <-lm(y~x)
>>
>>> model.frame(lmfit) ## OK
>>         y.1        y.2        y.3 x
>> 1 0.02159809 0.15593110 0.59007262 1
>> 2 0.91169201 0.30725236 0.41035328 2
>> 3 0.45079051 0.29174545 0.18771042 3
>> 4 0.07983415 0.37301448 0.70319143 4
>>
>>> get_all_vars(lmfit)## not OK ?
>>           y         x        NA NA
>> 1 0.02159809 0.1559311 0.5900726  1
>> 2 0.91169201 0.3072524 0.4103533  2
>> 3 0.45079051 0.2917455 0.1877104  3
>> 4 0.07983415 0.3730145 0.7031914  4
>>
>> ## model.frame() gives correct response variable names; get_all_vars() does not.
>>
>> R version 3.0.2 (2013-09-25)
>> Platform: i386-w64-mingw32/i386 (32-bit)
>>
>> Many thanks.
>>
>> -- Bert
>>
>> Bert Gunter
>> Genentech Nonclinical Biostatistics
>> (650) 467-7374
>>
>> "Data is not information. Information is not knowledge. And knowledge
>> is certainly not wisdom."
>> H. Gilbert Welch
>>
>> ______________________________________________
>> R-help at r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.




More information about the R-help mailing list