建立gauss library的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我想建立一个内建的程序叫
myprobit(y,x)裡面包含兩個proc
一個是probit的likelihood
一個是qnewton的程序

proc myprobit(y,x);
local b0,coef,g,g,rec;
b0=olsqr(y,x);
{coef,f,g,rec}=qnewton(&probll,b0);
rept(coef);
endp;

proc probll(b);
local prob;
prob= cdfn(x*b);
ll=-sum(y.*log(prob)+(1-y).*log(1-prob));
rept(ll);
endp;


1
把這個程序存成myprobit.src放在gauss\src裡頭

2
在lib\user.lcg加入
myprobit.src
    myprobit                          : proc
    prob_ll                          : proc

3
開始引用
library user;
e=rndn(100,1);
x1=rndu(100,1)
y1=1+2*x+e;
y1=y1 .>0;
coef=myprobit(y1,x1)

4
error: y is undefined symbol
想請教兩個問題:

a: 請問需要加入指令什麼才能讓myprobit的程序
自動認為y=y1,x=x1

b:如果引用qnewton，是不是要在myprobit.src加入
#include qnewton.ext

謝謝 ^^

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：Library BRARY GAUSS Aus bra GAUSS Library

a. gauss不用指令，错误说明定义有误
b. 不用