I think there is a typo in your reprex l(x^2) ?? mydt[1,2] contains a list. Which when unlisted contains a load of data.
I'm not sure what you are asking for? Are you trying to unlist that and have it as a row? Sort of pivot.wider if you like or unnest in tidyverse concepts? I think the data.table verbs are rbindlist But I haven't really understood what your end "product" looks like. I think it's a two row table, with 20 intercepts and 20 gradients, repeated a second time? I'd have made that a 20 row table with x, y and then 4 columns for the two versions? On Sat, 21 Sep 2024, 22:25 Naresh Gurbuxani, <naresh_gurbux...@hotmail.com> wrote: > I am trying to store regression objects in a data.table > > df <- data.frame(x = rnorm(20)) > df[, "y"] <- with(df, x + 0.1 * x^2 + 0.2 * rnorm(20)) > > mydt <- data.table(mypower = c(1, 2), myreg = list(lm(y ~ x, data = df), > lm(y ~ x + I(x^2), data = df))) > > mydt > # mypower myreg > # <num> <list> > #1: 1 <lm[12]> > #2: 2 <lm[12]> > > But mydt[1, 2] has only the coeffients of the first regression. mydt[2, > 2] has residuals of the first regression. These are the first two > components of "lm" object. > > mydt[1, myreg[[1]]] > #(Intercept) x > # 0.107245 1.034110 > > Is there a way to put full "lm" object in each row? > > Thanks, > Naresh > > ______________________________________________ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > https://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide https://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
