Corrected bug in select_weight formula

This commit is contained in:
2025-05-26 11:58:04 +02:00
parent 0a8437fc48
commit cdda1aab85

View File

@ -17,7 +17,12 @@ select_weight <- function(df=NULL,grp=NULL,X=NULL,instr=NULL) {
formu <- paste0(grp,"~") formu <- paste0(grp,"~")
formu2 <- paste(X,sep="+",collapse="+") formu2 <- paste(X,sep="+",collapse="+")
formu3 <- paste(instr,sep="+",collapse="+") formu3 <- paste(instr,sep="+",collapse="+")
formu2 <- paste(formu2,formu3,sep="+") if (!is.null(instr)) {
formu2 <- paste(formu2,formu3,sep="+")
}
if (is.null(X)) {
formu2 <- formu3
}
formu <- paste(formu,formu2) formu <- paste(formu,formu2)
logit_mod <- glm(formula = formu,data = df,family = binomial(link = "probit")) logit_mod <- glm(formula = formu,data = df,family = binomial(link = "probit"))
res <- residuals(logit_mod) res <- residuals(logit_mod)