我已經用我的資料擬合了一個 nlsList 模型，現在我需要使用我的模型估計對應于給定“y”的“x”值。
我曾嘗試使用投資包??中的“投資”，但它僅適用于來自類 lm、glm、nls 或 lme 的物件。
uj5u.com熱心網友回復：
發布到 SO 時，問題需要提供一個可重現的測驗示例，包括輸入、代碼和預期輸出。請參閱頂部的說明r標簽頁。這次我已經在最后的注釋中為您完成了。
創建一個pred預測Sepal.Width給定Sepal.Length和的函式Species。然后我們可以使用uniroot如圖：
library(nlme)
# run code in Note
pred <- function(Sepal.Length, Species) {
predict(fm, list(Sepal.Length = Sepal.Length, Species = Species))
}
uniroot(function(x) pred(x, "virginica") - 3, c(-10, 10))
這給出了以下內容，即Sepal.Length對應于Sepal.Width3 和弗吉尼亞州的 aSpecies是 6.7168。
$root
[1] 6.716803
$f.root
virginica
2.524802e-06
attr(,"label")
[1] "Predicted values"
$iter
[1] 5
$init.it
[1] NA
$estim.prec
[1] 6.103516e-05
筆記
使用內置=in iris 資料集的可重現示例如下。我們假設給定我們想要找到的 Sepal.Width = 3 和 Species = "virginica" Sepal.Length。
library(nlme)
fm <- nlsList(Sepal.Width ~ a * exp(Sepal.Length*b) | Species, data = iris,
start = list(a = 0.5, b = 0.5))
轉載請註明出處，本文鏈接：https://www.uj5u.com/shujuku/532086.html
標籤：r目的模型预言