-
Notifications
You must be signed in to change notification settings - Fork 23
Open
Description
Hello,
I am using ols_step_forward_aic function from olsrr package to do stepwise AIC forward regression but I am getting an error:
I have attached my datasets for reference:
dput(combine_df[1:5,1:52])
structure(list(ENSG00000130943.6 = c(-0.0814098051060948, 0.206379671030366,
0.321031350919992, 0.513621983895094, 0.659905091586951), rs1 = c(2L,
0L, 1L, 0L, 0L), rs2 = c(2L, 0L, 1L, 0L, 0L), rs3 = c(2L, 0L,
1L, 1L, 1L), rs4 = c(0L, 0L, 0L, 0L, 0L), rs5 = c(0L, 0L, 0L,
1L, 0L), rs6 = c(1L, 0L, 1L, 0L, 0L), rs7 = c(0L, 0L, 1L, 1L,
0L), rs8 = c(0L, 0L, 0L, 0L, 0L), rs9 = c(0L, 0L, 0L, 0L, 0L),
rs10 = c(0L, 2L, 2L, 2L, 1L), rs11 = c(2L, 1L, 1L, 1L, 2L
), rs12 = c(2L, 1L, 0L, 1L, 0L), rs13 = c(0L, 1L, 0L, 0L,
0L), rs14 = c(0L, 1L, 0L, 0L, 0L), rs15 = c(0L, 0L, 0L, 1L,
0L), rs16 = c(0L, 0L, 0L, 0L, 0L), rs17 = c(0L, 1L, 0L, 0L,
0L), rs18 = c(0L, 1L, 0L, 0L, 0L), rs19 = c(0L, 2L, 1L, 2L,
2L), rs20 = c(0L, 2L, 1L, 2L, 2L), rs21 = c(0L, 0L, 1L, 0L,
0L), rs22 = c(0L, 0L, 0L, 0L, 0L), rs23 = c(0L, 0L, 1L, 0L,
0L), rs24 = c(0L, 0L, 1L, 0L, 0L), rs25 = c(0L, 0L, 0L, 0L,
0L), rs26 = c(0L, 0L, 1L, 0L, 0L), rs27 = c(2L, 1L, 1L, 2L,
2L), rs28 = c(0L, 0L, 1L, 0L, 0L), rs29 = c(0L, 0L, 1L, 0L,
0L), rs30 = c(0L, 0L, 1L, 0L, 0L), rs31 = c(0L, 0L, 1L, 0L,
0L), rs32 = c(0L, 0L, 1L, 0L, 0L), rs33 = c(0L, 0L, 1L, 0L,
0L), rs34 = c(0L, 0L, 1L, 0L, 0L), rs35 = c(0L, 0L, 1L, 0L,
0L), rs36 = c(0L, 0L, 1L, 0L, 0L), rs37 = c(0L, 0L, 1L, 0L,
0L), rs38 = c(0L, 0L, 1L, 0L, 0L), rs39 = c(0L, 0L, 0L, 1L,
0L), rs40 = c(0L, 2L, 1L, 1L, 1L), rs41 = c(2L, 2L, 2L, 1L,
2L), rs42 = c(0L, 2L, 2L, 1L, 2L), rs43 = c(0L, 2L, 2L, 1L,
2L), rs44 = c(2L, 1L, 1L, 0L, 0L), rs45 = c(0L, 0L, 0L, 0L,
0L), rs46 = c(0L, 0L, 1L, 0L, 0L), rs47 = c(1L, 0L, 0L, 1L,
1L), rs48 = c(2L, 0L, 0L, 1L, 1L), rs49 = c(0L, 0L, 1L, 0L,
0L), rs50 = c(0L, 0L, 0L, 0L, 1L), rs51 = c(0L, 0L, 0L, 0L,
0L)), row.names = c("GTEX-1117F", "GTEX-111FC", "GTEX-1128S",
"GTEX-117XS", "GTEX-1192X"), class = "data.frame")
model <- lm(ENSG00000130943.6 ~ ., data = combine_df)
> ols_step_forward_aic(model)
Error in if (aics[minaic] < laic[lpreds]) { :
the condition has length > 1
Reactions are currently unavailable
Metadata
Metadata
Assignees
Labels
No labels