[回归分析求助] nlsur命令使用 [推广有奖]

11楼

735074371 发表于 2018-12-6 15:57:44

林学习- 发表于 2018-7-13 16:37
请问，楼主现在会了吗？我也在看这个方法但是很少有人用这个方法，请问您可否将nlsur的代码和对应的方程上传 ...

哥们，求分享。

12楼

王zhen 发表于 2019-7-2 11:28:46

你好可以私聊一下吗？扣：1023560168

13楼

jinggege0217 发表于 2019-8-18 11:54:37

王zhen 发表于 2019-7-2 11:28
你好可以私聊一下吗？扣：1023560168

你好，请问nlsur可以请教您么？

14楼

dxq2 发表于 2021-9-30 10:50:11

仪大姗发表于 2018-8-2 15:29
你找到相关资料了吗朋友！能不能分享一下！

请问您找到相关资料了吗？求指点，我的邮箱18706638587@163.com

15楼

dxq2 发表于 2021-9-30 10:51:03

林学习- 发表于 2018-7-13 16:37
请问，楼主现在会了吗？我也在看这个方法但是很少有人用这个方法，请问您可否将nlsur的代码和对应的方程上传 ...

请问您找到相关资料了吗？求指点。我的邮箱18706638587@163.com

16楼

momostudent 发表于 2021-10-29 17:01:08

林学习- 发表于 2018-7-13 16:37
请问，楼主现在会了吗？我也在看这个方法但是很少有人用这个方法，请问您可否将nlsur的代码和对应的方程上传 ...

你的问题处理完了吗，可不可以分享一下呀

17楼

tongjixue2005 发表于 2023-3-6 10:42:25

R 软件计算也很方便，例子如下，供参考
http://127.0.0.1:20985/library/nlsur/html/nlsur.html，这是帮助文件

## https://janmarvin.r-universe.dev/nlsur
## install.packages('nlsur', repos = c('https://janmarvin.r-universe.dev', 'https://cloud.r-project.org'))
library(nlsur)

data(costs)

dat <- costs
# apply a patch to create Greene Ed. 7 Data
dat$Sm[dat$Year == 1958] <- 0.61886
dat$Pe[dat$Year == 1950] <- 1.12442
dat$Pm[dat$Year == 1949] <- 1.06625

# model
model <- list(
  Sk ~ bk + dkk * log(Pk/Pm) + dkl * log(Pl/Pm) + dke * log(Pe/Pm),
  Sl ~ bl + dkl * log(Pk/Pm) + dll * log(Pl/Pm) + dle * log(Pe/Pm),
  Se ~ be + dke * log(Pk/Pm) + dle * log(Pl/Pm) + dee * log(Pe/Pm)
)

erg <- nlsur(eqns = model, data = dat, type = "FGNLS")
erg

## Not run:
# Greene Example 10.3
library(nlsur)

url <- "http://www.stern.nyu.edu/~wgreene/Text/Edition7/TableF10-2.txt"
dd <- read.table(url, header = T)

names(dd) <-
c("Year", "Cost", "Sk", "Sl", "Se", "Sm", "Pk", "Pl", "Pe", "Pm")

eqns <-
list( Sk ~ bk + dkk * log(Pk/Pm) + dkl * log(Pl/Pm) + dke * log(Pe/Pm),
   Sl ~ bl + dkl * log(Pk/Pm) + dll * log(Pl/Pm) + dle * log(Pe/Pm),
   Se ~ be + dke * log(Pk/Pm) + dle * log(Pl/Pm) + dee * log(Pe/Pm))

strtvls <- c(be = 0, bk = 0, bl = 0,
         dkk = 0, dkl = 0, dke = 0,
         dll = 0, dle = 0, dee = 0)

erg <- ?nlsur(eqns = eqns, data = dd, startvalues = strtvls, type ="FGNLS",
         trace = TRUE, eps = 1e-10)

erg