[rkward-cvs] SF.net SVN: rkward:[2717] trunk/rkward
m-eik at users.sourceforge.net
m-eik at users.sourceforge.net
Mon Nov 30 20:46:34 UTC 2009
Revision: 2717
http://rkward.svn.sourceforge.net/rkward/?rev=2717&view=rev
Author: m-eik
Date: 2009-11-30 20:46:11 +0000 (Mon, 30 Nov 2009)
Log Message:
-----------
irt: added gpcm estimation, re-worked a lot of the html output
Modified Paths:
--------------
trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_2pl.php
trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_2pl.rkh
trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_3pl.php
trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_3pl.rkh
trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_lltm.php
trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_lltm.rkh
trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_rasch.php
trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_rasch.rkh
trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_grm.php
trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_grm.rkh
trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_lpcm.php
trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_lpcm.rkh
trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_lrsm.php
trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_lrsm.rkh
trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_pcm.php
trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_pcm.rkh
trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_rsm.php
trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_rsm.rkh
trunk/rkward/rkward/plugins/analysis/irt/tests/ltm_unidimensional.php
trunk/rkward/rkward/plugins/analysis/irt/tests/ltm_unidimensional.xml
trunk/rkward/rkward/plugins/irt.pluginmap
trunk/rkward/tests/item_response_theory/RKTestStandard.2PL_parameter_estimation.rkcommands.R
trunk/rkward/tests/item_response_theory/RKTestStandard.2PL_parameter_estimation.rkout
trunk/rkward/tests/item_response_theory/RKTestStandard.3PL_parameter_estimation.rkcommands.R
trunk/rkward/tests/item_response_theory/RKTestStandard.3PL_parameter_estimation.rkout
trunk/rkward/tests/item_response_theory/RKTestStandard.GRM_parameter_estimation.rkcommands.R
trunk/rkward/tests/item_response_theory/RKTestStandard.GRM_parameter_estimation.rkout
trunk/rkward/tests/item_response_theory/RKTestStandard.LLTM_parameter_estimation.rkcommands.R
trunk/rkward/tests/item_response_theory/RKTestStandard.LLTM_parameter_estimation.rkout
trunk/rkward/tests/item_response_theory/RKTestStandard.LPCM_parameter_estimation.rkcommands.R
trunk/rkward/tests/item_response_theory/RKTestStandard.LPCM_parameter_estimation.rkout
trunk/rkward/tests/item_response_theory/RKTestStandard.LRSM_parameter_estimation.rkcommands.R
trunk/rkward/tests/item_response_theory/RKTestStandard.LRSM_parameter_estimation.rkout
trunk/rkward/tests/item_response_theory/RKTestStandard.PCM_parameter_estimation.rkcommands.R
trunk/rkward/tests/item_response_theory/RKTestStandard.PCM_parameter_estimation.rkout
trunk/rkward/tests/item_response_theory/RKTestStandard.RSM_parameter_estimation.rkcommands.R
trunk/rkward/tests/item_response_theory/RKTestStandard.RSM_parameter_estimation.rkout
trunk/rkward/tests/item_response_theory/RKTestStandard.Rasch_parameter_estimation.rkcommands.R
trunk/rkward/tests/item_response_theory/RKTestStandard.Rasch_parameter_estimation.rkout
trunk/rkward/tests/item_response_theory/RKTestStandard.goodnes-of-fit_Rasch.rkout
trunk/rkward/tests/item_response_theory/RKTestStandard.unidimensional.rkcommands.R
trunk/rkward/tests/item_response_theory/RKTestStandard.unidimensional.rkout
trunk/rkward/tests/item_response_theory.R
Modified: trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_2pl.php
===================================================================
--- trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_2pl.php 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_2pl.php 2009-11-30 20:46:11 UTC (rev 2717)
@@ -53,7 +53,10 @@
$save_name = getRK_val("save_name");
?>
rk.header ("2PL parameter estimation")
-rk.print (estimates.2pl)
+rk.print (paste("Call: <code>",deparse(estimates.2pl$call, width.cutoff=500),"</code>"))
+rk.print ("<h4>Coefficients:</h4>")
+rk.print (coef(estimates.2pl))
+rk.print (paste("Log-likelihood value at convergence:",round(estimates.2pl$log.Lik, digits=1)))
<?
// check if results are to be saved:
if ($save && $save_name) {
Modified: trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_2pl.rkh
===================================================================
--- trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_2pl.rkh 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_2pl.rkh 2009-11-30 20:46:11 UTC (rev 2717)
@@ -7,7 +7,7 @@
<usage>
To fit the two parameter logistic model with default settings, just the matrix of manifest variables (dichotomous item responses) must be selected. Further options are described in the sections below.
- The complete results are stored as a new object in the current workspace for further use (e.g. plotting item characteristic curves). This object will be called <code>estimates.2pl</code>.
+ The complete results can be stored as a new object in the current workspace for further use (e.g. plotting item characteristic curves).
</usage>
<settings>
Modified: trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_3pl.php
===================================================================
--- trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_3pl.php 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_3pl.php 2009-11-30 20:46:11 UTC (rev 2717)
@@ -59,7 +59,10 @@
$save_name = getRK_val("save_name");
?>
rk.header ("3PL parameter estimation")
-rk.print (estimates.3pl)
+rk.print (paste("Call: <code>",deparse(estimates.3pl$call, width.cutoff=500),"</code>"))
+rk.print ("<h4>Coefficients:</h4>")
+rk.print (coef(estimates.3pl))
+rk.print (paste("Log-likelihood value at convergence:",round(estimates.3pl$log.Lik, digits=1)))
<?
// check if results are to be saved:
if ($save && $save_name) {
Modified: trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_3pl.rkh
===================================================================
--- trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_3pl.rkh 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_3pl.rkh 2009-11-30 20:46:11 UTC (rev 2717)
@@ -7,7 +7,7 @@
<usage>
To fit the Birnbaum three parameter model with default settings, just the matrix of manifest variables (dichotomous item responses) must be selected. Further options are described in the sections below.
- The complete results are stored as a new object in the current workspace for further use (e.g. plotting item characteristic curves). This object will be called <code>estimates.3pl</code>.
+ The complete results can be stored as a new object in the current workspace for further use (e.g. plotting item characteristic curves).
</usage>
<settings>
Modified: trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_lltm.php
===================================================================
--- trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_lltm.php 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_lltm.php 2009-11-30 20:46:11 UTC (rev 2717)
@@ -34,7 +34,11 @@
$save_name = getRK_val("save_name");
?>
rk.header ("LLTM parameter estimation")
-rk.print (estimates.lltm)
+rk.print (paste("Call: <code>",deparse(estimates.lltm$call, width.cutoff=500),"</code>"))
+rk.print ("<h4>Coefficients:</h4>")
+rk.print(t(rbind(Eta=estimates.lltm$etapar,StdErr=estimates.lltm$se.eta)))
+rk.print (paste("Conditional log-likelihood:",round(estimates.lltm$loglik, digits=1),
+"<br />Number of iterations:",estimates.lltm$iter,"<br />Number of parameters:",estimates.lltm$npar))
<?
// check if results are to be saved:
if ($save && $save_name) {
Modified: trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_lltm.rkh
===================================================================
--- trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_lltm.rkh 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_lltm.rkh 2009-11-30 20:46:11 UTC (rev 2717)
@@ -7,7 +7,7 @@
<usage>
To fit the linear logstic test model with default settings, just the matrix of manifest variables (dichotomous item responses) must be selected. Further options are described in the sections below. For detailed information on the expected layout of matrices and vectors please refer to the manual for the <link href="rkward://rhelp/LLTM">LLTM()</link> function.
- The complete results are stored as a new object in the current workspace for further use (e.g. plotting item characteristic curves). This object will be called <code>estimates.lltm</code>.
+ The complete results can be stored as a new object in the current workspace for further use (e.g. plotting item characteristic curves).
</usage>
<settings>
Modified: trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_rasch.php
===================================================================
--- trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_rasch.php 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_rasch.php 2009-11-30 20:46:11 UTC (rev 2717)
@@ -46,9 +46,13 @@
// check whether parameter estimations should be kept in the global enviroment
$save = getRK_val("chk_save");
$save_name = getRK_val("save_name");
+ $irtparam = getRK_val("irtparam");
?>
rk.header ("Rasch parameter estimation")
-rk.print (estimates.rasch$coefficients)
+rk.print (paste("Call: <code>",deparse(estimates.rasch$call, width.cutoff=500),"</code>"))
+rk.print ("<h4>Coefficients:</h4>")
+rk.print (coef(estimates.rasch))
+rk.print (paste("Log-likelihood value at convergence:",round(estimates.rasch$log.Lik, digits=1)))
<?
// check if results are to be saved:
if ($save && $save_name) {
Modified: trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_rasch.rkh
===================================================================
--- trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_rasch.rkh 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/rkward/plugins/analysis/irt/dichotomous/par_est_rasch.rkh 2009-11-30 20:46:11 UTC (rev 2717)
@@ -7,7 +7,7 @@
<usage>
To fit the Rasch model with default settings, just the matrix of manifest variables (dichotomous item responses) must be selected. Further options are described in the sections below.
- The complete results are stored as a new object in the current workspace for further use (e.g. plotting item characteristic curves). This object will be called <code>estimates.rasch</code>.
+ The complete results can be stored as a new object in the current workspace for further use (e.g. plotting item characteristic curves).
</usage>
<settings>
Modified: trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_grm.php
===================================================================
--- trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_grm.php 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_grm.php 2009-11-30 20:46:11 UTC (rev 2717)
@@ -53,7 +53,10 @@
$save_name = getRK_val("save_name");
?>
rk.header ("GRM parameter estimation")
-rk.print (estimates.grm)
+rk.print (paste("Call: <code>",deparse(estimates.grm$call, width.cutoff=500),"</code>"))
+rk.print ("<h4>Coefficients:</h4>")
+rk.print (coef(estimates.grm))
+rk.print (paste("Log-likelihood value at convergence:",round(estimates.grm$log.Lik, digits=1)))
<?
// check if results are to be saved:
if ($save && $save_name) {
Modified: trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_grm.rkh
===================================================================
--- trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_grm.rkh 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_grm.rkh 2009-11-30 20:46:11 UTC (rev 2717)
@@ -7,7 +7,7 @@
<usage>
To fit the graded response model with default settings, just the matrix of manifest variables (polytomous item responses) must be selected. Further options are described in the sections below.
- The complete results are stored as a new object in the current workspace for further use (e.g. plotting item characteristic curves). This object will be called <code>estimates.grm</code>.
+ The complete results can be stored as a new object in the current workspace for further use (e.g. plotting item characteristic curves).
</usage>
<settings>
Modified: trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_lpcm.php
===================================================================
--- trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_lpcm.php 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_lpcm.php 2009-11-30 20:46:11 UTC (rev 2717)
@@ -34,7 +34,11 @@
$save_name = getRK_val("save_name");
?>
rk.header ("LPCM parameter estimation")
-rk.print (estimates.lpcm)
+rk.print (paste("Call: <code>",deparse(estimates.lpcm$call, width.cutoff=500),"</code>"))
+rk.print ("<h4>Coefficients:</h4>")
+rk.print(t(rbind(Eta=estimates.lpcm$etapar,StdErr=estimates.lpcm$se.eta)))
+rk.print (paste("Conditional log-likelihood:",round(estimates.lpcm$loglik, digits=1),
+"<br />Number of iterations:",estimates.lpcm$iter,"<br />Number of parameters:",estimates.lpcm$npar))
<?
// check if results are to be saved:
if ($save && $save_name) {
Modified: trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_lpcm.rkh
===================================================================
--- trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_lpcm.rkh 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_lpcm.rkh 2009-11-30 20:46:11 UTC (rev 2717)
@@ -7,7 +7,7 @@
<usage>
To fit the linear partial credit model with default settings, just the matrix of manifest variables (item responses) must be selected. Further options are described in the sections below. For detailed information on the expected layout of matrices and vectors please refer to the manual for the <link href="rkward://rhelp/LPCM">LPCM()</link> function.
- The complete results are stored as a new object in the current workspace for further use (e.g. plotting item characteristic curves). This object will be called <code>estimates.lpcm</code>.
+ The complete results can be stored as a new object in the current workspace for further use (e.g. plotting item characteristic curves).
</usage>
<settings>
Modified: trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_lrsm.php
===================================================================
--- trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_lrsm.php 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_lrsm.php 2009-11-30 20:46:11 UTC (rev 2717)
@@ -34,7 +34,11 @@
$save_name = getRK_val("save_name");
?>
rk.header ("LRSM parameter estimation")
-rk.print (estimates.lrsm)
+rk.print (paste("Call: <code>",deparse(estimates.lrsm$call, width.cutoff=500),"</code>"))
+rk.print ("<h4>Coefficients:</h4>")
+rk.print(t(rbind(Eta=estimates.lrsm$etapar,StdErr=estimates.lrsm$se.eta)))
+rk.print (paste("Conditional log-likelihood:",round(estimates.lrsm$loglik, digits=1),
+"<br />Number of iterations:",estimates.lrsm$iter,"<br />Number of parameters:",estimates.lrsm$npar))
<?
// check if results are to be saved:
if ($save && $save_name) {
Modified: trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_lrsm.rkh
===================================================================
--- trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_lrsm.rkh 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_lrsm.rkh 2009-11-30 20:46:11 UTC (rev 2717)
@@ -7,7 +7,7 @@
<usage>
To fit the linear rating scale model with default settings, just the matrix of manifest variables (item responses) must be selected. Further options are described in the sections below. For detailed information on the expected layout of matrices and vectors please refer to the manual for the <link href="rkward://rhelp/LRSM">LRSM()</link> function.
- The complete results are stored as a new object in the current workspace for further use (e.g. plotting item characteristic curves). This object will be called <code>estimates.lrsm</code>.
+ The complete results can be stored as a new object in the current workspace for further use (e.g. plotting item characteristic curves).
</usage>
<settings>
Modified: trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_pcm.php
===================================================================
--- trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_pcm.php 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_pcm.php 2009-11-30 20:46:11 UTC (rev 2717)
@@ -29,7 +29,11 @@
$save_name = getRK_val("save_name");
?>
rk.header ("PCM parameter estimation")
-rk.print (estimates.pcm)
+rk.print (paste("Call: <code>",deparse(estimates.pcm$call, width.cutoff=500),"</code>"))
+rk.print ("<h4>Coefficients:</h4>")
+rk.print(t(rbind(Eta=estimates.pcm$etapar,StdErr=estimates.pcm$se.eta)))
+rk.print (paste("Conditional log-likelihood:",round(estimates.pcm$loglik, digits=1),
+"<br />Number of iterations:",estimates.pcm$iter,"<br />Number of parameters:",estimates.pcm$npar))
<?
// check if results are to be saved:
if ($save && $save_name) {
Modified: trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_pcm.rkh
===================================================================
--- trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_pcm.rkh 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_pcm.rkh 2009-11-30 20:46:11 UTC (rev 2717)
@@ -7,7 +7,7 @@
<usage>
To fit the partial credit model with default settings, just the matrix of manifest variables (item responses) must be selected. Further options are described in the sections below. For detailed information on the expected layout of matrices and vectors please refer to the manual for the <link href="rkward://rhelp/PCM">PCM()</link> function.
- The complete results are stored as a new object in the current workspace for further use (e.g. plotting item characteristic curves). This object will be called <code>estimates.pcm</code>.
+ The complete results can be are stored as a new object in the current workspace for further use (e.g. plotting item characteristic curves).
</usage>
<settings>
Modified: trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_rsm.php
===================================================================
--- trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_rsm.php 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_rsm.php 2009-11-30 20:46:11 UTC (rev 2717)
@@ -29,7 +29,11 @@
$save_name = getRK_val("save_name");
?>
rk.header ("RSM parameter estimation")
-rk.print (estimates.rsm)
+rk.print (paste("Call: <code>",deparse(estimates.rsm$call, width.cutoff=500),"</code>"))
+rk.print ("<h4>Coefficients:</h4>")
+rk.print(t(rbind(Eta=estimates.rsm$etapar,StdErr=estimates.rsm$se.eta)))
+rk.print (paste("Conditional log-likelihood:",round(estimates.rsm$loglik, digits=1),
+"<br />Number of iterations:",estimates.rsm$iter,"<br />Number of parameters:",estimates.rsm$npar))
<?
// check if results are to be saved:
if ($save && $save_name) {
Modified: trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_rsm.rkh
===================================================================
--- trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_rsm.rkh 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/rkward/plugins/analysis/irt/polytomous/par_est_rsm.rkh 2009-11-30 20:46:11 UTC (rev 2717)
@@ -7,7 +7,7 @@
<usage>
To fit the rating scale model with default settings, just the matrix of manifest variables (item responses) must be selected. Further options are described in the sections below. For detailed information on the expected layout of matrices and vectors please refer to the manual for the <link href="rkward://rhelp/RSM">RSM()</link> function.
- The complete results are stored as a new object in the current workspace for further use (e.g. plotting item characteristic curves). This object will be called <code>estimates.rsm</code>.
+ The complete results can be stored as a new object in the current workspace for further use (e.g. plotting item characteristic curves).
</usage>
<settings>
Modified: trunk/rkward/rkward/plugins/analysis/irt/tests/ltm_unidimensional.php
===================================================================
--- trunk/rkward/rkward/plugins/analysis/irt/tests/ltm_unidimensional.php 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/rkward/plugins/analysis/irt/tests/ltm_unidimensional.php 2009-11-30 20:46:11 UTC (rev 2717)
@@ -9,16 +9,24 @@
// let's read all values into php variables for the sake of readable code
$spin_samples = getRK_val("spin_samples");
-?>unidim.res <<- unidimTest(<? getRK("x");
+?>unidim.res <- unidimTest(<? getRK("x");
// check if any options must be inserted
if($spin_samples != "100") echo(", B=$spin_samples") ;
?>)
<?}
function printout () {
+ $save = getRK_val("chk_save");
+ $save_name = getRK_val("save_name");
?>
rk.header ("Unidimensionality check (<? getRK("x"); ?>)")
rk.print (unidim.res)
<?
+// check if results are to be saved:
+if ($save && $save_name) {
+?>
+# keep results in current workspace
+<? echo($save_name); ?> <<- unidim.res
+<?}
}
?>
\ No newline at end of file
Modified: trunk/rkward/rkward/plugins/analysis/irt/tests/ltm_unidimensional.xml
===================================================================
--- trunk/rkward/rkward/plugins/analysis/irt/tests/ltm_unidimensional.xml 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/rkward/plugins/analysis/irt/tests/ltm_unidimensional.xml 2009-11-30 20:46:11 UTC (rev 2717)
@@ -4,6 +4,8 @@
<!--<help file="ltm_unidimensional.rkh" />-->
<logic>
+ <connect client="save_name.enabled" governor="chk_save.state" />
+ <connect client="save_name.required" governor="chk_save.state" />
</logic>
<dialog label="Unidimensionality check (Rasch, LTM, 3PM)">
@@ -18,6 +20,10 @@
<frame id="frame_options" label="Monte Carlo options">
<spinbox id="spin_samples" label="Number of samples:" type="integer" initial="100" min="1" size="small" />
</frame>
+ <frame id="frame_save">
+ <checkbox id="chk_save" label="Keep results in current workspace" value="save" checked="true"/>
+ <saveobject id="save_name" initial="unidim.res"/>
+ </frame>
</column>
</row>
</tab>
Modified: trunk/rkward/rkward/plugins/irt.pluginmap
===================================================================
--- trunk/rkward/rkward/plugins/irt.pluginmap 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/rkward/plugins/irt.pluginmap 2009-11-30 20:46:11 UTC (rev 2717)
@@ -10,19 +10,20 @@
<component type="standard" id="par_est_grm" file="analysis/irt/polytomous/par_est_grm.xml" label="Graded response model fit" />
<component type="standard" id="par_est_rsm" file="analysis/irt/polytomous/par_est_rsm.xml" label="Rating scale model fit" />
<component type="standard" id="par_est_pcm" file="analysis/irt/polytomous/par_est_pcm.xml" label="Partial credit model fit" />
+ <component type="standard" id="par_est_gpcm" file="analysis/irt/polytomous/par_est_gpcm.xml" label="Generalized partial credit model fit" />
<component type="standard" id="par_est_lrsm" file="analysis/irt/polytomous/par_est_lrsm.xml" label="Linear rating scale model fit" />
<component type="standard" id="par_est_lpcm" file="analysis/irt/polytomous/par_est_lpcm.xml" label="Linear partial credit model fit" />
- <component type="standard" id="ltm_cronbach_alpha" file="analysis/irt/tests/ltm_cronbach_alpha.xml" label="Cronbach's alpha" />
+ <component type="standard" id="ltm_cronbach_alpha" file="analysis/irt/tests/ltm_cronbach_alpha.xml" label="Cronbach's alpha" />
- <component type="standard" id="ltm_gof_rasch" file="analysis/irt/tests/ltm_gof_rasch.xml" label="Goodness of Fit (Rasch)" />
- <component type="standard" id="ltm_unidimensional" file="analysis/irt/tests/ltm_unidimensional.xml" label="Unidimensionality check (Rasch, LTM, 3PM)" />
- <component type="standard" id="ltm_item_fit" file="analysis/irt/tests/ltm_item_fit.xml" label="Item-fit statistics (Rasch, LTM, 3PM)" />
+ <component type="standard" id="ltm_gof_rasch" file="analysis/irt/tests/ltm_gof_rasch.xml" label="Goodness of Fit (Rasch)" />
+ <component type="standard" id="ltm_unidimensional" file="analysis/irt/tests/ltm_unidimensional.xml" label="Unidimensionality check (Rasch, LTM, 3PM)" />
+ <component type="standard" id="ltm_item_fit" file="analysis/irt/tests/ltm_item_fit.xml" label="Item-fit statistics (Rasch, LTM, 3PM)" />
<component type="standard" id="ltm_person_fit" file="analysis/irt/tests/ltm_person_fit.xml" label="Person-fit statistics (Rasch, LTM, 3PM)" />
<component type="standard" id="eRm_waldtest" file="analysis/irt/tests/eRm_waldtest.xml" label="Wald test (RSM, PCM)" />
<component type="standard" id="eRm_plotLR" file="plots/irt/tests/eRm_plotLR.xml" label="Andersen LR Plot (RSM, PCM)" />
- <component type="standard" id="plot_rasch" file="plots/irt/dichotomous/plot_rasch.xml" label="Plot fitted Rasch model" />
+ <component type="standard" id="plot_rasch" file="plots/irt/dichotomous/plot_rasch.xml" label="Plot fitted Rasch model" />
<component type="standard" id="plot_ltm" file="plots/irt/dichotomous/plot_ltm.xml" label="Plot fitted two parameter logistic model" />
<component type="standard" id="plot_tpm" file="plots/irt/dichotomous/plot_tpm.xml" label="Plot fitted Birnbaum three parameter model" />
@@ -50,6 +51,7 @@
<entry component="par_est_grm" />
<entry component="par_est_rsm" />
<entry component="par_est_pcm" />
+ <entry component="par_est_gpcm" />
<entry component="par_est_lrsm" />
<entry component="par_est_lpcm" />
</menu>
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.2PL_parameter_estimation.rkcommands.R
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.2PL_parameter_estimation.rkcommands.R 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.2PL_parameter_estimation.rkcommands.R 2009-11-30 20:46:11 UTC (rev 2717)
@@ -5,7 +5,10 @@
estimates.2pl <- ltm(WIRS ~ z1 * z2)
## Print result
rk.header ("2PL parameter estimation")
-rk.print (estimates.2pl)
+rk.print (paste("Call: <code>",deparse(estimates.2pl$call, width.cutoff=500),"</code>"))
+rk.print ("<h4>Coefficients:</h4>")
+rk.print (coef(estimates.2pl))
+rk.print (paste("Log-likelihood value at convergence:",round(estimates.2pl$log.Lik, digits=1)))
# keep results in current workspace
estimates.2pl <<- estimates.2pl
})
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.2PL_parameter_estimation.rkout
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.2PL_parameter_estimation.rkout 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.2PL_parameter_estimation.rkout 2009-11-30 20:46:11 UTC (rev 2717)
@@ -1,3 +1,19 @@
<h1>2PL parameter estimation</h1>
DATE<br>
-<p class='character'><br>Call:<br>ltm(formula = WIRS ~ z1 * z2)<br><br>Coefficients:<br> (Intercept) z1 z2 z1:z2<br>Item 1 -1.097 0.271 -3.317 -1.160<br>Item 2 0.944 1.091 2.173 2.955<br>Item 3 -1.458 1.811 -0.314 0.543<br>Item 4 -1.465 1.116 0.556 0.269<br>Item 5 -1.083 2.133 -0.466 1.147<br>Item 6 -2.753 1.698 -0.935 1.072<br><br>Log.Lik: -3293.2<br></p>
\ No newline at end of file
+<p class='character'>Call: <code> ltm(formula = WIRS ~ z1 * z2) </code></p><p class='character'><h4>Coefficients:</h4></p>
+<p align= center >
+<table cellspacing=0 border=1><caption align=bottom class=captiondataframe></caption><tr><td>
+ <table border=0 class=dataframe>
+ <tbody> <tr class= firstline > <th></th><th>(Intercept)</th><th>z1</th><th>z2</th><th>z1:z2</th> </tr>
+ <tr><td class=firstcolumn>Item 1</td><td class=cellinside>-1.10</td><td class=cellinside> 0.27</td><td class=cellinside>-3.32</td><td class=cellinside>-1.16</td></tr>
+ <tr><td class=firstcolumn>Item 2</td><td class=cellinside> 0.94</td><td class=cellinside> 1.09</td><td class=cellinside> 2.17</td><td class=cellinside> 2.95</td></tr>
+ <tr><td class=firstcolumn>Item 3</td><td class=cellinside>-1.46</td><td class=cellinside> 1.81</td><td class=cellinside>-0.31</td><td class=cellinside> 0.54</td></tr>
+ <tr><td class=firstcolumn>Item 4</td><td class=cellinside>-1.46</td><td class=cellinside> 1.12</td><td class=cellinside> 0.56</td><td class=cellinside> 0.27</td></tr>
+ <tr><td class=firstcolumn>Item 5</td><td class=cellinside>-1.08</td><td class=cellinside> 2.13</td><td class=cellinside>-0.47</td><td class=cellinside> 1.15</td></tr>
+ <tr><td class=firstcolumn>Item 6</td><td class=cellinside>-2.75</td><td class=cellinside> 1.70</td><td class=cellinside>-0.94</td><td class=cellinside> 1.07</td></tr>
+
+ </tbody>
+</table>
+ </td></table>
+ <br>
+<p class='character'>Log-likelihood value at convergence: -3293.2</p>
\ No newline at end of file
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.3PL_parameter_estimation.rkcommands.R
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.3PL_parameter_estimation.rkcommands.R 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.3PL_parameter_estimation.rkcommands.R 2009-11-30 20:46:11 UTC (rev 2717)
@@ -5,7 +5,10 @@
estimates.3pl <- tpm(LSAT)
## Print result
rk.header ("3PL parameter estimation")
-rk.print (estimates.3pl)
+rk.print (paste("Call: <code>",deparse(estimates.3pl$call, width.cutoff=500),"</code>"))
+rk.print ("<h4>Coefficients:</h4>")
+rk.print (coef(estimates.3pl))
+rk.print (paste("Log-likelihood value at convergence:",round(estimates.3pl$log.Lik, digits=1)))
# keep results in current workspace
estimates.3pl <<- estimates.3pl
})
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.3PL_parameter_estimation.rkout
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.3PL_parameter_estimation.rkout 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.3PL_parameter_estimation.rkout 2009-11-30 20:46:11 UTC (rev 2717)
@@ -1,3 +1,18 @@
<h1>3PL parameter estimation</h1>
DATE<br>
-<p class='character'><br>Call:<br>tpm(data = LSAT)<br><br>Coefficients:<br> Gussng Dffclt Dscrmn<br>Item 1 0.037 -3.296 0.829<br>Item 2 0.078 -1.145 0.760<br>Item 3 0.012 -0.249 0.902<br>Item 4 0.035 -1.766 0.701<br>Item 5 0.053 -2.990 0.666<br><br>Log.Lik: -2466.7<br></p>
\ No newline at end of file
+<p class='character'>Call: <code> tpm(data = LSAT) </code></p><p class='character'><h4>Coefficients:</h4></p>
+<p align= center >
+<table cellspacing=0 border=1><caption align=bottom class=captiondataframe></caption><tr><td>
+ <table border=0 class=dataframe>
+ <tbody> <tr class= firstline > <th></th><th>Gussng</th><th>Dffclt</th><th>Dscrmn</th> </tr>
+ <tr><td class=firstcolumn>Item 1</td><td class=cellinside> 0.037</td><td class=cellinside>-3.296</td><td class=cellinside> 0.829</td></tr>
+ <tr><td class=firstcolumn>Item 2</td><td class=cellinside> 0.078</td><td class=cellinside>-1.145</td><td class=cellinside> 0.760</td></tr>
+ <tr><td class=firstcolumn>Item 3</td><td class=cellinside> 0.012</td><td class=cellinside>-0.249</td><td class=cellinside> 0.902</td></tr>
+ <tr><td class=firstcolumn>Item 4</td><td class=cellinside> 0.035</td><td class=cellinside>-1.766</td><td class=cellinside> 0.701</td></tr>
+ <tr><td class=firstcolumn>Item 5</td><td class=cellinside> 0.053</td><td class=cellinside>-2.990</td><td class=cellinside> 0.666</td></tr>
+
+ </tbody>
+</table>
+ </td></table>
+ <br>
+<p class='character'>Log-likelihood value at convergence: -2466.7</p>
\ No newline at end of file
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.GRM_parameter_estimation.rkcommands.R
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.GRM_parameter_estimation.rkcommands.R 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.GRM_parameter_estimation.rkcommands.R 2009-11-30 20:46:11 UTC (rev 2717)
@@ -5,7 +5,10 @@
estimates.grm <- grm(Environment)
## Print result
rk.header ("GRM parameter estimation")
-rk.print (estimates.grm)
+rk.print (paste("Call: <code>",deparse(estimates.grm$call, width.cutoff=500),"</code>"))
+rk.print ("<h4>Coefficients:</h4>")
+rk.print (coef(estimates.grm))
+rk.print (paste("Log-likelihood value at convergence:",round(estimates.grm$log.Lik, digits=1)))
# keep results in current workspace
estimates.grm <<- estimates.grm
})
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.GRM_parameter_estimation.rkout
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.GRM_parameter_estimation.rkout 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.GRM_parameter_estimation.rkout 2009-11-30 20:46:11 UTC (rev 2717)
@@ -1,3 +1,19 @@
<h1>GRM parameter estimation</h1>
DATE<br>
-<p class='character'><br>Call:<br>grm(data = Environment)<br><br>Coefficients:<br> Extrmt1 Extrmt2 Dscrmn<br>LeadPetrol 0.487 2.584 1.378<br>RiverSea 1.058 2.499 2.341<br>RadioWaste 0.779 1.793 3.123<br>AirPollution 0.457 2.157 3.283<br>Chemicals 0.809 1.868 2.947<br>Nuclear 0.073 1.427 1.761<br><br>Log.Lik: -1090.4<br></p>
\ No newline at end of file
+<p class='character'>Call: <code> grm(data = Environment) </code></p><p class='character'><h4>Coefficients:</h4></p>
+<p align= center >
+<table cellspacing=0 border=1><caption align=bottom class=captiondataframe></caption><tr><td>
+ <table border=0 class=dataframe>
+ <tbody> <tr class= firstline > <th></th><th>Extrmt1</th><th>Extrmt2</th><th>Dscrmn</th> </tr>
+ <tr><td class=firstcolumn>LeadPetrol</td><td class=cellinside>0.487</td><td class=cellinside>2.584</td><td class=cellinside>1.378</td></tr>
+ <tr><td class=firstcolumn>RiverSea</td><td class=cellinside>1.058</td><td class=cellinside>2.499</td><td class=cellinside>2.341</td></tr>
+ <tr><td class=firstcolumn>RadioWaste</td><td class=cellinside>0.779</td><td class=cellinside>1.793</td><td class=cellinside>3.123</td></tr>
+ <tr><td class=firstcolumn>AirPollution</td><td class=cellinside>0.457</td><td class=cellinside>2.157</td><td class=cellinside>3.283</td></tr>
+ <tr><td class=firstcolumn>Chemicals</td><td class=cellinside>0.809</td><td class=cellinside>1.868</td><td class=cellinside>2.947</td></tr>
+ <tr><td class=firstcolumn>Nuclear</td><td class=cellinside>0.073</td><td class=cellinside>1.427</td><td class=cellinside>1.761</td></tr>
+
+ </tbody>
+</table>
+ </td></table>
+ <br>
+<p class='character'>Log-likelihood value at convergence: -1090.4</p>
\ No newline at end of file
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.LLTM_parameter_estimation.rkcommands.R
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.LLTM_parameter_estimation.rkcommands.R 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.LLTM_parameter_estimation.rkcommands.R 2009-11-30 20:46:11 UTC (rev 2717)
@@ -5,7 +5,11 @@
estimates.lltm <- LLTM(lltmdat1, mpoints=2)
## Print result
rk.header ("LLTM parameter estimation")
-rk.print (estimates.lltm)
+rk.print (paste("Call: <code>",deparse(estimates.lltm$call, width.cutoff=500),"</code>"))
+rk.print ("<h4>Coefficients:</h4>")
+rk.print(t(rbind(Eta=estimates.lltm$etapar,StdErr=estimates.lltm$se.eta)))
+rk.print (paste("Conditional log-likelihood:",round(estimates.lltm$loglik, digits=1),
+"<br />Number of iterations:",estimates.lltm$iter,"<br />Number of parameters:",estimates.lltm$npar))
# keep results in current workspace
estimates.lltm <<- estimates.lltm
})
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.LLTM_parameter_estimation.rkout
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.LLTM_parameter_estimation.rkout 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.LLTM_parameter_estimation.rkout 2009-11-30 20:46:11 UTC (rev 2717)
@@ -1,3 +1,28 @@
<h1>LLTM parameter estimation</h1>
DATE<br>
-<p class='character'><br>Results of LLTM estimation: <br><br>Call: LLTM(X = lltmdat1, mpoints = 2) <br><br>Conditional log-likelihood: -1511.2 <br>Number of iterations: 15 <br>Number of parameters: 15 <br><br>Basic Parameters eta:<br> eta 1 eta 2 eta 3 eta 4 eta 5 eta 6 eta 7 eta 8<br>Estimate 0.18215 0.32422 -0.70372 -0.24748 0.49126 0.81037 0.68577 0.25312<br>Std.Err 0.14900 0.14929 0.15799 0.15099 0.15025 0.15393 0.15220 0.14909<br> eta 9 eta 10 eta 11 eta 12 eta 13 eta 14 eta 15<br>Estimate -0.95333 0.25312 0.13485 -0.83996 -1.32511 0.51532 0.244881<br>Std.Err 0.16410 0.14909 0.14900 0.16111 0.17648 0.15044 0.082067<br><br></p>
\ No newline at end of file
+<p class='character'>Call: <code> LLTM(X = lltmdat1, mpoints = 2) </code></p><p class='character'><h4>Coefficients:</h4></p>
+<p align= center >
+<table cellspacing=0 border=1><caption align=bottom class=captiondataframe></caption><tr><td>
+ <table border=0 class=dataframe>
+ <tbody> <tr class= firstline > <th></th><th>Eta</th><th>StdErr</th> </tr>
+ <tr><td class=firstcolumn>eta 1</td><td class=cellinside> 0.182</td><td class=cellinside> 0.149</td></tr>
+ <tr><td class=firstcolumn>eta 2</td><td class=cellinside> 0.324</td><td class=cellinside> 0.149</td></tr>
+ <tr><td class=firstcolumn>eta 3</td><td class=cellinside>-0.704</td><td class=cellinside> 0.158</td></tr>
+ <tr><td class=firstcolumn>eta 4</td><td class=cellinside>-0.247</td><td class=cellinside> 0.151</td></tr>
+ <tr><td class=firstcolumn>eta 5</td><td class=cellinside> 0.491</td><td class=cellinside> 0.150</td></tr>
+ <tr><td class=firstcolumn>eta 6</td><td class=cellinside> 0.810</td><td class=cellinside> 0.154</td></tr>
+ <tr><td class=firstcolumn>eta 7</td><td class=cellinside> 0.686</td><td class=cellinside> 0.152</td></tr>
+ <tr><td class=firstcolumn>eta 8</td><td class=cellinside> 0.253</td><td class=cellinside> 0.149</td></tr>
+ <tr><td class=firstcolumn>eta 9</td><td class=cellinside>-0.953</td><td class=cellinside> 0.164</td></tr>
+ <tr><td class=firstcolumn>eta 10</td><td class=cellinside> 0.253</td><td class=cellinside> 0.149</td></tr>
+ <tr><td class=firstcolumn>eta 11</td><td class=cellinside> 0.135</td><td class=cellinside> 0.149</td></tr>
+ <tr><td class=firstcolumn>eta 12</td><td class=cellinside>-0.840</td><td class=cellinside> 0.161</td></tr>
+ <tr><td class=firstcolumn>eta 13</td><td class=cellinside>-1.325</td><td class=cellinside> 0.176</td></tr>
+ <tr><td class=firstcolumn>eta 14</td><td class=cellinside> 0.515</td><td class=cellinside> 0.150</td></tr>
+ <tr><td class=firstcolumn>eta 15</td><td class=cellinside> 0.245</td><td class=cellinside> 0.082</td></tr>
+
+ </tbody>
+</table>
+ </td></table>
+ <br>
+<p class='character'>Conditional log-likelihood: -1511.2 <br />Number of iterations: 15 <br />Number of parameters: 15</p>
\ No newline at end of file
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.LPCM_parameter_estimation.rkcommands.R
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.LPCM_parameter_estimation.rkcommands.R 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.LPCM_parameter_estimation.rkcommands.R 2009-11-30 20:46:11 UTC (rev 2717)
@@ -5,7 +5,11 @@
estimates.lpcm <- LPCM(lpcmdat, mpoints=2, groupvec=G)
## Print result
rk.header ("LPCM parameter estimation")
-rk.print (estimates.lpcm)
+rk.print (paste("Call: <code>",deparse(estimates.lpcm$call, width.cutoff=500),"</code>"))
+rk.print ("<h4>Coefficients:</h4>")
+rk.print(t(rbind(Eta=estimates.lpcm$etapar,StdErr=estimates.lpcm$se.eta)))
+rk.print (paste("Conditional log-likelihood:",round(estimates.lpcm$loglik, digits=1),
+"<br />Number of iterations:",estimates.lpcm$iter,"<br />Number of parameters:",estimates.lpcm$npar))
# keep results in current workspace
estimates.lpcm <<- estimates.lpcm
})
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.LPCM_parameter_estimation.rkout
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.LPCM_parameter_estimation.rkout 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.LPCM_parameter_estimation.rkout 2009-11-30 20:46:11 UTC (rev 2717)
@@ -1,3 +1,23 @@
<h1>LPCM parameter estimation</h1>
DATE<br>
-<p class='character'><br>Results of LPCM estimation: <br><br>Call: LPCM(X = lpcmdat, mpoints = 2, groupvec = G) <br><br>Conditional log-likelihood: -103.72 <br>Number of iterations: 28 <br>Number of parameters: 10 <br><br>Basic Parameters eta:<br> eta 1 eta 2 eta 3 eta 4 eta 5 eta 6 eta 7 eta 8<br>Estimate 0.41216 -0.29898 -0.13451 0.034894 -0.42887 -0.28640 0.16412 0.10074<br>Std.Err 0.39659 0.49423 0.40918 0.394328 0.47845 0.44882 0.39920 0.43511<br> eta 9 eta 10<br>Estimate -0.20149 1.09404<br>Std.Err 0.26082 0.38704<br><br></p>
\ No newline at end of file
+<p class='character'>Call: <code> LPCM(X = lpcmdat, mpoints = 2, groupvec = G) </code></p><p class='character'><h4>Coefficients:</h4></p>
+<p align= center >
+<table cellspacing=0 border=1><caption align=bottom class=captiondataframe></caption><tr><td>
+ <table border=0 class=dataframe>
+ <tbody> <tr class= firstline > <th></th><th>Eta</th><th>StdErr</th> </tr>
+ <tr><td class=firstcolumn>eta 1</td><td class=cellinside> 0.412</td><td class=cellinside> 0.397</td></tr>
+ <tr><td class=firstcolumn>eta 2</td><td class=cellinside>-0.299</td><td class=cellinside> 0.494</td></tr>
+ <tr><td class=firstcolumn>eta 3</td><td class=cellinside>-0.135</td><td class=cellinside> 0.409</td></tr>
+ <tr><td class=firstcolumn>eta 4</td><td class=cellinside> 0.035</td><td class=cellinside> 0.394</td></tr>
+ <tr><td class=firstcolumn>eta 5</td><td class=cellinside>-0.429</td><td class=cellinside> 0.478</td></tr>
+ <tr><td class=firstcolumn>eta 6</td><td class=cellinside>-0.286</td><td class=cellinside> 0.449</td></tr>
+ <tr><td class=firstcolumn>eta 7</td><td class=cellinside> 0.164</td><td class=cellinside> 0.399</td></tr>
+ <tr><td class=firstcolumn>eta 8</td><td class=cellinside> 0.101</td><td class=cellinside> 0.435</td></tr>
+ <tr><td class=firstcolumn>eta 9</td><td class=cellinside>-0.201</td><td class=cellinside> 0.261</td></tr>
+ <tr><td class=firstcolumn>eta 10</td><td class=cellinside> 1.094</td><td class=cellinside> 0.387</td></tr>
+
+ </tbody>
+</table>
+ </td></table>
+ <br>
+<p class='character'>Conditional log-likelihood: -103.7 <br />Number of iterations: 28 <br />Number of parameters: 10</p>
\ No newline at end of file
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.LRSM_parameter_estimation.rkcommands.R
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.LRSM_parameter_estimation.rkcommands.R 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.LRSM_parameter_estimation.rkcommands.R 2009-11-30 20:46:11 UTC (rev 2717)
@@ -5,7 +5,11 @@
estimates.lrsm <- LRSM(lrsmdat, mpoints=2, se=FALSE, sum0=FALSE)
## Print result
rk.header ("LRSM parameter estimation")
-rk.print (estimates.lrsm)
+rk.print (paste("Call: <code>",deparse(estimates.lrsm$call, width.cutoff=500),"</code>"))
+rk.print ("<h4>Coefficients:</h4>")
+rk.print(t(rbind(Eta=estimates.lrsm$etapar,StdErr=estimates.lrsm$se.eta)))
+rk.print (paste("Conditional log-likelihood:",round(estimates.lrsm$loglik, digits=1),
+"<br />Number of iterations:",estimates.lrsm$iter,"<br />Number of parameters:",estimates.lrsm$npar))
# keep results in current workspace
estimates.lrsm <<- estimates.lrsm
})
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.LRSM_parameter_estimation.rkout
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.LRSM_parameter_estimation.rkout 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.LRSM_parameter_estimation.rkout 2009-11-30 20:46:11 UTC (rev 2717)
@@ -1,3 +1,18 @@
<h1>LRSM parameter estimation</h1>
DATE<br>
-<p class='character'><br>Results of LRSM estimation: <br><br>Call: LRSM(X = lrsmdat, mpoints = 2, se = FALSE, sum0 = FALSE) <br><br>Conditional log-likelihood: -109.32 <br>Number of iterations: 11 <br>Number of parameters: 5 <br><br>Basic Parameters eta:<br> eta 1 eta 2 eta 3 eta 4 eta 5<br>Estimate -0.045897 0.11393 0.19745 -0.041557 0.34602<br>Std.Err NA NA NA NA NA<br><br></p>
\ No newline at end of file
+<p class='character'>Call: <code> LRSM(X = lrsmdat, mpoints = 2, se = FALSE, sum0 = FALSE) </code></p><p class='character'><h4>Coefficients:</h4></p>
+<p align= center >
+<table cellspacing=0 border=1><caption align=bottom class=captiondataframe></caption><tr><td>
+ <table border=0 class=dataframe>
+ <tbody> <tr class= firstline > <th></th><th>Eta</th><th>StdErr</th> </tr>
+ <tr><td class=firstcolumn>eta 1</td><td class=cellinside>-0.046</td><td class=cellinside> NA</td></tr>
+ <tr><td class=firstcolumn>eta 2</td><td class=cellinside> 0.114</td><td class=cellinside> NA</td></tr>
+ <tr><td class=firstcolumn>eta 3</td><td class=cellinside> 0.197</td><td class=cellinside> NA</td></tr>
+ <tr><td class=firstcolumn>eta 4</td><td class=cellinside>-0.042</td><td class=cellinside> NA</td></tr>
+ <tr><td class=firstcolumn>eta 5</td><td class=cellinside> 0.346</td><td class=cellinside> NA</td></tr>
+
+ </tbody>
+</table>
+ </td></table>
+ <br>
+<p class='character'>Conditional log-likelihood: -109.3 <br />Number of iterations: 11 <br />Number of parameters: 5</p>
\ No newline at end of file
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.PCM_parameter_estimation.rkcommands.R
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.PCM_parameter_estimation.rkcommands.R 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.PCM_parameter_estimation.rkcommands.R 2009-11-30 20:46:11 UTC (rev 2717)
@@ -5,7 +5,11 @@
estimates.pcm <- PCM(pcmdat)
## Print result
rk.header ("PCM parameter estimation")
-rk.print (estimates.pcm)
+rk.print (paste("Call: <code>",deparse(estimates.pcm$call, width.cutoff=500),"</code>"))
+rk.print ("<h4>Coefficients:</h4>")
+rk.print(t(rbind(Eta=estimates.pcm$etapar,StdErr=estimates.pcm$se.eta)))
+rk.print (paste("Conditional log-likelihood:",round(estimates.pcm$loglik, digits=1),
+"<br />Number of iterations:",estimates.pcm$iter,"<br />Number of parameters:",estimates.pcm$npar))
# keep results in current workspace
estimates.pcm <<- estimates.pcm
})
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.PCM_parameter_estimation.rkout
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.PCM_parameter_estimation.rkout 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.PCM_parameter_estimation.rkout 2009-11-30 20:46:11 UTC (rev 2717)
@@ -1,3 +1,32 @@
<h1>PCM parameter estimation</h1>
DATE<br>
-<p class='character'><br>Results of PCM estimation: <br><br>Call: PCM(X = pcmdat) <br><br>Conditional log-likelihood: -121.10 <br>Number of iterations: 16 <br>Number of parameters: 19 <br><br>Basic Parameters eta:<br> eta 1 eta 2 eta 3 eta 4 eta 5 eta 6 eta 7 eta 8<br>Estimate 0.93685 -0.33032 0.25086 0.065013 -1.4316 -0.22394 -0.074476 -0.87828<br>Std.Err 0.56521 0.80941 0.52546 0.579792 1.0239 0.56099 0.567946 0.78393<br> eta 9 eta 10 eta 11 eta 12 eta 13 eta 14 eta 15 eta 16<br>Estimate 0.42298 -0.26350 0.13957 -0.98036 0.096206 0.40303 -0.58611 0.50414<br>Std.Err 0.54673 0.68881 0.64970 0.78801 0.584619 0.57883 0.82641 0.63450<br> eta 17 eta 18 eta 19<br>Estimate 1.09952 0.84933 -0.29809<br>Std.Err 0.58696 0.50443 0.70054<br><br></p>
\ No newline at end of file
+<p class='character'>Call: <code> PCM(X = pcmdat) </code></p><p class='character'><h4>Coefficients:</h4></p>
+<p align= center >
+<table cellspacing=0 border=1><caption align=bottom class=captiondataframe></caption><tr><td>
+ <table border=0 class=dataframe>
+ <tbody> <tr class= firstline > <th></th><th>Eta</th><th>StdErr</th> </tr>
+ <tr><td class=firstcolumn>eta 1</td><td class=cellinside> 0.937</td><td class=cellinside> 0.565</td></tr>
+ <tr><td class=firstcolumn>eta 2</td><td class=cellinside>-0.330</td><td class=cellinside> 0.809</td></tr>
+ <tr><td class=firstcolumn>eta 3</td><td class=cellinside> 0.251</td><td class=cellinside> 0.525</td></tr>
+ <tr><td class=firstcolumn>eta 4</td><td class=cellinside> 0.065</td><td class=cellinside> 0.580</td></tr>
+ <tr><td class=firstcolumn>eta 5</td><td class=cellinside>-1.432</td><td class=cellinside> 1.024</td></tr>
+ <tr><td class=firstcolumn>eta 6</td><td class=cellinside>-0.224</td><td class=cellinside> 0.561</td></tr>
+ <tr><td class=firstcolumn>eta 7</td><td class=cellinside>-0.074</td><td class=cellinside> 0.568</td></tr>
+ <tr><td class=firstcolumn>eta 8</td><td class=cellinside>-0.878</td><td class=cellinside> 0.784</td></tr>
+ <tr><td class=firstcolumn>eta 9</td><td class=cellinside> 0.423</td><td class=cellinside> 0.547</td></tr>
+ <tr><td class=firstcolumn>eta 10</td><td class=cellinside>-0.263</td><td class=cellinside> 0.689</td></tr>
+ <tr><td class=firstcolumn>eta 11</td><td class=cellinside> 0.140</td><td class=cellinside> 0.650</td></tr>
+ <tr><td class=firstcolumn>eta 12</td><td class=cellinside>-0.980</td><td class=cellinside> 0.788</td></tr>
+ <tr><td class=firstcolumn>eta 13</td><td class=cellinside> 0.096</td><td class=cellinside> 0.585</td></tr>
+ <tr><td class=firstcolumn>eta 14</td><td class=cellinside> 0.403</td><td class=cellinside> 0.579</td></tr>
+ <tr><td class=firstcolumn>eta 15</td><td class=cellinside>-0.586</td><td class=cellinside> 0.826</td></tr>
+ <tr><td class=firstcolumn>eta 16</td><td class=cellinside> 0.504</td><td class=cellinside> 0.634</td></tr>
+ <tr><td class=firstcolumn>eta 17</td><td class=cellinside> 1.100</td><td class=cellinside> 0.587</td></tr>
+ <tr><td class=firstcolumn>eta 18</td><td class=cellinside> 0.849</td><td class=cellinside> 0.504</td></tr>
+ <tr><td class=firstcolumn>eta 19</td><td class=cellinside>-0.298</td><td class=cellinside> 0.701</td></tr>
+
+ </tbody>
+</table>
+ </td></table>
+ <br>
+<p class='character'>Conditional log-likelihood: -121.1 <br />Number of iterations: 16 <br />Number of parameters: 19</p>
\ No newline at end of file
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.RSM_parameter_estimation.rkcommands.R
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.RSM_parameter_estimation.rkcommands.R 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.RSM_parameter_estimation.rkcommands.R 2009-11-30 20:46:11 UTC (rev 2717)
@@ -5,7 +5,11 @@
estimates.rsm <- RSM(rsmdat)
## Print result
rk.header ("RSM parameter estimation")
-rk.print (estimates.rsm)
+rk.print (paste("Call: <code>",deparse(estimates.rsm$call, width.cutoff=500),"</code>"))
+rk.print ("<h4>Coefficients:</h4>")
+rk.print(t(rbind(Eta=estimates.rsm$etapar,StdErr=estimates.rsm$se.eta)))
+rk.print (paste("Conditional log-likelihood:",round(estimates.rsm$loglik, digits=1),
+"<br />Number of iterations:",estimates.rsm$iter,"<br />Number of parameters:",estimates.rsm$npar))
# keep results in current workspace
estimates.rsm <<- estimates.rsm
})
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.RSM_parameter_estimation.rkout
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.RSM_parameter_estimation.rkout 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.RSM_parameter_estimation.rkout 2009-11-30 20:46:11 UTC (rev 2717)
@@ -1,3 +1,20 @@
<h1>RSM parameter estimation</h1>
DATE<br>
-<p class='character'><br>Results of RSM estimation: <br><br>Call: RSM(X = rsmdat) <br><br>Conditional log-likelihood: -107.56 <br>Number of iterations: 11 <br>Number of parameters: 7 <br><br>Basic Parameters eta:<br> eta 1 eta 2 eta 3 eta 4 eta 5 eta 6 eta 7<br>Estimate -0.30515 -0.25586 -0.067308 0.16013 0.44424 0.16739 -0.13791<br>Std.Err 0.20537 0.20282 0.196507 0.19652 0.20849 0.45966 0.73249<br><br></p>
\ No newline at end of file
+<p class='character'>Call: <code> RSM(X = rsmdat) </code></p><p class='character'><h4>Coefficients:</h4></p>
+<p align= center >
+<table cellspacing=0 border=1><caption align=bottom class=captiondataframe></caption><tr><td>
+ <table border=0 class=dataframe>
+ <tbody> <tr class= firstline > <th></th><th>Eta</th><th>StdErr</th> </tr>
+ <tr><td class=firstcolumn>eta 1</td><td class=cellinside>-0.305</td><td class=cellinside> 0.205</td></tr>
+ <tr><td class=firstcolumn>eta 2</td><td class=cellinside>-0.256</td><td class=cellinside> 0.203</td></tr>
+ <tr><td class=firstcolumn>eta 3</td><td class=cellinside>-0.067</td><td class=cellinside> 0.197</td></tr>
+ <tr><td class=firstcolumn>eta 4</td><td class=cellinside> 0.160</td><td class=cellinside> 0.197</td></tr>
+ <tr><td class=firstcolumn>eta 5</td><td class=cellinside> 0.444</td><td class=cellinside> 0.208</td></tr>
+ <tr><td class=firstcolumn>eta 6</td><td class=cellinside> 0.167</td><td class=cellinside> 0.460</td></tr>
+ <tr><td class=firstcolumn>eta 7</td><td class=cellinside>-0.138</td><td class=cellinside> 0.732</td></tr>
+
+ </tbody>
+</table>
+ </td></table>
+ <br>
+<p class='character'>Conditional log-likelihood: -107.6 <br />Number of iterations: 11 <br />Number of parameters: 7</p>
\ No newline at end of file
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.Rasch_parameter_estimation.rkcommands.R
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.Rasch_parameter_estimation.rkcommands.R 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.Rasch_parameter_estimation.rkcommands.R 2009-11-30 20:46:11 UTC (rev 2717)
@@ -5,7 +5,10 @@
estimates.rasch <- rasch(LSAT)
## Print result
rk.header ("Rasch parameter estimation")
-rk.print (estimates.rasch$coefficients)
+rk.print (paste("Call: <code>",deparse(estimates.rasch$call, width.cutoff=500),"</code>"))
+rk.print ("<h4>Coefficients:</h4>")
+rk.print (coef(estimates.rasch))
+rk.print (paste("Log-likelihood value at convergence:",round(estimates.rasch$log.Lik, digits=1)))
# keep results in current workspace
estimates.rasch <<- estimates.rasch
})
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.Rasch_parameter_estimation.rkout
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.Rasch_parameter_estimation.rkout 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.Rasch_parameter_estimation.rkout 2009-11-30 20:46:11 UTC (rev 2717)
@@ -1,17 +1,18 @@
<h1>Rasch parameter estimation</h1>
DATE<br>
-
+<p class='character'>Call: <code> rasch(data = LSAT) </code></p><p class='character'><h4>Coefficients:</h4></p>
<p align= center >
<table cellspacing=0 border=1><caption align=bottom class=captiondataframe></caption><tr><td>
<table border=0 class=dataframe>
- <tbody> <tr class= firstline > <th></th><th>beta.i</th><th>beta</th> </tr>
- <tr><td class=firstcolumn>Item 1</td><td class=cellinside>2.73</td><td class=cellinside>0.76</td></tr>
- <tr><td class=firstcolumn>Item 2</td><td class=cellinside>1.00</td><td class=cellinside>0.76</td></tr>
- <tr><td class=firstcolumn>Item 3</td><td class=cellinside>0.24</td><td class=cellinside>0.76</td></tr>
- <tr><td class=firstcolumn>Item 4</td><td class=cellinside>1.31</td><td class=cellinside>0.76</td></tr>
- <tr><td class=firstcolumn>Item 5</td><td class=cellinside>2.10</td><td class=cellinside>0.76</td></tr>
+ <tbody> <tr class= firstline > <th></th><th>Dffclt</th><th>Dscrmn</th> </tr>
+ <tr><td class=firstcolumn>Item 1</td><td class=cellinside>-3.62</td><td class=cellinside> 0.76</td></tr>
+ <tr><td class=firstcolumn>Item 2</td><td class=cellinside>-1.32</td><td class=cellinside> 0.76</td></tr>
+ <tr><td class=firstcolumn>Item 3</td><td class=cellinside>-0.32</td><td class=cellinside> 0.76</td></tr>
+ <tr><td class=firstcolumn>Item 4</td><td class=cellinside>-1.73</td><td class=cellinside> 0.76</td></tr>
+ <tr><td class=firstcolumn>Item 5</td><td class=cellinside>-2.78</td><td class=cellinside> 0.76</td></tr>
</tbody>
</table>
</td></table>
<br>
+<p class='character'>Log-likelihood value at convergence: -2466.9</p>
\ No newline at end of file
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.goodnes-of-fit_Rasch.rkout
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.goodnes-of-fit_Rasch.rkout 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.goodnes-of-fit_Rasch.rkout 2009-11-30 20:46:11 UTC (rev 2717)
@@ -1,3 +1,3 @@
<h1>Goodness of Fit for Rasch Models (estimates.rasch)</h1>
DATE<br>
-<p class='character'><br>Bootstrap Goodness-of-Fit using Pearson chi-squared<br><br>Call:<br>rasch(data = LSAT)<br><br>Tobs: 18.34 <br># data-sets: 50 <br>p-value: 0.86 <br></p>
\ No newline at end of file
+<p class='character'><br>Bootstrap Goodness-of-Fit using Pearson chi-squared<br><br>Call:<br>rasch(data = LSAT)<br><br>Tobs: 18.34 <br># data-sets: 50 <br>p-value: 0.94 <br></p>
\ No newline at end of file
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.unidimensional.rkcommands.R
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.unidimensional.rkcommands.R 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.unidimensional.rkcommands.R 2009-11-30 20:46:11 UTC (rev 2717)
@@ -2,10 +2,12 @@
## Prepare
require(ltm)
## Compute
-unidim.res <<- unidimTest(estimates.rasch)
+unidim.res <- unidimTest(estimates.rasch)
## Print result
rk.header ("Unidimensionality check (estimates.rasch)")
rk.print (unidim.res)
+# keep results in current workspace
+unidim.res <<- unidim.res
})
-.rk.rerun.plugin.link(plugin="rkward::ltm_unidimensional", settings="spin_samples.real=100.00\nx.available=estimates.rasch", label="Run again")
+.rk.rerun.plugin.link(plugin="rkward::ltm_unidimensional", settings="chk_save.state=save\nsave_name.selection=unidim.res\nspin_samples.real=100.00\nx.available=estimates.rasch", label="Run again")
.rk.make.hr()
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.unidimensional.rkout
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.unidimensional.rkout 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.unidimensional.rkout 2009-11-30 20:46:11 UTC (rev 2717)
@@ -1,3 +1,3 @@
<h1>Unidimensionality check (estimates.rasch)</h1>
DATE<br>
-<p class='character'><br>Unidimensionality Check using Modified Parallel Analysis<br><br>Call:<br>rasch(data = LSAT)<br><br>Matrix of tertachoric correlations<br> Item 1 Item 2 Item 3 Item 4 Item 5<br>Item 1 1.000 0.170 0.228 0.107 0.067<br>Item 2 0.170 1.000 0.189 0.111 0.172<br>Item 3 0.228 0.189 1.000 0.187 0.105<br>Item 4 0.107 0.111 0.187 1.000 0.201<br>Item 5 0.067 0.172 0.105 0.201 1.000<br><br>Alternative hypothesis: the second eigenvalue of the observed data is substantially larger <br> than the second eigenvalue of data under the assumed IRT model<br><br>Second eigenvalue in the observed data: 0.225<br>Average of second eigenvalues in Monte Carlo samples: 0.249<br>Monte Carlo samples: 100<br>p-value: 0.614<br></p>
\ No newline at end of file
+<p class='character'><br>Unidimensionality Check using Modified Parallel Analysis<br><br>Call:<br>rasch(data = LSAT)<br><br>Matrix of tertachoric correlations<br> Item 1 Item 2 Item 3 Item 4 Item 5<br>Item 1 1.000 0.170 0.228 0.107 0.067<br>Item 2 0.170 1.000 0.189 0.111 0.172<br>Item 3 0.228 0.189 1.000 0.187 0.105<br>Item 4 0.107 0.111 0.187 1.000 0.201<br>Item 5 0.067 0.172 0.105 0.201 1.000<br><br>Alternative hypothesis: the second eigenvalue of the observed data is substantially larger <br> than the second eigenvalue of data under the assumed IRT model<br><br>Second eigenvalue in the observed data: 0.225<br>Average of second eigenvalues in Monte Carlo samples: 0.254<br>Monte Carlo samples: 100<br>p-value: 0.673<br></p>
\ No newline at end of file
Modified: trunk/rkward/tests/item_response_theory.R
===================================================================
--- trunk/rkward/tests/item_response_theory.R 2009-11-30 18:27:43 UTC (rev 2716)
+++ trunk/rkward/tests/item_response_theory.R 2009-11-30 20:46:11 UTC (rev 2717)
@@ -22,6 +22,7 @@
# polytomous data:
data("rsmdat") # rating scale model
data("pcmdat") # partial credit model
+ data("Science") # generalized partial credit model
data("lrsmdat") # linear rating scale model
data("lpcmdat") # linear partial credit model
},
@@ -61,6 +62,10 @@
rk.call.plugin ("rkward::par_est_pcm", chk_save.state="save", design.string="auto", etastart.string="NULL", stderr.state="se", sumnull.state="sum0", x.available="pcmdat", submit.mode="submit")
estimates$pcm <<- estimates.pcm
}),
+ new ("RKTest", id="GPCM_parameter_estimation", call=function () {
+ rk.call.plugin ("rkward::par_est_gpcm", chk_save.state="save", chk_select.state="select", constraint.string="rasch", epshess.real="0.000001", ghk_gpcm.real="21.00", inp_items.available="Science[[\"Work\"]]\nScience[[\"Industry\"]]\nScience[[\"Future\"]]\nScience[[\"Benefit\"]]", irtparam.state="TRUE", iterqn_gpcm.real="150.00", naaction.state="", numrderiv.string="fd", optimeth.string="BFGS", optimizer.string="optim", save_name.selection="estimates.gpcm", startval.string="NULL", verbose.state="", x.available="Science", submit.mode="submit")
+ estimates$gpcm <<- estimates.gpcm
+ }),
new ("RKTest", id="LRSM_parameter_estimation", call=function () {
rk.call.plugin ("rkward::par_est_lrsm", chk_save.state="save", design.string="auto", etastart.string="NULL", groups.string="1", mpoints.real="2.00", stderr.state="", sumnull.state="", x.available="lrsmdat", submit.mode="submit")
}),
@@ -73,7 +78,7 @@
## testing cronbach's alpha
new ("RKTest", id="Cronbach_alpha", call=function () {
rk.call.plugin ("rkward::ltm_cronbach_alpha", chk_bsci.state="bsci", chk_na.state="", chk_select.state="select", chk_standard.state="", inp_items.available="LSAT[[\"Item 1\"]]\nLSAT[[\"Item 2\"]]\nLSAT[[\"Item 3\"]]\nLSAT[[\"Item 4\"]]\nLSAT[[\"Item 5\"]]", spin_ci.real="0.95", spin_samples.real="1000.00", x.available="LSAT", submit.mode="submit")
- }, fuzzy_output=TRUE), # calculating the bootstrap confidence interval isn't always that accurate
+ }, fuzzy_output=TRUE), # bootstrap, varies a little
## now that our estimates are calculated, let's have a look at goodnes of fit
new ("RKTest", id="goodnes-of-fit_Rasch", call=function () {
@@ -147,7 +152,7 @@
# like initCalls: run after all tests to clean up.
postCalls = list (
function(){
- rm(list=c("LSAT","WIRS","lltmdat1","Environment","pcmdat","rsmdat","lrsmdat","lpcmdat","estimates"), envir=globalenv())
+ rm(list=c("LSAT","WIRS","lltmdat1","Environment","Science","pcmdat","rsmdat","lrsmdat","lpcmdat","estimates"), envir=globalenv())
}
)
)
This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.
More information about the rkward-tracker
mailing list