Skip to content

error while using olsrr function:  #206

Description

@kaygmz

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

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Fields

    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions