next up previous
Next: About this document

UNSW Astrophysics Postgraduate Course
Module 3, Session 2
Optical and Infrared Techniques
Lecture 5: Data Reduction and Analysis
Tutorial Answers
http://www.phys.unsw.edu.au/ tex2html_wrap_inline61 spd/pgc_astro/part3.html

  1. Lecture Notes
  2. Tutorial Questions

1. (a) Binomial distribution, so expected number of failures tex2html_wrap_inline63 .

(b) Probability that x students fail = tex2html_wrap_inline65 .

So P( tex2html_wrap_inline67 fail) = 1 - P(0 fail) - P(1 fail) - P(2 fail) - P(3 fail) - P(4 fail)
= 1 - 0.0884 - 0.2228 - 0.270 - 0.2142 - 0.1223 = 0.0803.

ie 8% chance of 5 or more failing.

2. (a) Poisson distribution tex2html_wrap_inline69 .

P(detecting 8 or more events in a day) = 1 - tex2html_wrap_inline71 P(i events seen)
= 1 - 0.1353 - 0.2707 - 0.2707 - 0.1804 - 0.0902 - 0.0301 - 0.0120 - 0.0034 = 0.0011.

(b) Considering just a 10 minute period then tex2html_wrap_inline73 events expected. Do above sum using this value of tex2html_wrap_inline75 and you get 0 to about 20 decimal places!

3. (a) Binomial distribution with p= P(R) = 670/1000 and q = P(L) = 330/1000. So uncertainty is tex2html_wrap_inline77 . Applies to both tex2html_wrap_inline79 and tex2html_wrap_inline81 .

(b) tex2html_wrap_inline83 . tex2html_wrap_inline85 . So tex2html_wrap_inline87 .

(c) Assume A=0.400. So tex2html_wrap_inline89 . Thus 0.4 = p - (1-p) so p = 0.7 and q = 1-p = 0.3. So tex2html_wrap_inline97 . Then tex2html_wrap_inline99 .

4. The frequency shift tex2html_wrap_inline101 .
So tex2html_wrap_inline103 , assuming the variables are independent.

We find that tex2html_wrap_inline105 , tex2html_wrap_inline107 and tex2html_wrap_inline109 . So that tex2html_wrap_inline111 . So tex2html_wrap_inline113 . Note that if we were measuring the actual frequency, instead of calculating what its shift would be, then the uncertainty in its measurement is greater than the actual shift in frequency!

5. Gaussian distribution. Let tex2html_wrap_inline115 and look for probability that tex2html_wrap_inline117 . We seek tex2html_wrap_inline119 . From tables we find that tex2html_wrap_inline121 so that tex2html_wrap_inline123 , respectively.

6. Bin into 10 mark ranges. Expected number per bin is tex2html_wrap_inline125 where tex2html_wrap_inline127 is bin width and N total number of observations, with tex2html_wrap_inline131 . Events within each bin are distributed according to the Poisson distribution as we are extracting a random sample from a parent population for that bin value. For the Poisson distribution tex2html_wrap_inline133 as tex2html_wrap_inline135 . So the error estimate per bin is simply tex2html_wrap_inline137 .

Thus we can calculate tex2html_wrap_inline139 , the observed number per bin, tex2html_wrap_inline141 , the expected number per bin and tex2html_wrap_inline143 , their spread for that bin. We find that tex2html_wrap_inline145 . We have 10 bins so 9 degrees of freedom. ie tex2html_wrap_inline147 . From tables tex2html_wrap_inline149 and tex2html_wrap_inline151 . Thus in repeated experiments the probability of obtaining tex2html_wrap_inline153 . This is a high probability, so the data is consistent with coming from a Gaussian parent population with tex2html_wrap_inline155 and tex2html_wrap_inline157 .

Anyone who answered 29 or 38 to the final part is failed for impertinence!

image not found or accessible

Worksheet for calculating the number of samples per bin, the expected number per bin, the uncertainity and the Chi-square statisitic from the deviations.

image not found or accessible

Comparison between binned data and best Gaussian fit, whose parameters are determined from the mean and standard deviation of the data. The Chi-square test says the data are consistent with having come from this parent distribution.



next up previous
Next: About this document

Michael Burton
Wed Sep 3 22:13:59 EST 1997