想请问大家如何得到该模型的P值呢？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

输入代码如下：

RData<-read.delim(file="D:/SI data/STG/7 south.txt",header=TRUE)#打开数据文件
> y<-RData$SI
> stg<-RData$STG
> stgi<-log(stg)
> RS<-nls(y~exp(-a*(stgi-b)^2),data=RData,start=list(a=10,b=0.2))
>summary(RS)

得到结果如下：
Formula: y ~ exp(-a * (stgi - b)^2)

Parameters:
  Estimate Std. Error t value Pr(>|t|)   
a  1.15335    0.17014   6.779 3.04e-05 ***
b -0.22398    0.06931  -3.232  0.00799 **
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 0.08747 on 11 degrees of freedom

Number of iterations to convergence: 8
Achieved convergence tolerance: 1.129e-06

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Convergence iterations parameters Parameter Tolerance

1. summary(RS)[["parameters"]][,4]

复制代码

试试看。

把summary(fit)当成一个list,先用name查看有哪些元素名称，然后用美元符号提取，然后再按照数字顺序去找p值

llb_321 发表于 2022-4-19 07:23
试试看。

谢谢老师，非常有用！

jiangbeilu 发表于 2022-4-19 07:33
把summary(fit)当成一个list,先用name查看有哪些元素名称，然后用美元符号提取，然后再按照数字顺序去找p值 ...

谢谢老师的帮助！

用函数str(summary(RS))查看提取p值得结构，summary(RS)$p.value

用嵌套