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)) }