私はRの新しいユーザーであり、私の質問がばかげている場合は、あなたが私に耐えてくれることを願っています。Rの最尤推定量を使用して次のモデルを推定したいと思います。
y= a+b*(lnx-α)
ここで、a、b、およびαは推定されるパラメーターであり、XおよびYは私のデータセットです。Webから取得した次のコードを使用しようとしました。
library(foreign)
maindata <- read.csv("C:/Users/NUNU/Desktop/maindata/output2.csv")
h <- subset(maindata, cropid==10)
library(likelihood)
modelfun <- function (a, b, x) { b *(x-a)}
par <- list(a = 0, b = 0)
var<-list(x = "x")
par_lo <- list(a = 0, b = 0)
par_hi <- list(a = 50, b = 50)
var$y <- "y"
var$mean <- "predicted"
var$sd <- 0.815585
var$log <- TRUE
results <- anneal(model = modelfun, par = par, var = var,
source_data = h, par_lo = par_lo, par_hi = par_hi,
pdf = dnorm, dep_var = "y", max_iter = 20000)
データは異なりますが、つまり、cropidを変更した場合でも、得られる結果は似ています。同様に、生成される予測値は、yではなくxに対するものです。何を逃したのか、何がうまくいかなかったのかわかりません。あなたの助けは大歓迎です。