[问答] 请高人指点程序 [推广有奖]

21楼

xingzhaoh 发表于 2012-3-28 15:52:46

epoh 发表于 2012-3-27 13:41
library( systemfit )
rd

老师，我刚刚用了您给的程序，我怎么还没做出来，还是出现
library( systemfit )
rd<-read.csv("data.csv")
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=2.7, a2=3.8, a3=0.7, k1=1.6, k2=0.0003,k3=0.14, c=0.14,d0=8.7, d1=7.9)
model <- list( BAS.formula, VOL.formula )
model.ols <- nlsystemfit( "OLS", model, start.values, data=rd, eqnlabels=labels )

错误于qr.default(jacobian) : 外接函数调用时不能有NA/NaN/Inf(arg1)

难道是我的版本不对，我的R是2.14.0，还是别的原因，好困惑呀

22楼

xingzhaoh 发表于 2012-3-28 15:54:32

xingzhaoh 发表于 2012-3-28 15:52
老师，我刚刚用了您给的程序，我怎么还没做出来，还是出现
library( systemfit )
rd

您都做出来了，我照抄还是没出做出来，好受打击
刚看到您回复很高兴，现在自己照抄都不能做出来，到底是什么原因呢

23楼

epoh 发表于 2012-3-28 16:10:40

xingzhaoh 发表于 2012-3-28 15:54
您都做出来了，我照抄还是没出做出来，好受打击
刚看到您回复很高兴，现在自己照抄都不能做出来，到底是 ...

你的数据文件,有没混淆?
rd<-read.csv("sysfit_1346.csv")
sysfit_1346.csv

sysfit_1346.rar (19.9 KB) 本附件包括：
sysfit_1346.csv

24楼

xingzhaoh 发表于 2012-3-28 16:11:50

xingzhaoh 发表于 2012-3-28 15:52
老师，我刚刚用了您给的程序，我怎么还没做出来，还是出现
library( systemfit )
rd

如果您方便的话，可以告诉一下QQ号，想让您详细帮我指导一下到底什么原因
谢谢

25楼

xingzhaoh 发表于 2012-3-28 16:48:11

epoh 发表于 2012-3-28 16:10
你的数据文件,有没混淆?
rd

老师这个参数初始值的设置有什么好的方法吗？您能教教我吗？
非常希望能够得到您的指点

26楼

epoh 发表于 2012-3-28 18:55:37

xingzhaoh 发表于 2012-3-28 16:48
老师这个参数初始值的设置有什么好的方法吗？您能教教我吗？
非常希望能够得到您的指点

呵呵
initial values应该是研究者最清楚的,
是随着模型及数据而变,
凭着经验及对模型的了解,try and error而设定的

27楼

xingzhaoh 发表于 2012-3-28 19:04:11

epoh 发表于 2012-3-28 18:55
呵呵
initial values应该是研究者最清楚的,
是随着模型及数据而变,

再向您请教一下数据文本大概有多大，才能在Memory Limits in R 限制之内，不会产生 Error: cannot allocate vector of size 528.4 Mb

另外有没有方法可以对数据大小不受限制？

刚刚学会了，却因为数据量太大受限制，有点不甘心，想请教老师看看有没有解决的方法