R error, ## Predicted probability of being assigned to RHC rhc$pRhc <- predict(psModel, type = "response")

[wwh199459]

rhc$pRhc <- predict(psModel, type = "response")
Error in $<-.data.frame(*tmp*, pRhc, value = c(1100109 = 0.82757271771222, :
replacement has 418 rows, data has 448

[ndevln]

Hi,

do you want to tell us that the SMD vignette does not build for you? What do you want to do and when does this problem occur? Do you have problems installing the package?

This minimal example works for me:

rhc <- read.csv("http://biostat.mc.vanderbilt.edu/wiki/pub/Main/DataSets/rhc.csv")

rhc$swang1 <- factor(rhc$swang1, levels = c("No RHC", "RHC"))
## Fit model
psModel <- glm(formula = swang1 ~ age + sex + race + edu + income + ninsclas +
                 cat1 + das2d3pc + dnr1 + ca + surv2md1 + aps1 + scoma1 +
                 wtkilo1 + temp1 + meanbp1 + resp1 + hrt1 + pafi1 +
                 paco21 + ph1 + wblc1 + hema1 + sod1 + pot1 + crea1 +
                 bili1 + alb1 + resp + card + neuro + gastr + renal +
                 meta + hema + seps + trauma + ortho + cardiohx + chfhx +
                 dementhx + psychhx + chrpulhx + renalhx + liverhx + gibledhx +
                 malighx + immunhx + transhx + amihx,
               family  = binomial(link = "logit"),
               data    = rhc)
## Predicted probability of being assigned to RHC
rhc$pRhc <- predict(psModel, type = "response")

[wwh199459]

Dear Dr. Alexander Hello, I am Wang Wenhao, Department of Obstetrics and Gynecology, Second Hospital of Shanxi Medical University. Thank you for your suggestion. Thank you for the R package, let us process data easier. Yours sincerely student Wang Wenhao. Wenhao Wang 2020.09.07 | | 王文豪 | | wwh199459@163.com | 签名由网易邮箱大师定制 On 9/5/2020 23:21，Alexander Bartel<notifications@github.com> wrote： Hi, do you want to tell us that the SMD vignette does not build for you? What do you want to do and when does this problem occur? Do you have problems installing the package? This minimal example works for me: rhc <- read.csv("http://biostat.mc.vanderbilt.edu/wiki/pub/Main/DataSets/rhc.csv") rhc$swang1 <- factor(rhc$swang1, levels = c("No RHC", "RHC")) ## Fit model psModel <- glm(formula = swang1 ~ age + sex + race + edu + income + ninsclas + cat1 + das2d3pc + dnr1 + ca + surv2md1 + aps1 + scoma1 + wtkilo1 + temp1 + meanbp1 + resp1 + hrt1 + pafi1 + paco21 + ph1 + wblc1 + hema1 + sod1 + pot1 + crea1 + bili1 + alb1 + resp + card + neuro + gastr + renal + meta + hema + seps + trauma + ortho + cardiohx + chfhx + dementhx + psychhx + chrpulhx + renalhx + liverhx + gibledhx + malighx + immunhx + transhx + amihx, family = binomial(link = "logit"), data = rhc) ## Predicted probability of being assigned to RHC rhc$pRhc <- predict(psModel, type = "response") — You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub, or unsubscribe.