2

nlrq()私は（quantregパッケージから）非線形分位点回帰モデルを当てはめようとしています。しかし、関数が開始値に非常に敏感であることがわかりました。そこでセルフスターターSSexp()を使ってみました。セルフスターター関数が作成したことを読んだことnls()も、それで機能します。しかし、常に次のエラーメッセージが表示されます。

"Error in getInitial(formula, mf) : unused argument (mf)"

私の入力コードは次のとおりです。

fit12 <- nlrq(visit.rate ~ SSexp(het.total, y0, b), 
                               data = flower_mat,tau = 0.75)

getInitialセルフスターター機能を作成するために何らかの形で接続されている漠然とした考えがあります。しかし、率直に言って、どこに問題があるのか​​ わかりません。誰でもこれに関して私を助けることができますか?

追加情報

# R version 3.1.2 -- "Pumpkin Helmet"
# quantreg version: 5.05
# nlrwr version : 1.1-0

# toy dataset and code
x <- c(0.3,0.6,0.9,1,1.5,2, 2.1, 2.5,3, 3,5,10,11,12,14,13,17,21,23,27,30, 50)
y <- c(0,0.1,0.3,0.4,0.6,0.2,0.27,0.2,0.25,0.4,0.15,0.05,0.25,0.2,0.3,0.35,0.1,0.15,0.1,0.14,0,0)
dat <- data.frame(x,y,stringsAsFactors=FALSE)
plot(y~x)
library(quantreg)
library(nlrwr)
fit.1 <- nlrq(y ~ SSexp(x, y0, b), 
              data = dat,tau = 0.75)
# above mentioned error message
with(dat, SSexp(x, 0.2, 4))
# this is working
getInitial((y ~ SSexp(x, y0, b)), dat)
# showing following error:
# Error in getInitial((y ~ SSexp(x, y0, b)), data = dat) : unused argument (dat)
4

2 に答える 2

0

問題は、ロード時にロードされる依存パッケージの 1 つが でnlrwrありdlc、独自のバージョンのgetInitial.

試してみてください stats::getInitial(etc) -あなたのデータとコードで試してみると、「Error in lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) : NA/NaN/Inf in 'y'"

于 2014-12-27T14:57:55.197 に答える