Мне нужно вычислить коэффициенты в уравнении y = 1000 * ln (Ax + 1)
Я передал данные таблицы с помощью x и y и пытаюсь провести их через R для вычисления коэффициентов.
Вот мой код.
TestModel=read.csv("bdata.csv",header=TRUE)
plot(TestModel$mon,TestModel$count)
TestModel$P = 1000*(log(abs(TestModel$mon+1)))
#TestModel$Q = 1000*(log1p(TestModel$mon))
head(TestModel)
#summary(1000*(log10(1+TestModel$mon)))
output1 = lm(TestModel$count~TestModel$P)
summary(output1)
Кажется, что значение A рассчитывается как 0,47508, но для этого уравнения Excel дает постоянное значение как 0,1432, что является правильным.Как я могу получить это значение в R?что мне не хватает?любая помощь высоко ценится.
Результат:
> head(TestModel)
count month P Q
1 134 1 693.1472 0.6931472
2 252 2 1098.6123 1.0986123
3 357 3 1386.2944 1.3862944
4 453 4 1609.4379 1.6094379
5 540 5 1791.7595 1.7917595
6 620 6 1945.9101 1.9459101
Residuals:
Min 1Q Median 3Q Max
-41.773 -32.478 -5.554 23.566 74.534
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -269.83135 40.65118 -6.638 5.80e-05 ***
TestModel$P 0.47508 0.02073 22.919 5.64e-10 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 40.25 on 10 degrees of freedom
Multiple R-squared: 0.9813, Adjusted R-squared: 0.9795
F-statistic: 525.3 on 1 and 10 DF, p-value: 5.643e-10