cancel
Showing results for 
Show  only  | Search instead for 
Did you mean: 
JMP is taking Discovery online, April 16 and 18. Register today and join us for interactive sessions featuring popular presentation topics, networking, and discussions with the experts.
Choose Language Hide Translation Bar
Mickyboy
Level V

95% confidence intervals

Hi All,

 

l am trying to get a summary from a data table, l would like the mean of a variable l10t, the standard error, and 95% confidence interval. It all seemed fairly straight forward, but am having trouble with the 95% confidence intervals, have tried different ways to no avail, currently have 

dt2<< Summary(
	Group( :Sample ),
	Mean( :dilution factor, :l10t ),
	Freq( "None" ),
	Weight( "None" ),
	Output Table( "log10t" ),
	Std Err( :l10t ),
	Confidence Interval( :l10t, 0.95),
	Link to original data table(0)

l am obviously doing something incorrectly, can anyone offer some advice.

 

Thanks Heaps

4 REPLIES 4
mzwald
Staff

Re: 95% confidence intervals

You can accomplish something similar like this:

Names Default To Here( 1 );
dt = Open( "$SAMPLE_DATA/Big Class.jmp" );
platform = dt << Oneway(
	Y( :height ),
	X( :sex ),
	Means and Std Dev( 1 ),
	Box Plots( 1 ),
	Mean Diamonds( 1 ),
	Mean Error Bars( 1 ),
	Std Dev Lines( 1 ),
	Connect Means( 1 ),
	Points Jittered( 1 ),
	Grand Mean( 0 )
);
Wait( 0 );
Report( platform )[Outline Box( "Oneway Analysis of height By sex" )][
Outline Box( "Means and Std Deviations" )][Table Box( 1 )] << Make Into Data Table;
Report( platform ) << Close Window;

Place your group variable in the X() role.

Mickyboy
Level V

Re: 95% confidence intervals

Hi mzwald,

 

Thanks very much for your help and your suggestion, i tried this code

platform = dt2 << Oneway(
	Y( :l10t ),
	X( :fold dilution ),
	Means and Std Dev( 1 ),
	Box Plots( 1 ),
	Mean Diamonds( 1 ),
	Mean Error Bars( 1 ),
	Std Dev Lines( 1 ),
	Connect Means( 1 ),
	Points Jittered( 1 ),
	Grand Mean( 0 )
);
Wait( 0 );
Report( platform )[Outline Box( "Oneway Analysis of height By sex" )][
Outline Box( "Means and Std Deviations" )][Table Box( 1 )] << Make Into Data Table;
Report( platform ) << Close Window;

and received the following error message

Cannot subscript Display Box in access or evaluation of 'Report(platform)[Outline Box("Oneway Analysis of height By sex")]' , Report( platform )[/*###*/Outline Box( "Oneway Analysis of height By sex" )]

however, l am almost there as to what i want to achieve for my data table with this,

dt2<< Summary(
	Group( :Sample ),
	Mean( :dilution factor, :l10t ),
	Freq( "None" ),
	Weight( "None" ),
	Output Table( "log10t" ),
	Std Err( :l10t ),
	Confidence Interval( :Name ("Mean(l10t)")),
	Link to original data table(0)
);

Loq = Data Table( "log10t" );

Loq << New Column( "Back T'sform Mean", Numeric, Continuous, Formula( Round(10^:Name ("Mean(l10t)"),0)));

Loq << New Column( "Back T'sform Std error", Numeric, Continuous, Formula(10^:Name ("Std Err(l10t)"), Format( "fixed dec", 6, 0 ) ) );

giving me this,

Data Table.PNG so l am almost there, l was hoping the confidence intervals would be a single simple line of code.

 

Thanks

Mickyboy

mzwald
Staff

Re: 95% confidence intervals

Well it looks like Summary doesn't currently support CI's, but you can always make a formula for them yourself.
The formulas for a 95% confidence interval (assuming a normal distribution) would be:
Upper95CI = mean + 1.96*stdev
Lower95CI = mean - 1.96*stdev
You can generate those two columns from the fields returned by the Summary platform. Those formulas will get you pretty close to what you are looking for.

Mickyboy
Level V

Re: 95% confidence intervals

Hi mzwald,

 

i have used the following and has worked a treat

New Column( "Lower 95% CI", Numeric, Continuous, Formula ( Round(:Name ("Mean(l10t)") - (:Name("Std Err(l10t)") * t Quantile (0.975,23)), 4)));

if i have used the following to create a variable

dt2<< Summary(
	Group( :Sample ),
	Mean( :dilution factor, :l10t ),
	Freq( "None" ),
	Weight( "None" ),
	Output Table( "log10t" ),
	Std Err( :l10t ),
	Link to original data table(0)
);

and i would like to round the variable Mean(l10t), to four decimal places, can you do that from the above block of syntax???

 

thanks

Mickyboy