내 실험 데이터에 비선형 모델 적합을보고하는 가장 좋은 방법을 찾으려고합니다.
필자는 비선형 모델에 맞는 스크립트를 작성하여 손실 함수를 사용하여 매개 변수 추정치를 찾습니다. 그런 다음 각 매개 변수 추정치와 비선형 오류를보고합니다.
그러나 결과 fit의 RMSE 또는 R-square도 얻고 싶습니다. 비선형 모델에서는 사용할 수 없습니까? 또는 추가 분석 단계를 추가해야합니까?
비선형 해석 플랫폼을 대화식으로 실행할 때 결과 플롯은 실험 데이터와 모델 커브 모두를 보여줍니다. 그러나 스크립트를 자동으로 실행하면 실험 데이터 만 표시됩니다. 왜 이런 일이 일어나는 지 아는 사람이 있습니까? 아래 스크립트에서 무엇인가 놓치고 있습니까?
누군가 자료 / 그래프 / 책 / 링크를 제공하여 데이터 그래프를 얻고 결과 모델 적합성을 확인하고 적합성의 R-square를보고 할 수 있습니까?
비선형 모델 적합을 위해 사용하고있는 스크립트는 다음과 같습니다.
col1 = New Column( "Model" );
col1 << set formula(
Parameter(
{Y1_Y1_0 = 0.2, K1_K1_0 = 0.01, P1_P1_0 = 0.004, Y2_Y2_0 = 0.15, K2_K2_0 = 0.06, P2_P2_0 = 0.3},
If( :Days > 0,
Y1_Y1_0 * Exp( -Exp( ((K1_K1_0 * Exp()) / Y1_Y1_0) * (P1_P1_0 - :Days) ) ) + Y2_Y2_0 *
Exp( -Exp( ((K2_K2_0 * Exp()) / Y2_Y2_0) * (P2_P2_0 - :Days) ) ),
0
)
)
);
col2 = New Column( "Err" );
col2 << set Formula( Sum( (:Average - :Model) ^ 2 ) );
col2 << Evalformula;
col2 << Get formula;
obj = Nonlinear( Y( :Methane Yield ), X( :Model ), Loss( :Err ), Iteration Limit( 600 ), By( :Sample ) );
obj << go;
//returns a list of vectors where each vector, contains the parameters estimates
obj << get estimates;
//also see Get SSE, Get Parameter Names, Get CI
//Alternately, make a table
_xx = obj << Xpath( "//OutlineBox[@helpKey='Nonlinear Solution']" );
//only 2 groups so only 2 are found; each contains two sub tables
fit_dt = (_xx[1] << Find( Table Box( 1 ) )) << Make Combined Data Table;
est_dt = (_xx[1] << Find( Table Box( 2 ) )) << Make Combined Data Table;
fit_dt << set Name( "NonLin Fit Errors" );
est_dt << set Name( "NonLin Fit Estimates" );
고마워.
내가 제안한대로 손실 기능을 변경했습니다. 피팅을 실행했지만 여전히 플롯에서 모델 커브를 볼 수 없습니다.
내 솔루션에서 나는 평균 손실과 Sqrt 평균 손실을 얻습니다. 이것이 나의 RSME가 그 다음의 무한한 손실이라는 것을 의미합니까?
RMSE는 솔루션 개요에보고됩니다.
적합 모델과 축소 모델의 두 모델이 필요하기 때문에 R 스퀘어를 사용할 수 없습니다. 비선형 모델의 경우 감소 된 모델이 없습니다. (OLS 회귀 분석에서 축소 된 모델은 상수 절편 만있는 모델입니다.)
Nonlinear에서 기본 손실 함수를 잘못 정의하고 있습니다. 기본 손실은 최소 제곱 오류입니다. 스크립트에서 최소 제곱 오류를 사용하려고합니다. 열 수식은 전체 표본이 아닌 각 관측치의 손실을 가져와야합니다. 따라서 손실 공식은 제곱 차이 만 가져야합니다. (즉, 제곱 된 차이의 합을 계산하지 마십시오.) 기본 손실 대신이 방법을 사용하면 얻을 수 있습니다.
따라서 사용자 지정 손실을 사용하면 솔루션 개요의 정보가 전체 손실 (이 경우 SSE), 평균 손실 (이 경우 MSE) 및 Sqrt 평균 손실 (이 경우 RMSE)으로 변경됩니다.
Nonlinear에서 스크립트를 저장하고 실행했습니다. 데이터 플롯은 여전히 비선형 함수의 그래프를 나타내는 곡선을 보여줍니다.
