Ошибка в пакете Predict () с использованием macOS - PullRequest
0 голосов
/ 29 апреля 2019

У меня проблема с запуском функции Predict() под MacOS 10.14.4 и R 3.5.3. Я использую тот же синтаксис и DB работает нормально в Windows10:

Predict(lvef, LVEF, ModSev_MRAS_overall, fun=exp)

Я всегда получаю сообщение об ошибке:

Error in match.arg(type) : 'arg' must be NULL or a character vector

Я также попросил других коллег, использующих Mac, и у них возникла та же проблема.

Затем я попытался прочитать справку Predict {rms} и не смог запустить код примеров, я всегда получаю сообщение об ошибке, и он не запускается.

n <- 1000    # define sample size
set.seed(17) # so can reproduce the results
age            <- rnorm(n, 50, 10)
blood.pressure <- rnorm(n, 120, 15)
cholesterol    <- rnorm(n, 200, 25)
sex            <- factor(sample(c('female','male'), n,TRUE))
label(age)            <- 'Age'      # label is in Hmisc
label(cholesterol)    <- 'Total Cholesterol'
label(blood.pressure) <- 'Systolic Blood Pressure'
label(sex)            <- 'Sex'
units(cholesterol)    <- 'mg/dl'   # uses units.default in Hmisc
units(blood.pressure) <- 'mmHg'

# Specify population model for log odds that Y=1
L <- .4*(sex=='male') + .045*(age-50) +
  (log(cholesterol - 10)-5.2)*(-2*(sex=='female') + 2*(sex=='male'))
# Simulate binary y to have Prob(y=1) = 1/[1+exp(-L)]
y <- ifelse(runif(n) < plogis(L), 1, 0)

ddist <- datadist(age, blood.pressure, cholesterol, sex)
options(datadist='ddist')

fit <- lrm(y ~ blood.pressure + sex * (age + rcs(cholesterol,4)))
Predict(fit, age, cholesterol, np=4)

Я получаю сообщение об ошибке:

Error in model.frame.default(Terms, newdata, na.action = na.action, ...) : 
  variable lengths differ (found for '(np)')
Inoltre: Warning messages:
1: In if (conf.int) { :
  la condizione la lunghezza > 1 e solo il promo elemento verrà utilizzato
2: In qnorm((1 + conf.int)/2) : Si è prodotto un NaN

Пожалуйста, помогите! G

...