| class.numer {MiscPsycho} | R Documentation |
~~ A concise (1-5 lines) description of what the function does. ~~
class.numer(x, b, prof_cut, mu = 0, sigma = 1, aboveQ)
x |
~~Describe x here~~ |
b |
~~Describe b here~~ |
prof_cut |
~~Describe prof_cut here~~ |
mu |
~~Describe mu here~~ |
sigma |
~~Describe sigma here~~ |
aboveQ |
~~Describe aboveQ here~~ |
~~ If necessary, more details than the description above ~~
~Describe the value returned If it is a LIST, use
comp1 |
Description of 'comp1' |
comp2 |
Description of 'comp2' |
...
....
~~further notes~~
~Make other sections like Warning with section{Warning }{....} ~
~~who you are~~
~put references to the literature/web site here ~
~~objects to See Also as help, ~~~
##---- Should be DIRECTLY executable !! ----
##-- ==> Define data, use random,
##-- or do help(data=index) for the standard data sets.
## The function is currently defined as
function(x,b, prof_cut, mu=0, sigma=1, aboveQ){
gauss_numer <- gauss.quad(49,kind="laguerre")
if(aboveQ==FALSE){
mat <- rbind(like.mat(x,b, (prof_cut-gauss_numer$nodes)), dnorm(prof_cut-gauss_numer$nodes, mean=mu, sd=sigma))
} else { mat <- rbind(like.mat(x,b, (gauss_numer$nodes+prof_cut)), dnorm(gauss_numer$nodes+prof_cut, mean=mu, sd=sigma))
}
f_y <- rbind(apply(mat, 2, prod), exp(gauss_numer$nodes), gauss_numer$weights)
sum(apply(f_y,2,prod))
}