lib/statsample/srs.rb in statsample-0.5.1 vs lib/statsample/srs.rb in statsample-0.6.0

- old
+ new

@@ -33,17 +33,20 @@ # Sample size estimation for proportions, finite poblation. def estimation_n(d,prop,n_pobl,margin=0.95) n0=estimation_n0(d,prop,margin) n0.quo( 1 + ((n0 - 1).quo(n_pobl))) end + + # Proportion confidence interval with t values # Uses estimated proportion, sample without replacement. def proportion_confidence_interval_t(prop, n_sample, n_population, margin=0.95) t = Distribution::T.p_value(1-((1-margin).quo(2)) , n_sample-1) proportion_confidence_interval(prop,n_sample,n_population, t) end + # Proportion confidence interval with z values # Uses estimated proportion, sample without replacement. def proportion_confidence_interval_z(p, n_sample, n_population, margin=0.95) z=Distribution::Normal.p_value(1-((1-margin).quo(2))) proportion_confidence_interval(p,n_sample,n_population, z) @@ -51,10 +54,10 @@ # Proportion confidence interval with x value # Uses estimated proportion, sample without replacement def proportion_confidence_interval(p, sam,pop , x) f=sam.quo(pop) - one_range=x * Math::sqrt((qf(sam, pop) * p * (1-p)) / (sam-1)) + (1.quo(sam * 2.0)) + one_range=x * Math::sqrt((qf(sam, pop) * p * (1-p)).quo (sam-1)) + (1.quo(sam * 2.0)) [p-one_range, p+one_range] end # Standard deviation for sample distribution of a proportion # Know proportion, sample with replacement. # Based on http://stattrek.com/Lesson6/SRS.aspx