気候学的な気圧の値と、気圧が測定された観測所の高さのデータがあります。それらに指数モデルを当てはめたいと思います。データは次のようになります。
x
[1] 539 575 1320 438 1840 496 316 552 325 386 1599 1073 556 1029 1661
[16] 2472 1594 1197 910 1035 596 646 420 516 1980 1045 2287 440 419 1611
[31] 577 3580 484 1018 1669 745 1974 366 273 454 203 588 1427 405 1403
[46] 485 490 2106 990 3305 1078 455 300 1638 1708 438 1303 482 775 2502
[61] 457 2690 422 1638 555 426
y
[1] 954.1 951.4 867.2 964.0 813.3 958.8 979.7 950.8 978.4 971.3 835.1 894.1
[13] 952.0 904.4 833.3 751.5 839.0 882.5 912.0 899.4 947.1 942.3 968.5 961.9
[25] 801.3 893.6 769.8 965.6 965.1 836.9 949.2 653.6 959.8 900.2 830.6 928.6
[37] 800.3 971.1 983.5 963.4 992.6 947.5 848.3 969.4 858.2 959.9 959.3 787.2
[49] 900.4 677.6 893.2 962.7 981.5 834.9 827.0 966.0 870.1 961.1 925.2 749.3
[61] 962.8 734.0 968.0 836.3 950.4 966.5
私は最初に、データの対数を取り、それらに適合させようとしlmました:
log.p=log(y)
log.height=log(x)
lmlog=lm(log.p~log.height)
しかし、これはまったく適合しないモデルを提供したため、nls他の投稿 (「開始」など) からさまざまなヒントを得て、関数を使用することにしました。
f <- function(x,a,b) {a* exp(b * x)}
dat <- as.list(x, y)
start <- coef(nls(log(y) ~ log(f(x, a, b)), dat, start = c(a = 1, b = -1)))
nls=nls(y~ f(x,a,b), data=dat, start=start)
残念ながら、「開始」しても以下のようなエラーが出てしまい、本当にどうしたらいいのかわかりません…
Error in numericDeriv(form[[3L]], names(ind), env) :
Missing value or an infinity produced when evaluating the model
In addition: Warning messages:
1: In min(x) : no non-missing arguments to min; returning Inf
2: In max(x) : no non-missing arguments to max; returning -Inf
誰でも助けることができますか？前もって感謝します！！！