help biglm.big.matrix; problem with weights

classic Classic list List threaded Threaded
2 messages Options
Reply | Threaded
Open this post in threaded view
|

help biglm.big.matrix; problem with weights

Rahim Hajji
Hello colleagues,

I  have tried to use the package biglm. I want to specify a
multivariate regression with a weight.

I have imported a large dataset with the library(bigmemory). I load
the library (biglm) and specified a regression with a weight. But I
get everytime a error message like “object not found” or “`weights'
must be a formula” or "error in eval(expr, envir, enclos)". I don’t
understand this messages. Can sombody probably explain me what this
means.

I have read the documentation. I don’t understand the documentation
for biglm.big.matrix. I found a sentence to weights that is as follow:
A one-sided, “weight  = single term formula specifying weights”. Who
knows what that means?

My Coding is:


library(bigmemory)
library(biganalytics)
library(biglm)

help(read.big.matrix)
NIKA <- read.big.matrix("G:\\VAR.dat",  header=TRUE, type="double",
sep="\t")

Reg <- biglm.big.matrix(formula = LAannualisiert ~ 0  + VAR01_2008 +
VAR02_2008 + VAR03_2008 + VAR04_2008 + VAR05_2008 + VAR06_2008 +
VAR07_2008 + VAR08_2008 + VAR09_2008 + VAR10_2008 + VAR11_2008 +
VAR12_2008 + VAR13_2008 +
VAR14_2008 + VAR15_2008 + VAR16_2008 + VAR17_2008 + VAR18_2008 +
VAR19_2008 +
VAR21_2008 + VAR22_2008 + VAR23_2008 + VAR24_2008 + VAR25_2008 +
VAR26_2008 + VAR27_2008 +
VAR28_2008 + VAR29_2008 + VAR30_2008 + VAR31_2008 + VAR32_2008 +
VAR33_2008 + VAR34_2008 +
VAR35_2008 + VAR36_2008 + VAR37_2008 + VAR38_2008 + VAR39_2008,
weights = ~Gewicht , data = NIKA)

summary(Reg)

Has sombody an idea, what my error is?

Thanks in advance for your help!


Best regards form germany, berlin

______________________________________________
[hidden email] 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.
Reply | Threaded
Open this post in threaded view
|

Re: help biglm.big.matrix; problem with weights

Thomas Lumley
On Tue, 15 Jun 2010, Rahim Hajji wrote:

> Hello colleagues,
>
> I  have tried to use the package biglm. I want to specify a
> multivariate regression with a weight.
>
> I have imported a large dataset with the library(bigmemory). I load
> the library (biglm) and specified a regression with a weight. But I
> get everytime a error message like “object not found” or “`weights'
> must be a formula” or "error in eval(expr, envir, enclos)". I don’t
> understand this messages. Can sombody probably explain me what this
> means.
>
> I have read the documentation. I don’t understand the documentation
> for biglm.big.matrix. I found a sentence to weights that is as follow:
> A one-sided, “weight  = single term formula specifying weights”. Who
> knows what that means?
>
> My Coding is:
>
>
> library(bigmemory)
> library(biganalytics)
> library(biglm)
>
> help(read.big.matrix)
> NIKA <- read.big.matrix("G:\\VAR.dat",  header=TRUE, type="double",
> sep="\t")
>
> Reg <- biglm.big.matrix(formula = LAannualisiert ~ 0  + VAR01_2008 +
> VAR02_2008 + VAR03_2008 + VAR04_2008 + VAR05_2008 + VAR06_2008 +
> VAR07_2008 + VAR08_2008 + VAR09_2008 + VAR10_2008 + VAR11_2008 +
> VAR12_2008 + VAR13_2008 +
> VAR14_2008 + VAR15_2008 + VAR16_2008 + VAR17_2008 + VAR18_2008 +
> VAR19_2008 +
> VAR21_2008 + VAR22_2008 + VAR23_2008 + VAR24_2008 + VAR25_2008 +
> VAR26_2008 + VAR27_2008 +
> VAR28_2008 + VAR29_2008 + VAR30_2008 + VAR31_2008 + VAR32_2008 +
> VAR33_2008 + VAR34_2008 +
> VAR35_2008 + VAR36_2008 + VAR37_2008 + VAR38_2008 + VAR39_2008,
> weights = ~Gewicht , data = NIKA)
>
> summary(Reg)
>
> Has sombody an idea, what my error is?

No, but it would help if you said exactly what the error message was for this code.

         -thomas

Thomas Lumley Assoc. Professor, Biostatistics
[hidden email] University of Washington, Seattle

______________________________________________
[hidden email] 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.