test_reg.R

In stp4/stp25APA2: Resultate Formatiern

require(stpvers)
#library(lmerTest )
#library(lmerTest)
 fit1 <- lm(chol0 ~  ak + rrs0 + med + g, hyper)


Ordnen(fit1)
Ordnen(fit1,
include.b = TRUE,
include.se = TRUE,
include.beta = TRUE,
# include.eta = TRUE,
include.ci = TRUE,
# include.variance = TRUE,
include.r = TRUE,
include.ftest = TRUE)


fit2 <-  glm(chol0 ~ med +   ak +   g + rrs0 , hyper, family = poisson())
fit3 <- lmerTest::lmer(chol0 ~ rrs0 + med +  ak  +  (1 | g) , hyper)
fits <- list(fit1, fit2, fit3)

x <- APA_Table(fit1, fit2, fit3,  
               test.my.fun = TRUE, output = FALSE)

str(x)
# x <-  APA_Table(fit1, fit2, fit3,
#                 type = "long",
#                 test.my.fun = TRUE, output = FALSE)

summary(fit3)
extract_param(fit3 ,  effects = "fixed" )
# effects = c("ran_pars","fixed"),
# scales = NULL, ## c("sdcor",NA),
# ran_prefix=NULL,
# conf.int = FALSE,
# conf.level = 0.95,
# conf.method = "Wald",

stop()
x1<-APA_Table(fit1,
              type="long",
              include.b = TRUE,
              include.se = TRUE,
              include.beta = TRUE,
              include.t =TRUE,
              #  
              # # include.odds = TRUE,
              include.ci = TRUE,
              #  
              #  
              #  #Fehler abfangeb
              include.p = TRUE,
              include.stars =  TRUE,
              #  
              #  include.variance = FALSE,
              include.r = FALSE,
              #  include.pseudo = TRUE,
              #  # noch nicht fertig
              # # include.ftest = TRUE,
              # # include.loglik = TRUE,
              #  # noch nicht fertig
              #  
              #  include.custom = TRUE,
              #  
              include.aic = FALSE,
              include.bic = FALSE,

              include.gof = FALSE,
              
              output = FALSE
)
 


# 
#  #  x2
#  
#   #RMSE( fit3)
#    
#   
#    
#  #  Extract2(fit1) 
#  
#    stp25APA2:::Ordnen.lm(fit1)
   broom::tidy(fit1)
    broom::tidy(fit2)
   broom::tidy(fit3)
#    
#  
#   # stp25APA2:::Extract2.default(fit1)
#    
#    
#    stp25APA2:::test_regression(fit1)   
   x2<-APA_Table(fit1,
                 type="long",
                 include.b = TRUE,
                  include.se = TRUE,
                  include.beta = TRUE,
                  include.t =TRUE,

                 include.ci = TRUE,

                #  #Fehler abfangeb
                 include.p = TRUE,
                 include.stars =  TRUE,
                #  
                 #  include.variance = FALSE,
                  include.r = FALSE,
               
                #  
                   include.aic = FALSE,
                   include.bic = FALSE,
                 
             
                include.param = FALSE,
                
                 output = FALSE,
                test.my.fun=TRUE
   )