SAS - Proc NLIN

   
Vista:

Proc NLIN

Publicado por Rodrigo Abad Guamán rodmabthebest@yahoo.com (1 intervención) el 22/05/2014 10:12:20
Hola Necesito correr el siguiente proc pero me sale error

Me podrías ayudar a ver que esta mal

1
2
3
4
5
6
7
8
9
10
11
PROC NLIN DATA=GASPRO BEST=1 ITER=100 CONVERGENCE=.000001
METHOD=	MARQUARDT;
parms b0=40 b1=.03 b2=1 b3=200 b4=.05;
model GAS=(b0/(1+exp(2-4*b1*(t-b2))))+(b3/(1+exp(2-4*b4*(t-b2))));
der.b0=1/(1+exp(2-4*b1*(t-b2)));
der.b1=(4*b0*(t-b2)*(exp(2-4*b1*(t-b2)))/((1+exp(2-4*b1*(t-b2)))*(1+exp(2-4*b1*(t-b2))));
der.b2=((-4*b0*b1*(exp(2-4*b1*(t-b2))))/((1+exp(2-4*b1*(t-b2)))*(1+exp(2-4*b1*(t-b2)))))-((4*b3*b4*(exp(2-4*b4*(t-b2))))/((1+exp(2-4*b4*(t-b2)))*(1+exp(2-4*b4*(t-b2)))));
der.b3=1/(1+exp(2-4*b4*(t-b2)));
der.b4= (4*b3*(t-b2)*exp(2-4*b4*(t-b2)))/((1+exp(2-4*b4*(t-b2)))*(1+exp(2-4*b4*(t-b2))));
BY BOTELLA;
RUN;
Valora esta pregunta
Me gusta: Está pregunta es útil y esta claraNo me gusta: Está pregunta no esta clara o no es útil
0
Responder