Stat 5810 Applied Spatial Statistics                                 Fall 2000


                 Introduction to S and S-Plus - Answers
                 ======================================

           Wed 9/13/00, Mon 9/18/00, Wed 9/20/00 & Mon 9/25/00


Dr. Juergen Symanzik

------------------------------------------------------------------------------

Q1:

> quest1 _ 0:100 * 2 * pi / 100
> quest1
  [1] 0.00000000 0.06283185 0.12566371 0.18849556 0.25132741 0.31415927
  [7] 0.37699112 0.43982297 0.50265482 0.56548668 0.62831853 0.69115038
 [13] 0.75398224 0.81681409 0.87964594 0.94247780 1.00530965 1.06814150
 [19] 1.13097336 1.19380521 1.25663706 1.31946891 1.38230077 1.44513262
 [25] 1.50796447 1.57079633 1.63362818 1.69646003 1.75929189 1.82212374
 [31] 1.88495559 1.94778745 2.01061930 2.07345115 2.13628300 2.19911486
 [37] 2.26194671 2.32477856 2.38761042 2.45044227 2.51327412 2.57610598
 [43] 2.63893783 2.70176968 2.76460154 2.82743339 2.89026524 2.95309709
 [49] 3.01592895 3.07876080 3.14159265 3.20442451 3.26725636 3.33008821
 [55] 3.39292007 3.45575192 3.51858377 3.58141563 3.64424748 3.70707933
 [61] 3.76991118 3.83274304 3.89557489 3.95840674 4.02123860 4.08407045
 [67] 4.14690230 4.20973416 4.27256601 4.33539786 4.39822972 4.46106157
 [73] 4.52389342 4.58672527 4.64955713 4.71238898 4.77522083 4.83805269
 [79] 4.90088454 4.96371639 5.02654825 5.08938010 5.15221195 5.21504380
 [85] 5.27787566 5.34070751 5.40353936 5.46637122 5.52920307 5.59203492
 [91] 5.65486678 5.71769863 5.78053048 5.84336234 5.90619419 5.96902604
 [97] 6.03185789 6.09468975 6.15752160 6.22035345 6.28318531
> plot (quest1, cos (quest1), type = "o")


Q2:

> seq (-40, 100, by = 10)
 [1]  -40  -30  -20  -10    0   10   20   30   40   50   60   70   80   90  100
> (seq (-40, 100, by = 10) - 32) / 1.8
 [1] -40.000000 -34.444444 -28.888889 -23.333333 -17.777778 -12.222222
 [7]  -6.666667  -1.111111   4.444444  10.000000  15.555556  21.111111
[13]  26.666667  32.222222  37.777778


Q3:

par (mfrow = c(2,2), oma = c (0, 0, 4, 0))
plot (sort (rbeta(100, 1, 1)), qbeta( (1:100 - 0.5) / 100, 1, 1), main = "Beta")
plot (sort (rf(100, 5, 8)), qf( (1:100 - 0.5) / 100, 5, 8), main = "F")
plot (sort (rt(100, 5)), qt( (1:100 - 0.5) / 100, 5), main = "t")
plot (sort (rgamma(100, 5)), qgamma( (1:100 - 0.5) / 100, 5), main = "Gamma")
mtext (side = 3, line = 0, cex = 2, outer = TRUE, "Several Distributions")