453819_en-US

liliconcom

need script that computes the statistical parameters shown below in my post

Hi - I am still new at jsl scripting and am trying to write a script that creates and saves four new columns based on my linear regression fit:

- lower 95% mean confidence interval

- upper 95% mean confidence interval

- lower 95% individual confidence interval

- upper 95% individual confidence interval

Finally, each of the above are computed in terms of log 10 and I would instead like them converted to antilog terms (e.g. My x-y fit produces 2.78, but I need to show the antilog of 2.78 which is 2489).  A portion of my current jsl code is shown below.

Thank you!

obj = Fit Model(
 Y( :LOG_ECOLI_PEAR ),
 Effects( :LOG_TURB_PEAR_FNU),
 Personality( "Standard Least Squares" ),
 Emphasis( "Effect Leverage" ),
 Run(:LOG_ECOLI_PEAR 
 )) ;
 
  obj << Prediction Formula;

 

ScriptingRe: need script that computes the statistical parameters shown below in my post

Thank you for the response Mark.  I am working through the solution now.  As I am still learning JSL, much of what I do is trila and error so it may take me a while to sort through this.  In any case, I have gotten a portion of it to work and will circle back once I am successful or have a follow up question.  Thanks!

Re: need script that computes the statistical parameters shown below in my post

If you do not need the platform for interaction, then you might use the Minimal Report for Emphasis and the Invisible option in the Fit Model call. (Then remember to close it.

 

I suggest that you use the Scripting Index to explore the other messages for Fit Least Squares platform instead of << Prediction Formula. Here is one example:

 

 

The column formulas are actually JSL expressions. As such, they may be manipulated as data. That is to say, you can apply the anti-log transformation after these expressions are saved to result in an updated formula. So use the << Get Formula message with the new column to obtain the current expression. Insert it into the anti-log function expression as its argument,  and then save the new expression with the << Set Formula message to replace the original formula.

Re: need script that computes the statistical parameters shown below in my post

I got the confidence intervals to work and I am now trying to get the anti log (i.e.back transform) portion working.  Can you briefly explain the following lines of code?  (I think it is looking for a certain column in my data table for back transforming?)

 

n = N Col( dt );
For( col = n, col > n - 3, col--, // start at end and work back