lab2 - DATASETS NEEDED coalash.dt(Coal Ash dataset...

Info icon This preview shows pages 1–3. Sign up to view the full content.

View Full Document Right Arrow Icon
Sheet1 Page 1 DATASETS NEEDED: coalash.dt (Coal Ash dataset) www4.stat.ncsu.edu/~fuentes/coalash.dat davis.dat (topographic heights) www4.stat.ncsu.edu/~fuentes/davis.dat ozone.dat (ozone values) www4.stat.ncsu.edu/~fuentes/ozone.dat software R with library geoR, fields and akima PART II Spatial prediction: ##EXERCISE 1: Estimating covariance parameters ##EXERCISE 2: Moving neighborhood kriging ##EXERCISE 3: ANISTROPY ##EXERCISE 4: Maps and Images ####################################################################### # download your data using read.table # to learn more about read.table type # > ? read.table coal.ash<-read.table('coalash.dat') coal.m<-as.matrix(coal.ash) davis.dat<-read.table('davis.dat') davis.m<-as.matrix(davis.dat) ################################### ##EXERCISE 1: Estimating model parameters ###A. Bayesian estimates of model parameters: #Function: krige.bayes
Image of page 1

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full Document Right Arrow Icon
Sheet1 Page 2 # You obtain posteriors for the trend, range, partial sill and nugget # Recommended (default) priors: # uniform for beta, 1/sigma^2 for the partial sill # and discrete uniform for the range and nugget. # # aniso.pars are the parameters for anisotropy: stretching and rotation angle. # # the default is with nugget 0, but you # can also get a posterior for the nugget when you give a prior to the nugget, # by saying nugget.prior="uniform" # #Additional information: #---------------------- # This function computes (Bayesian) estimates of a spatial linear #model and performs Bayesian and/or kriging prediction in a set of #locations specified by the user. # Priors options for mean and/or covariance parameters # coords : data coordinates (vector for 1d or matrix for 2d data) # data : vector with data values # locations : coordinates of points to be estimated (vector for 1d or # matrix for 2d data). If locations='no' only model parameters estimates are #returned if the full bayesian model is considered. # trend.d : trend data in data locations. Default is constant trend. #The options '1st' or '2nd' builts a first or second degree polinomial.
Image of page 2
Image of page 3
This is the end of the preview. Sign up to access the rest of the document.

{[ snackBarMessage ]}

What students are saying

  • Left Quote Icon

    As a current student on this bumpy collegiate pathway, I stumbled upon Course Hero, where I can find study resources for nearly all my courses, get online help from tutors 24/7, and even share my old projects, papers, and lecture notes with other students.

    Student Picture

    Kiran Temple University Fox School of Business ‘17, Course Hero Intern

  • Left Quote Icon

    I cannot even describe how much Course Hero helped me this summer. It’s truly become something I can always rely on and help me. In the end, I was not only able to survive summer classes, but I was able to thrive thanks to Course Hero.

    Student Picture

    Dana University of Pennsylvania ‘17, Course Hero Intern

  • Left Quote Icon

    The ability to access any university’s resources through Course Hero proved invaluable in my case. I was behind on Tulane coursework and actually used UCLA’s materials to help me move forward and get everything together on time.

    Student Picture

    Jill Tulane University ‘16, Course Hero Intern