Difference between revisions of "SMHS Estimation"

From SOCR
Jump to: navigation, search
(Scientific Methods for Health Sciences - Parameter Estimation)
Line 71: Line 71:
 
MSE(δ)=var(δ)+(bias(δ))^2. The MVUE minimizes MSE among unbiased estimators. In some cases biased estimators have lower MSE because they have a smaller variance than does any unbiased estimator.  
 
MSE(δ)=var(δ)+(bias(δ))^2. The MVUE minimizes MSE among unbiased estimators. In some cases biased estimators have lower MSE because they have a smaller variance than does any unbiased estimator.  
 
   
 
   
 +
 
4) Applications
 
4) Applications
  
Line 84: Line 85:
 
http://socr.ucla.edu/htmls/SOCR_Distributions.html  
 
http://socr.ucla.edu/htmls/SOCR_Distributions.html  
 
http://socr.ucla.edu/htmls/SOCR_Experiments.html http://socr.ucla.edu/htmls/SOCR_Charts.html  
 
http://socr.ucla.edu/htmls/SOCR_Experiments.html http://socr.ucla.edu/htmls/SOCR_Charts.html  
 +
 +
  
 
6) Problems
 
6) Problems
Line 131: Line 134:
 
c. Mean = 11, SE = 2 


 
c. Mean = 11, SE = 2 


 
d. Mean = 13, SE = 3 

 
d. Mean = 13, SE = 3 

 +
  
  
Line 137: Line 141:
 
http://mirlyn.lib.umich.edu/Record/002457053  
 
http://mirlyn.lib.umich.edu/Record/002457053  
 
http://mirlyn.lib.umich.edu/Record/005059626  
 
http://mirlyn.lib.umich.edu/Record/005059626  
 +
  
  

Revision as of 10:25, 24 July 2014

Scientific Methods for Health Sciences - Parameter Estimation

IV. HS 850: Fundamentals

Estimation

1) Overview: Estimation is an important concept in the field of statistics and application of estimation is widely applied in various areas. It deals with estimating values of parameters of the population based on the sample data. And the parameters describe an underlying physical setting and their value would affect the distribution of the measured data. Two major approaches are commonly used in estimation: (1) the probabilistic approach assumes that the measured data is random with probability distribution dependent on the parameters; (2) the set-membership approach assumes that the measured data vector belongs to a set which depends on the parameter vector. The purpose of estimation is to find an estimator that is interpretable, accurate and exhibits some form of optimality. Indicators like minimum variance unbiased estimator is usually applied to measure estimator optimality, although it is possible that an optimal estimator don’t always exist. Here we present the fundamentals of estimation theory and illustrate how to apply estimation in real studies.


2) Motivation: To obtain a desired estimator, or estimation, we need to first determine a probability distribution with parameters of interest based on the data. After deciding the probabilistic model, we need to find the theoretically achievable precision available to any estimator based on the model and then develop an estimator based on this model. There are variety of methods and criteria to develop and choose between estimators based on their performance: maximum likelihood estimators, Bayes estimators, method of moments estimators, minimum mean square error estimators, minimum variance unbiased estimator, best linear unbiased estimator, etc. Experiment or simulations can also be run to test estimators’ performance.


3) Theory

3.1) An estimate of a population parameter may be expressed in two ways:

  • Point estimate: a single value of estimate. For example, sample mean is a point estimate of the population mean.
  • Interval estimate: an interval estimate is defined by two numbers, between which a population parameter is said to lie.

3.2) Confidence Intervals: describes the uncertainly of a sampling method and contains a confidence level, a statistic and a margin of error. The statistic and the margin of error define an interval estimate, which represent the precision of the method. Confidence Interval is expressed as sample statistic ± margin of error. Interpretation of a confidence interval at 95% confidence level is that we have 95% confidence that the parameter will fall within the margin of the interval.

3.3) Confidence level: the probability part of a confidence interval. It describes the likelihood that a particular sampling method will produce a confidence interval that includes the true population parameter.

3.4) Margin of error: range of the values above and below the sample statistic in confidence interval. margin of error=critical value*standard deviation of the statistic.

3.5) Critical value: The central limit theorem states that the sampling distribution of a statistic will be normal or nearly normal and the critical value can be expressed as a t score or as a z score, if ANY of the following conditions apply:

  • The population distribution is normal;
  • The sampling distribution is symmetric, unimodal, without outliers, and the sample size is 15 or less;
  • The sampling distribution is moderately skewed, unimodal, without outliers, and the sample size is between 16 and 40;
  • The sample size is greater than 40, without outliers.

To find the critical value, follow these steps.

  • Compute alpha (α): α = 1 - (confidence level / 100)
  • Find the critical probability (p*): p* = 1 - α/2
  • To express the critical value as a z score, find the z score having a cumulative probability equal to the critical probability (p*).
  • To express the critical value as a t score, follow these steps. Find the degree of freedom (DF): when estimating a mean score or a proportion from a single sample, DF is equal to the sample size minus one. For other applications, the degrees of freedom may be calculated differently. We will describe those computations as they come up.

The critical t score (t*) is the t score having degrees of freedom equal to DF and a cumulative probability equal to the critical probability (p*).

Should you express the critical value as a t score or as a z score? As a practical matter, when the sample size is large (greater than 40), it doesn't make much difference. Both approaches yield similar results. Strictly speaking, when the population standard deviation is unknown or when the sample size is small, the t score is preferred. Nevertheless, many introductory statistics texts use the z score exclusively.

3.6) Standard error: an estimate of the standard deviation of a statistic. When the values of population parameters are unknown, it is valuable to compute the standard error as an unbiased estimate of the standard deviation of a statistic. It is computed form known sample statistic. The table below shows how to compute the standard error for simple random samples assuming that the population size is at least 10 times larger than the sample size.

Sample proportion, p || 〖SE〗_p=sqrt(p(1-p)/n)
Statistic Standard error
Sample mean, x ̅ 〖SE〗_x ̅ =s/sqrt(n)
Difference between means, (x_1 ) ̅-(x_2 ) ̅ 〖SE〗_( (x_1 ) ̅-(x_2 ) ̅ )=sqrt(s_1^2/n_1+s_2^2/n_2)
Difference between proportions, (p_1 ) ̅-(p_2 ) ̅ 〖SE〗_( (p_1 ) ̅-(p_2 ) ̅ )=sqrt(p_1 (1-p_1)/n_1+p_2 (1-p_2)/n_2)


3.7) Degrees of freedom: the number of independent pieces of information on which the estimate is based

In general, the degrees of freedom for an estimate is equal to the number of values minus the number of parameters estimated to the estimate in question. Suppose we have sampled 20 data points then our estimate of the variance has 20 – 1 = 19 degree of freedom.

3.8) Characteristics of Estimators: Bias: refers to whether an estimator tends to either overestimate or underestimate the parameter. We say an estimator is biased if the mean of the sampling distribution of the statistic is not equal to the parameter. For example, σ^2=(x-μ)^2/N is a biased estimator of the population variance and sample variance s^2=(x-x ̅ )^2/(N-1) is unbiased estimate of the population variance. Sampling variability: refers to how much the estimate varies from sample to sample. It is usually measured by its standard error: the smaller the standard error, the less the sampling variability. For example, the standard error of the mean is σ_M=σ/√N. So the larger the sample size (N), the smaller the standard error of the mean, hence the smaller the sample variability. Unbiased estimate: η(X_1,X_2,…,X_n )=E[δ(X_1,X_2,…,X_n )|T] then δ(X_1,X_2,…,X_n ) is unbiased estimate for g(θ) and T is a complete sufficient statistic for the family of densities. (Uniformly) Minimum-variance unbiased estimator (UMVUE, or MVUE) is an unbiased estimator that has lower variance than any other unbiased estimator for all possible values of the parameter. It may not exist. Consider estimation of g(θ) based on data X_1,X_2,…,X_n independent and identically distributed from some member of a family with density p_θ, θ ϵ Ω , an unbiased estimator δ(X_1,X_2,…,X_n) of g(θ) is UMVUE if ∀θ ϵ Ω, var(δ(X_1,X_2,…,X_n ))≤var(δ ̃(X_1,X_2,…,X_n )) for any other unbiased estimator δ ̃. MSE(δ)=var(δ)+(bias(δ))^2. The MVUE minimizes MSE among unbiased estimators. In some cases biased estimators have lower MSE because they have a smaller variance than does any unbiased estimator.


4) Applications

4.1) This article (http://wiki.stat.ucla.edu/socr/index.php/AP_Statistics_Curriculum_2007_Estim_MOM_MLE) presents the MOM and MLE methods of estimation. It illustrates the MOM method in detailed examples and attached several exercise for students to practice. MOM, which is short for Method Of Moments, is one of the most commonly used methods to estimate population parameters using observed data from the specific process. The idea is to use the sample data to calculate sample moments and then set these equal to their corresponding population counterparts. Steps: (1) determine the k parameters of interest and specific distribution for this process; (2) compute the first k (or more) sample-moments; (3) set the sample-moments equal to the population moments and solve for a system of k equations with k unknowns. Let’s look at a simple example as application of the MOM method. Consider we want to estimate the true probability of a head by flipping the coins (assume a unfair coin). Suppose we flip the coin 10 times and observe the following outcome: {H,T,H,H,T,T,T,H,T,T}. With MOM: (1) the parameter of interest is p=P(H) and it follows a Bernoulli distribution, (2) np=E[Y]=4,p=2/5, where Y is the number of heads for one experiment and it follows a Binomial distribution. (3) estimate of true probability of flipping a head in one experiment equals 2/5. This is an easy example of MOM proportion example.

4.1) This article (http://onlinestatbook.com/2/estimation/estimation.html) presents a fundamental introduction to estimation theory and illustrated on basic concepts and application of estimation. It offers specific examples and exercises on each concept and application and works as a good start of introduction to estimation theory.

4.2) This article (http://digital-library.theiet.org/content/journals/10.1049/ip-f-2.1993.0015) proposed an algorithm, the bootstrap filter, for implementing recursive Bayesian filters. The required density of the state vector is represented as a set of random samples, which are updated and propagated by the algorithm. The method presented is not restricted by assumptions of linearity or Gaussian noise and it may be applied to any state transition or measurement model. It presents a simulation example of the bearings only tracking problems and includes schemes for improving the efficiency of the basic algorithm.


5) Software http://socr.ucla.edu/htmls/SOCR_Distributions.html http://socr.ucla.edu/htmls/SOCR_Experiments.html http://socr.ucla.edu/htmls/SOCR_Charts.html


6) Problems

6.1) Which of the following statements is true. a. When the margin of error is small, the confidence level is high. 
 b. When the margin of error is small, the confidence level is low. 
 c. A confidence interval is a type of point estimate. 
 d. A population mean is an example of a point estimate. e. None of the above.

6.2) Which of the following statements is true. a. The standard error is computed solely from sample attributes. 
 b. The standard deviation is computed solely from sample attributes. 
 c. The standard error is a measure of central tendency. d. All of the above. e. None of the above.

6.3) 900 students were randomly selected for a national survey. Among survey participants, the mean grade-point average (GPA) was 2.7, and the standard deviation was 0.4. What is the margin of error, assuming a 95% confidence level? a. 0.013 
 b. 0.025 
 c. 0.500 
 d. 1.960 


6.4) Suppose we want to estimate the average weight of an adult male in Dekalb County, Georgia. We draw a random sample of 1,000 men from a population of 1,000,000 men and weigh them. We find that the average man in our sample weighs 180 pounds, and the standard deviation of the sample is 30 pounds. What is the 95% confidence interval? a. 180±1.86 
 b. 180±3.0 
 c. 180±5.88 d. 180±30

6.5) Suppose that simple random samples of seniors are selected from two colleges: 15 students from school A and 20 students from school B. On a standardized test, the sample from school A has an average score of 1000 with a standard deviation of 100. The sample from school B has an average score of 950 with a standard deviation of 90. What is the 90% confidence interval for the difference in test scores at the two schools, assuming that test scores came from normal distributions in both schools? (Hint: Since the sample sizes are small, use a t score as the critical value.) a. 50 + 1.70 
 b. 50 + 28.49 
 c. 50 + 32.74 
 d. 50 + 55.66

6.6) You know the population mean for a certain test score. You select 10 people from the population to estimate the standard deviation. How many degrees of freedom does your estimation of the standard deviation have? a. 8 b. 9 c. 10 d. 11

6.6) In the population, a parameter has a value of 10. Based on the means and standard errors of their sampling distributions, which of these statistics estimates this parameter with the least sampling variability? 

 a. Mean = 10, SE = 5 

 b. Mean = 9, SE = 4

 c. Mean = 11, SE = 2 

 d. Mean = 13, SE = 3 



7) References http://mirlyn.lib.umich.edu/Record/004199238 http://mirlyn.lib.umich.edu/Record/002457053 http://mirlyn.lib.umich.edu/Record/005059626


Answers: e, a, b, a, d, c, c





Translate this page:

(default)
Uk flag.gif

Deutsch
De flag.gif

Español
Es flag.gif

Français
Fr flag.gif

Italiano
It flag.gif

Português
Pt flag.gif

日本語
Jp flag.gif

България
Bg flag.gif

الامارات العربية المتحدة
Ae flag.gif

Suomi
Fi flag.gif

इस भाषा में
In flag.gif

Norge
No flag.png

한국어
Kr flag.gif

中文
Cn flag.gif

繁体中文
Cn flag.gif

Русский
Ru flag.gif

Nederlands
Nl flag.gif

Ελληνικά
Gr flag.gif

Hrvatska
Hr flag.gif

Česká republika
Cz flag.gif

Danmark
Dk flag.gif

Polska
Pl flag.png

România
Ro flag.png

Sverige
Se flag.gif