nls函数做非线性回归

hongyashan123 发表于 2017-9-4 16:39
你把数据的和模型发过来

data <- read.csv("shuju1.csv",header = T)
data <- data[-which(dataB2grWindSpeed=="****"),]&nbsp;##&nbsp;剔除缺失值&nbsp;
Speed&nbsp;<-&nbsp;as.numeric(as.character(data
B2grWindSpeed=="****"),]&nbsp;##&nbsp;剔除缺失值&nbsp;Speed&nbsp;<-&nbsp;as.numeric(as.character(data
B2grWindSpeed)) ##　风速转换为数值型
Process <- as.numeric(as.character(dataB2grGenPowerForProcess))&nbsp;##&nbsp;功率转换为数值型&nbsp;
Direction&nbsp;<-&nbsp;as.numeric(as.character(data
B2grGenPowerForProcess))&nbsp;##&nbsp;功率转换为数值型&nbsp;Direction&nbsp;<-&nbsp;as.numeric(as.character(data
B2grWindDirection)) ## 风向转换为数值型
Speed = Speed[-which(is.na (Process))]
Direction = Direction[-which(is.na (Process))]
Process = Process[-which(is.na (Process))]
p_max = max(p)
start = list(a = 1, b=4)
Mode2 <- nls(p~p_max*(1-exp(-(v/b)^a)), start = start)
M2 <- summary(Mode2)
M2$coefficients ## 显示拟合参数

你的数据中只有两个变量啊

hongyashan123 发表于 2017-9-6 13:59
你的数据中只有两个变量啊

恩，v和P

hongyashan123 发表于 2017-9-6 13:59
你的数据中只有两个变量啊

后面没有用到，我就没有弄另外的一个变量。因为数据太多了

hongyashan123 发表于 2017-9-4 14:27
是R还是matlab中的函数？

R

你的数据读出来有问题啊