errorObs {SCEPtER} | R Documentation |
Perturbe a set of observational data
Description
The function returns a perturbed object, starting from observational data and a vector of uncertainties.
Usage
errorObs(sigma, STAR, parallel=FALSE, corr=0)
Arguments
sigma |
a vector of 7 elements containing the uncertainties on observational data. The error on seismic parameters must be expressed as percentage. |
STAR |
a vector of 7 elements containing the observational data. |
parallel |
logical flag specifying if the computation should be performed in parallel on a multi-core architecture. |
corr |
the value of assumed the correlation between asteroseismic parameters. Default to zero (no correlation). |
Details
The vector of observations STAR
must contain, in the order:
The effective temperature (in K) of the star.
The logarithm of the surface gravity of the star (in cm s^-2).
The metallicity [Fe/H] of the star.
The value of large frequency separation of the star, divided by the solar value of this quantity.
The value of the frequency of maximum oscillation power of the star, divided by the solar value of this quantity.
The mass (in solar units) of the star (if known).
The radius (in solar unit) of the star (if known).
The vector of the uncertainties on the observation sigma
must
contain the uncertainty on the above quantity. The
uncertainty on the last
two values must be expressed as relative error.
Value
errorObs
returns an object obtained sampling from a
multivariate normal distribution with vector of mean star
and
covariance matrix diag(sigma^2)
, with elements (4,5) and (5,4)
modified according to the value of corr
.