View source: R/TG_density.fd.R
TG_density.fd | R Documentation |
Like the regular S-PLUS function density
, this function
computes a probability density function for a sample of values of a
random variable. However, in this case the density function is
defined by a functional parameter object logdensfdPar
along with a
normalizing constant C
.
The density function $p(indexdens)$ has the form
p(indexdens) = C exp[W(indexdens)]
where function $W(indexdens)$ is defined by the functional data object
logdensfdPar
.
## S3 method for class 'fd'
TG_density(indexdens, logdensfd, conv=0.0001, iterlim=20,
active=1:nbasis, dbglev=0)
indexdens |
a set observations, which may be one of two forms:
The first option corresponds to all $f_i = 1$. |
logdensfd |
a functional data object specifying the initial value, basis object, roughness penalty and smoothing parameter defining function $W(t).$ |
conv |
a positive constant defining the convergence criterion. |
iterlim |
the maximum number of iterations allowed. |
active |
a logical vector of length equal to the number of coefficients
defining |
dbglev |
either 0, 1, or 2. This controls the amount information printed out on each iteration, with 0 implying no output, 1 intermediate output level, and 2 full output. If levels 1 and 2 are used, it is helpful to turn off the output buffering option in S-PLUS. |
The goal of the function is provide a smooth density function
estimate that approaches some target density by an amount that is
controlled by the linear differential operator Lfdobj
and
the penalty parameter. For example, if the second derivative of
$W(t)$ is penalized heavily, this will force the function to
approach a straight line, which in turn will force the density function
itself to be nearly normal or Gaussian. Similarly, to each textbook
density function there corresponds a $W(t)$, and to each of these
in turn their corresponds a linear differential operator that will, when
apply to $W(t)$, produce zero as a result.
To plot the density function or to evaluate it, evaluate Wfdobj
,
exponentiate the resulting vector, and then divide by the normalizing
constant C
.
a named list of length 4 containing:
Wfdobj |
a functional data object defining function $W(indexdens)$ that that optimizes the fit to the data of the monotone function that it defines. |
C |
the normalizing constant. |
Flist |
a named list containing three results for the final converged solution: (1) f: the optimal function value being minimized, (2) grad: the gradient vector at the optimal solution, and (3) norm: the norm of the gradient vector at the optimal solution. |
iternum |
the number of iterations. |
iterhist |
a |
Juan Li and James Ramsay
Ramsay, J. O., Li J. and Wiberg, M. (2020) Full information optimal scoring. Journal of Educational and Behavioral Statistics, 45, 297-315.
Ramsay, J. O., Li J. and Wiberg, M. (2020) Better rating scale scores with information-based psychometrics. Psych, 2, 347-360.
intensity.fd
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.