请高人指点程序

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我写了这么一段代码：
library( systemfit )
data( rd )
bas.formula <- bas ~(a1*S1+a2*S2)*SCI^a3*(1-exp(-(k1*S1+k2*S2)*(SDI/10000)^k3*AGE))^c
vol.formula <- vol ~ bas*H(d0/(H+d1))
labels <- list( "bas.growth", "vol.growth" )
inst <- ~ SCI + SDI + AGE + H
start.values <- c(a1=20, a2=20, a3=0.1, k1=3, k2=3,
+ k3=4, c=0.1, d0=20, d1=30 )
model <- list( bas.formula, vol.formula )
model.ols <- nlsystemfit( "OLS", model, start.values, data=rd, eqnlabels=labels )
提示错误：
The following object(s) are masked from 'data (position 3)':

    AGE, DBH, H, NHA, S1, S2
The following object(s) are masked from 'rd':

    AGE, DBH, H, NHA, S1, S2
错误于as.matrix(eval(as.formula(eqns[])[[2]])) :
  在为'as.matrix'函数选择方法时评估'x'参数出了错: 错误于eval(expr, envir, enclos) : 找不到对象'bas'

上面的bas.formula和vol.formula为两模型，数据里面只有AGE, DBH, H, NHA, S1, S2，
bas是方程1的因变量作为方程2的自变量，是不是用到非线性度量误差联立方程组，请教您怎么写程序才能正确。希望得到您的指点。

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：请高人指点 高人指点 Library Formula System library 程序

就像范例一样
data( ppine )
names(ppine)
#[1] "elev" "smi"  "dbh"  "tht"  "cr"   "ba"   "dg"   "hg"
hg.formula <- hg ~ exp( h0 + h1*log(tht) + h2*tht^2 + h3*elev + h4*cr)
dg.formula <- dg ~ exp( d0 + d1*log(dbh) + d2*hg + d3*cr + d4*ba )

你的数据里要有AGE,DBH,H,NHA,S1,S2,SCI,SDI,bas,vol
bas.formula <- bas ~(a1*S1+a2*S2)*SCI^a3*(1-exp(-(k1*S1+k2*S2)*(SDI/10000)^k3*AGE))^c
vol.formula <- vol ~ bas*H(d0/(H+d1))

epoh 发表于 2012-3-22 18:23
就像范例一样
data( ppine )
names(ppine)

谢谢您的指点，我还是不能清除，bas 和vol需要AGE,DBH,H,NHA,S1,S2,SCI,SDI求出来，数据中没有呀？是不是数据只有表头有bas,vol，bas和vol行下面都是空格？

xingzhaoh 发表于 2012-3-22 18:59
谢谢您的指点，我还是不能清除，bas 和vol需要AGE,DBH,H,NHA,S1,S2,SCI,SDI求出来，数据中没有呀？是不是 ...

那bas和vol,你就要事先求出,不能空白

xingzhaoh 发表于 2012-3-22 18:59
谢谢您的指点，我还是不能清除，bas 和vol需要AGE,DBH,H,NHA,S1,S2,SCI,SDI求出来，数据中没有呀？是不是 ...

您能说的详细一点吗？初学，还是不怎么懂？谢谢

xingzhaoh 发表于 2012-3-22 19:40
您能说的详细一点吗？初学，还是不怎么懂？谢谢

麻烦你你把数据传上来
并把你求response bas,vol的公式列出来

谢谢

xingzhaoh 发表于 2012-3-22 20:39
公式为：bas ~(a1*S1+a2*S2)*SCI^a3*(1-exp(-(k1*S1+k2*S2)*(SDI/10000)^k3*AGE))^c
vol ~ bas*H(d0/(H+d1 ...

Linear Regression
  Model : Y = beta0+ beta1*A + beta2*B
  Syntax : Y ~ A + B
要求beta0,beta1,beta2,就要给出Y,A,B
仅给A,B如何求出beta0,beta1,beta2

同理你给的是bas的regression formula
bas ~(a1*S1+a2*S2)*SCI^a3*(1-exp(-(k1*S1+k2*S2)*(SDI/10000)^k3*AGE))^c
但没给出bas,如何求a1,a2,a3,k1,k2,k3

" <- "," ~ "两者是不同的

epoh 发表于 2012-3-22 21:09
Linear Regression
  Model : Y = beta0+ beta1*A + beta2*B
  Syntax : Y ~ A + B

对呀，哎糊涂了，谢谢你的指点

xingzhaoh 发表于 2012-3-22 21:16
对呀，哎糊涂了，谢谢你的指点

今天试了一下，又有问题向您询问了老师，程序：
library( systemfit )
rd<-read.delim("甘肃黑龙江分析02.txt")
attach(rd)
BAS.formula <- BAS ~(a1*S1+a2*S2)*SCI^a3*(1-exp(-(k1*S1+k2*S2)*(SDI/10000)^k3*AGE))^c
VOL.formula <- VOL~ BAS *H*(d0/(H+d1))
labels <- list( "bas.growth", "vol.growth" )
inst <- ~ SCI + SDI + AGE +H + BAS+VOL
start.values <- c(a1=10, a2=10, a3=0.05, k1=1, k2=1,k3=0.1, c=0.01,
d0=10, d1=20 )
model <- list( BAS.formula, VOL.formula )
model.ols <- nlsystemfit( "OLS", model, start.values, data=rd, eqnlabels=labels )
到此处出现
错误于qr.default(jacobian) : 外接函数调用时不能有NA/NaN/Inf(arg1)


数据见附件