Displaying 4 results from an estimated 4 matches for "aboveq".
Did you mean:
above
2006 Jul 20
2
Timing benefits of mapply() vs. for loop was: Wrap a loop inside a function
...at[i, ] <- pcm(theta, items[[i]],
score[[i]])
like.mat
}
# turn this off for now
#like.mat <- function(score, items, theta){
#matrix(mapply(pcm,rep(theta,length(items)),items,score),ncol=length(the
ta),byrow=TRUE)
#}
class.numer <- function(score,items, prof_cut, mu=0, sigma=1, aboveQ){
gauss_numer <- gauss.quad(49,kind="laguerre")
if(aboveQ==FALSE){
mat <- rbind(like.mat(score,items, (prof_cut-gauss_numer$nodes)),
dnorm(prof_cut-gauss_numer$nodes, mean=mu, sd=sigma))
} else { mat <- rbind(like.mat(score,items,
(gauss_numer$nodes+prof_cut)...
2006 May 05
0
Spline integration & Gaussian quadrature (was: gauss.quad.prob)
...any suggestions would be very much appreciated.
library(statmod)
rasch <- function(b,theta){
1 / (1 + exp(outer(b,theta,'-')))
}
like.mat <- function(x,b,theta){
rasch(b, theta)^x * (1 - rasch(b,theta))^(1-x)
}
class.numer <- 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_...
2007 May 08
0
MiscPsycho Package 1.0
...meters
scoreCon(aa$params)
The student in row 1 has a number correct score of 3 out of 10. From the
score conversion table, his theta estimate is -1.36. So, we can ask,
what is the probability that his true score is above, say, .5 on the
theta scale.
class.acc(as.numeric(xx$d[1,]), aa$params, .5, aboveQ=T)
Any recommendations, bugs, hints, comments, are welcome.
Harold Doran
[[alternative HTML version deleted]]
_______________________________________________
R-packages mailing list
R-packages at stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-packages
2007 May 08
0
MiscPsycho Package 1.0
...meters
scoreCon(aa$params)
The student in row 1 has a number correct score of 3 out of 10. From the
score conversion table, his theta estimate is -1.36. So, we can ask,
what is the probability that his true score is above, say, .5 on the
theta scale.
class.acc(as.numeric(xx$d[1,]), aa$params, .5, aboveQ=T)
Any recommendations, bugs, hints, comments, are welcome.
Harold Doran
[[alternative HTML version deleted]]
_______________________________________________
R-packages mailing list
R-packages at stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-packages