Question

我需要计算方程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

如何使用R计算log（Ax + 1）中的常数A

0 个答案: