235574_es-ES

liliconcom

Intervalos de confianza del 95%

Hola a todos,

 

Estoy tratando de obtener un resumen de una tabla de datos, me gustaría la media de una variable l10t, el error estándar y el intervalo de confianza del 95%. Todo parecía bastante sencillo, pero estoy teniendo problemas con los intervalos de confianza del 95%, he intentado diferentes formas en vano, actualmente tengo 

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)

Obviamente estoy haciendo algo incorrectamente, ¿alguien puede ofrecer algún consejo?

 

Montones de gracias

Re: intervalos de confianza del 95%

Puedes lograr algo similar a esto:

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;

Coloque su variable de grupo en el rol X ().

Re: intervalos de confianza del 95%

Hola mzwald

 

Muchas gracias por su ayuda y su sugerencia, probé este código

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;

y recibí el siguiente mensaje de error

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" )]

Sin embargo, estoy casi allí en cuanto a lo que quiero lograr para mi tabla de datos con esto,

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 ) ) );

dándome esto

 así que estoy casi allí, esperaba que los intervalos de confianza fueran una sola línea de código simple.

 

Gracias

Mickyboy

Re: intervalos de confianza del 95%

Bueno, parece que Summary no admite actualmente CI's, pero siempre puedes hacer una fórmula para ellos.
Las fórmulas para un intervalo de confianza del 95% (suponiendo una distribución normal) serían:
Upper95CI = media + 1.96 * stdev
Lower95CI = media - 1.96 * stdev
Puede generar esas dos columnas a partir de los campos devueltos por la plataforma Resumen. Esas fórmulas lo acercarán bastante a lo que está buscando.

Re: intervalos de confianza del 95%

Hola mzwald

 

He usado lo siguiente y he trabajado de maravilla

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

si he usado lo siguiente para crear una 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)
);

y me gustaría redondear la variable Media (l10t) a cuatro decimales, ¿puedes hacer eso desde el bloque de sintaxis anterior?

 

Gracias

Mickyboy