[rkward-cvs] SF.net SVN: rkward:[2959] trunk/rkward
m-eik at users.sourceforge.net
m-eik at users.sourceforge.net
Fri Aug 6 15:06:46 UTC 2010
Revision: 2959
http://rkward.svn.sourceforge.net/rkward/?rev=2959&view=rev
Author: m-eik
Date: 2010-08-06 15:06:45 +0000 (Fri, 06 Aug 2010)
Log Message:
-----------
irt: updated several dialogs and test outputs, fixed Wald test inconsistency
Modified Paths:
--------------
trunk/rkward/rkward/plugins/analysis/irt/tests/eRm_waldtest.js
trunk/rkward/rkward/plugins/analysis/irt/tests/eRm_waldtest.xml
trunk/rkward/rkward/plugins/analysis/irt/tests/ltm_gof_rasch.js
trunk/rkward/rkward/plugins/analysis/irt/tests/ltm_item_fit.js
trunk/rkward/rkward/plugins/analysis/irt/tests/ltm_person_fit.js
trunk/rkward/rkward/plugins/analysis/irt/tests/ltm_unidimensional.js
trunk/rkward/tests/item_response_theory/RKTestStandard.2PL_parameter_estimation.rkout
trunk/rkward/tests/item_response_theory/RKTestStandard.3PL_parameter_estimation.rkout
trunk/rkward/tests/item_response_theory/RKTestStandard.Cronbach_alpha.rkout
trunk/rkward/tests/item_response_theory/RKTestStandard.GPCM_parameter_estimation.rkout
trunk/rkward/tests/item_response_theory/RKTestStandard.GRM_parameter_estimation.rkout
trunk/rkward/tests/item_response_theory/RKTestStandard.LLTM_parameter_estimation.rkout
trunk/rkward/tests/item_response_theory/RKTestStandard.LPCM_parameter_estimation.rkout
trunk/rkward/tests/item_response_theory/RKTestStandard.LRSM_parameter_estimation.rkout
trunk/rkward/tests/item_response_theory/RKTestStandard.PCM_parameter_estimation.rkout
trunk/rkward/tests/item_response_theory/RKTestStandard.RSM_parameter_estimation.rkout
trunk/rkward/tests/item_response_theory/RKTestStandard.Rasch_parameter_estimation.rkout
trunk/rkward/tests/item_response_theory/RKTestStandard.Wald_test.messages.txt
trunk/rkward/tests/item_response_theory/RKTestStandard.Wald_test.rkcommands.R
trunk/rkward/tests/item_response_theory/RKTestStandard.Wald_test.rkout
trunk/rkward/tests/item_response_theory/RKTestStandard.goodnes-of-fit_Rasch.rkcommands.R
trunk/rkward/tests/item_response_theory/RKTestStandard.goodnes-of-fit_Rasch.rkout
trunk/rkward/tests/item_response_theory/RKTestStandard.item_fit_statistics.rkcommands.R
trunk/rkward/tests/item_response_theory/RKTestStandard.item_fit_statistics.rkout
trunk/rkward/tests/item_response_theory/RKTestStandard.person_fit_statistics.rkcommands.R
trunk/rkward/tests/item_response_theory/RKTestStandard.person_fit_statistics.rkout
trunk/rkward/tests/item_response_theory/RKTestStandard.unidimensional.rkcommands.R
trunk/rkward/tests/item_response_theory/RKTestStandard.unidimensional.rkout
Modified: trunk/rkward/rkward/plugins/analysis/irt/tests/eRm_waldtest.js
===================================================================
--- trunk/rkward/rkward/plugins/analysis/irt/tests/eRm_waldtest.js 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/rkward/plugins/analysis/irt/tests/eRm_waldtest.js 2010-08-06 15:06:45 UTC (rev 2959)
@@ -7,15 +7,21 @@
// let's read all values into php variables for the sake of readable code
var rad_splitcr = getValue("rad_splitcr");
var splitvector = getValue("splitvector");
+ var optimizer = getValue("drop_optimizer");
+ if (optimizer != "nlm") echo('.GlobalEnv$fitctrl <- \"'+optimizer+'\"\n');
echo ('waldtest.res <- Waldtest(' + getValue("x"));
// check if any advanced control options must be inserted
- if (rad_splitcr == "mean") echo(", splitcr=\"mean\"") ;
- if (rad_splitcr == "vector" && splitvector) echo(", splitcr="+splitvector) ;
+ if (rad_splitcr == "mean") echo(", splitcr=\"mean\"");
+ if (rad_splitcr == "vector" && splitvector) echo(", splitcr="+splitvector);
echo (')\n');
+ if (optimizer != "nlm") echo('rm(fitctrl, envir=.GlobalEnv)\n');
}
function printout () {
echo ('rk.header ("Wald test (' + getValue("x") + ')")\n');
- echo ('rk.print (waldtest.res)\n');
+ echo ('rk.print ("Call:")\n');
+ echo ('rk.print.literal (deparse(waldtest.res$call, width.cutoff=500))\n');
+ echo ('rk.header ("Wald test on item level (z-values):", level=4)\n');
+ echo ('rk.print(waldtest.res$coef.table)\n');
}
Modified: trunk/rkward/rkward/plugins/analysis/irt/tests/eRm_waldtest.xml
===================================================================
--- trunk/rkward/rkward/plugins/analysis/irt/tests/eRm_waldtest.xml 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/rkward/plugins/analysis/irt/tests/eRm_waldtest.xml 2010-08-06 15:06:45 UTC (rev 2959)
@@ -7,25 +7,39 @@
<convert id="split_vec" mode="equals" sources="rad_splitcr.string" standard="vector" />
<connect client="splitvector.enabled" governor="split_vec" />
+ <connect client="vars2.enabled" governor="split_vec" />
<connect client="splitvector.required" governor="split_vec" />
</logic>
<dialog label="Wald test">
<tabbook>
- <tab id="tab_var_options" label="Data">
+ <tab id="tab_var_data" label="Data">
<row id="row_vars">
<varselector id="vars" />
<column>
<varslot type="numeric" id="x" classes="Rm" source="vars" required="true" label="Chose fitted model to test:"/>
<text>Only objects of class Rm are valid!</text>
<stretch />
+ </column>
+ </row>
+ </tab>
+ <tab id="tab_var_options" label="Options">
+ <row id="row_vars">
+ <varselector id="vars2" />
+ <column>
<radio id="rad_splitcr" label="Split criterion for subject raw scores">
<option label="Median" value="median" checked="true" />
<option label="Mean" value="mean" />
<option label="Use dichotomous vector to define groups" value="vector" />
</radio>
- <varslot type="numeric" id="splitvector" source="vars" label="Select grouping vector:"/>
- </column>
+ <varslot type="numeric" id="splitvector" source="vars2" label="Select grouping vector:"/>
+ <frame id="optimizer" label="Optimizer">
+ <dropdown id="drop_optimizer" label="Used optimizer:">
+ <option label="nlm()" value="nlm" checked="true" />
+ <option label="optim()" value="optim" />
+ </dropdown>
+ </frame>
+ </column>
</row>
</tab>
</tabbook>
Modified: trunk/rkward/rkward/plugins/analysis/irt/tests/ltm_gof_rasch.js
===================================================================
--- trunk/rkward/rkward/plugins/analysis/irt/tests/ltm_gof_rasch.js 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/rkward/plugins/analysis/irt/tests/ltm_gof_rasch.js 2010-08-06 15:06:45 UTC (rev 2959)
@@ -15,5 +15,10 @@
function printout () {
echo ('rk.header ("Goodness of Fit for Rasch Models (' + getValue("x") + ')")\n');
- echo ('rk.print (GoFRasch.res)\n');
+ echo ('rk.print ("Call:")\n');
+ echo ('rk.print.literal (deparse(GoFRasch.res$call, width.cutoff=500))\n');
+ echo ('rk.header ("Parametric Bootstrap test:", level=4)\n');
+ echo ('rk.print (paste("Chi-squared statistic:", round(GoFRasch.res$Tobs, digits=3)))\n');
+ echo ('rk.print (paste("Bootstrap samples:", GoFRasch.res$B))\n');
+ echo ('rk.print (paste("p-value:", GoFRasch.res$p.value))\n');
}
Modified: trunk/rkward/rkward/plugins/analysis/irt/tests/ltm_item_fit.js
===================================================================
--- trunk/rkward/rkward/plugins/analysis/irt/tests/ltm_item_fit.js 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/rkward/plugins/analysis/irt/tests/ltm_item_fit.js 2010-08-06 15:06:45 UTC (rev 2959)
@@ -30,6 +30,15 @@
}
function printout () {
+ var rad_pvalue = getValue("rad_pvalue");
+
echo ('rk.header ("Item-fit statistics (' + getValue("x") + ')")\n');
- echo ('rk.print (itemfit.res)\n');
+ echo ('rk.print ("Call:")\n');
+ echo ('rk.print.literal (deparse(itemfit.res$call, width.cutoff=500))\n');
+ echo ('rk.print ("Alternative: Items do not fit the model")\n');
+ echo ('rk.print (paste("Ability Categories:", itemfit.res$G))\n');
+ if (rad_pvalue == "montecarlo")
+ echo ('rk.print (paste("Monte Carlo samples:", itemfit.res$B))\n');
+ echo ('rk.header ("Item-Fit Statistics and P-values:", level=4)\n');
+ echo ('rk.print(cbind("X^2"=round(itemfit.res$Tobs, digits=3), "Pr (>X^2)"=format(round(itemfit.res$p.values, digits=3), nsmall=3)))\n');
}
Modified: trunk/rkward/rkward/plugins/analysis/irt/tests/ltm_person_fit.js
===================================================================
--- trunk/rkward/rkward/plugins/analysis/irt/tests/ltm_person_fit.js 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/rkward/plugins/analysis/irt/tests/ltm_person_fit.js 2010-08-06 15:06:45 UTC (rev 2959)
@@ -32,5 +32,8 @@
function printout () {
echo ('rk.header ("Person-fit statistics (' + getValue("x") + ')")\n');
- echo ('rk.print (personfit.res)\n');
+ echo ('rk.print ("Call:")\n');
+ echo ('rk.print.literal (deparse(personfit.res$call, width.cutoff=500))\n');
+ echo ('rk.header ("Response patterns, person-fit statistics (L0, Lz) and probabilities for each response pattern (Pr):", level=4)\n');
+ echo ('rk.print(cbind(format(personfit.res$resp.patterns, nsmall=0), round(personfit.res$Tobs, digits=3), "Pr (<Lz)"=round(c(personfit.res$p.values), digits=3)))\n');
}
Modified: trunk/rkward/rkward/plugins/analysis/irt/tests/ltm_unidimensional.js
===================================================================
--- trunk/rkward/rkward/plugins/analysis/irt/tests/ltm_unidimensional.js 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/rkward/plugins/analysis/irt/tests/ltm_unidimensional.js 2010-08-06 15:06:45 UTC (rev 2959)
@@ -18,7 +18,16 @@
var save_name = getValue("save_name");
echo ('rk.header ("Unidimensionality check (' + getValue("x") + ')")\n');
- echo ('rk.print (unidim.res)\n');
+ echo ('rk.print ("Call:")\n');
+ echo ('rk.print.literal (deparse(unidim.res$call, width.cutoff=500))\n');
+ echo ('rk.header ("Matrix of tertachoric correlations:", level=4)\n');
+ echo ('rk.print (unidim.res$Rho)\n');
+ echo ('rk.header ("Unidimensionality Check using Modified Parallel Analysis:", level=4)\n');
+ echo ('rk.print ("Alternative hypothesis: <em>The second eigenvalue of the observed data is substantially larger than the second eigenvalue of data under the assumed IRT model</em>")\n');
+ echo ('rk.print (paste("Second eigenvalue in the observed data:", round(unidim.res$Tobs[2], digits=3)))\n');
+ echo ('rk.print (paste("Average of second eigenvalues in Monte Carlo samples:", round(mean(unidim.res$T.boot[,2]), digits=3)))\n');
+ echo ('rk.print (paste("Monte Carlo samples:", dim(unidim.res$T.boot)[1]))\n');
+ echo ('rk.print (paste("p-value:", round(unidim.res$p.value, digits=3)))\n');
// check if results are to be saved:
if (save && save_name) {
echo ('# keep results in current workspace\n');
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.2PL_parameter_estimation.rkout
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.2PL_parameter_estimation.rkout 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.2PL_parameter_estimation.rkout 2010-08-06 15:06:45 UTC (rev 2959)
@@ -1,22 +1,57 @@
<h1>2PL parameter estimation</h1>
DATE<br>
-<p class='character'>Call:</p><pre>ltm(formula = WIRS ~ z1 * z2)</pre>
+
+<p class='character'>Call:</p>
+<pre>ltm(formula = WIRS ~ z1 * z2)</pre>
<h4>Coefficients:</h4>
<br>
+
<p align= center >
-<table cellspacing=0 border=1><caption align=bottom class=captiondataframe></caption><tr><td>
+<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>
+<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
+
+<p class='character'>Log-likelihood value at convergence: -3293.2</p>
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.3PL_parameter_estimation.rkout
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.3PL_parameter_estimation.rkout 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.3PL_parameter_estimation.rkout 2010-08-06 15:06:45 UTC (rev 2959)
@@ -1,21 +1,46 @@
<h1>3PL parameter estimation</h1>
DATE<br>
-<p class='character'>Call:</p><pre>tpm(data = LSAT)</pre>
+
+<p class='character'>Call:</p>
+<pre>tpm(data = LSAT)</pre>
<h4>Coefficients:</h4>
<br>
+
<p align= center >
-<table cellspacing=0 border=1><caption align=bottom class=captiondataframe></caption><tr><td>
+<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>
+<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
+
+<p class='character'>Log-likelihood value at convergence: -2466.7</p>
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.Cronbach_alpha.rkout
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.Cronbach_alpha.rkout 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.Cronbach_alpha.rkout 2010-08-06 15:06:45 UTC (rev 2959)
@@ -1,29 +1,56 @@
<h1>Cronbach's alpha</h1>
DATE<br>
-<p class='character'>for the 'LSAT' data-set (subset: Item 1, Item 2, Item 3, Item 4, Item 5)</p><p class='character'>Items: 5 <br />Sample units: 1000 <br /><strong>alpha: 0.29 </strong></p><p class='character'>Effects on alpha if items are removed:</p>
+
+<p class='character'>for the 'LSAT' data-set (subset: Item 1, Item 2, Item 3, Item 4, Item 5)</p>
+
+<p class='character'>Items: 5 <br />Sample units: 1000 <br /><strong>alpha: 0.29 </strong></p>
+
+<p class='character'>Effects on alpha if items are removed:</p>
+
+
<p align= center >
-<table cellspacing=0 border=1><caption align=bottom class=captiondataframe></caption><tr><td>
+<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>value</th> </tr>
- <tr><td class=firstcolumn>All Items</td><td class=cellinside>0.29</td></tr>
- <tr><td class=firstcolumn>Excluding Item 1</td><td class=cellinside>0.28</td></tr>
- <tr><td class=firstcolumn>Excluding Item 2</td><td class=cellinside>0.24</td></tr>
- <tr><td class=firstcolumn>Excluding Item 3</td><td class=cellinside>0.22</td></tr>
- <tr><td class=firstcolumn>Excluding Item 4</td><td class=cellinside>0.25</td></tr>
- <tr><td class=firstcolumn>Excluding Item 5</td><td class=cellinside>0.27</td></tr>
+<tr><td class=firstcolumn>All Items</td>
+<td class=cellinside>0.29</td></tr>
+
+<tr><td class=firstcolumn>Excluding Item 1</td>
+<td class=cellinside>0.28</td></tr>
+
+<tr><td class=firstcolumn>Excluding Item 2</td>
+<td class=cellinside>0.24</td></tr>
+
+<tr><td class=firstcolumn>Excluding Item 3</td>
+<td class=cellinside>0.22</td></tr>
+
+<tr><td class=firstcolumn>Excluding Item 4</td>
+<td class=cellinside>0.25</td></tr>
+
+<tr><td class=firstcolumn>Excluding Item 5</td>
+<td class=cellinside>0.27</td></tr>
+
</tbody>
</table>
</td></table>
<br>
+
<p class='character'>95% Confidence interval:</p>
+
+
<p align= center >
-<table cellspacing=0 border=1><caption align=bottom class=captiondataframe></caption><tr><td>
+<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>2.5%</th><th>97.5%</th> </tr>
- <tr><td class=firstcolumn></td><td class=cellinside>0.22</td><td class=cellinside>0.36</td></tr>
+<tr><td class=firstcolumn></td>
+<td class=cellinside>0.22</td>
+<td class=cellinside>0.36</td></tr>
+
</tbody>
</table>
</td></table>
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.GPCM_parameter_estimation.rkout
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.GPCM_parameter_estimation.rkout 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.GPCM_parameter_estimation.rkout 2010-08-06 15:06:45 UTC (rev 2959)
@@ -1,20 +1,45 @@
<h1>GPCM parameter estimation</h1>
DATE<br>
-<p class='character'>Call:</p><pre>gpcm(data = subset(Science, select = c("Work", "Industry", "Future", "Benefit")), constraint = "rasch")</pre>
+
+<p class='character'>Call:</p>
+<pre>gpcm(data = subset(Science, select = c("Work", "Industry", "Future", "Benefit")), constraint = "rasch")</pre>
<h4>Coefficients:</h4>
<br>
+
<p align= center >
-<table cellspacing=0 border=1><caption align=bottom class=captiondataframe></caption><tr><td>
+<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>Catgr.1</th><th>Catgr.2</th><th>Catgr.3</th><th>Dscrmn</th> </tr>
- <tr><td class=firstcolumn>Work</td><td class=cellinside>-1.80</td><td class=cellinside>-0.90</td><td class=cellinside> 1.80</td><td class=cellinside> 1.00</td></tr>
- <tr><td class=firstcolumn>Industry</td><td class=cellinside>-2.54</td><td class=cellinside>-1.76</td><td class=cellinside> 0.20</td><td class=cellinside> 1.00</td></tr>
- <tr><td class=firstcolumn>Future</td><td class=cellinside>-2.50</td><td class=cellinside>-1.39</td><td class=cellinside> 1.09</td><td class=cellinside> 1.00</td></tr>
- <tr><td class=firstcolumn>Benefit</td><td class=cellinside>-2.33</td><td class=cellinside>-0.88</td><td class=cellinside> 1.31</td><td class=cellinside> 1.00</td></tr>
+<tr><td class=firstcolumn>Work</td>
+<td class=cellinside>-1.80</td>
+<td class=cellinside>-0.90</td>
+<td class=cellinside> 1.80</td>
+<td class=cellinside> 1.00</td></tr>
+
+<tr><td class=firstcolumn>Industry</td>
+<td class=cellinside>-2.54</td>
+<td class=cellinside>-1.76</td>
+<td class=cellinside> 0.20</td>
+<td class=cellinside> 1.00</td></tr>
+
+<tr><td class=firstcolumn>Future</td>
+<td class=cellinside>-2.50</td>
+<td class=cellinside>-1.39</td>
+<td class=cellinside> 1.09</td>
+<td class=cellinside> 1.00</td></tr>
+
+<tr><td class=firstcolumn>Benefit</td>
+<td class=cellinside>-2.33</td>
+<td class=cellinside>-0.88</td>
+<td class=cellinside> 1.31</td>
+<td class=cellinside> 1.00</td></tr>
+
</tbody>
</table>
</td></table>
<br>
-<p class='character'>Log-likelihood value at convergence: -1756.8</p>
\ No newline at end of file
+
+<p class='character'>Log-likelihood value at convergence: -1756.8</p>
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.GRM_parameter_estimation.rkout
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.GRM_parameter_estimation.rkout 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.GRM_parameter_estimation.rkout 2010-08-06 15:06:45 UTC (rev 2959)
@@ -1,22 +1,51 @@
<h1>GRM parameter estimation</h1>
DATE<br>
-<p class='character'>Call:</p><pre>grm(data = Environment)</pre>
+
+<p class='character'>Call:</p>
+<pre>grm(data = Environment)</pre>
<h4>Coefficients:</h4>
<br>
+
<p align= center >
-<table cellspacing=0 border=1><caption align=bottom class=captiondataframe></caption><tr><td>
+<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>
+<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
+
+<p class='character'>Log-likelihood value at convergence: -1090.4</p>
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.LLTM_parameter_estimation.rkout
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.LLTM_parameter_estimation.rkout 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.LLTM_parameter_estimation.rkout 2010-08-06 15:06:45 UTC (rev 2959)
@@ -1,31 +1,81 @@
<h1>LLTM parameter estimation</h1>
DATE<br>
-<p class='character'>Call:</p><pre>LLTM(X = lltmdat1, mpoints = 2)</pre>
+
+<p class='character'>Call:</p>
+<pre>LLTM(X = lltmdat1, mpoints = 2)</pre>
<h4>Coefficients:</h4>
<br>
+
<p align= center >
-<table cellspacing=0 border=1><caption align=bottom class=captiondataframe></caption><tr><td>
+<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>
+<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
+
+<p class='character'>Conditional log-likelihood: -1511.2 <br />Number of iterations: 15 <br />Number of parameters: 15</p>
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.LPCM_parameter_estimation.rkout
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.LPCM_parameter_estimation.rkout 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.LPCM_parameter_estimation.rkout 2010-08-06 15:06:45 UTC (rev 2959)
@@ -1,26 +1,61 @@
<h1>LPCM parameter estimation</h1>
DATE<br>
-<p class='character'>Call:</p><pre>LPCM(X = lpcmdat, mpoints = 2, groupvec = G)</pre>
+
+<p class='character'>Call:</p>
+<pre>LPCM(X = lpcmdat, mpoints = 2, groupvec = G)</pre>
<h4>Coefficients:</h4>
<br>
+
<p align= center >
-<table cellspacing=0 border=1><caption align=bottom class=captiondataframe></caption><tr><td>
+<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>
+<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
+
+<p class='character'>Conditional log-likelihood: -103.7 <br />Number of iterations: 28 <br />Number of parameters: 10</p>
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.LRSM_parameter_estimation.rkout
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.LRSM_parameter_estimation.rkout 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.LRSM_parameter_estimation.rkout 2010-08-06 15:06:45 UTC (rev 2959)
@@ -1,21 +1,41 @@
<h1>LRSM parameter estimation</h1>
DATE<br>
-<p class='character'>Call:</p><pre>LRSM(X = lrsmdat, mpoints = 2, se = FALSE, sum0 = FALSE)</pre>
+
+<p class='character'>Call:</p>
+<pre>LRSM(X = lrsmdat, mpoints = 2, se = FALSE, sum0 = FALSE)</pre>
<h4>Coefficients:</h4>
<br>
+
<p align= center >
-<table cellspacing=0 border=1><caption align=bottom class=captiondataframe></caption><tr><td>
+<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>
+<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
+
+<p class='character'>Conditional log-likelihood: -109.3 <br />Number of iterations: 11 <br />Number of parameters: 5</p>
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.PCM_parameter_estimation.rkout
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.PCM_parameter_estimation.rkout 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.PCM_parameter_estimation.rkout 2010-08-06 15:06:45 UTC (rev 2959)
@@ -1,35 +1,97 @@
<h1>PCM parameter estimation</h1>
DATE<br>
-<p class='character'>Call:</p><pre>PCM(X = pcmdat)</pre>
+
+<p class='character'>Call:</p>
+<pre>PCM(X = pcmdat)</pre>
<h4>Coefficients:</h4>
<br>
+
<p align= center >
-<table cellspacing=0 border=1><caption align=bottom class=captiondataframe></caption><tr><td>
+<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>I1.c2</td><td class=cellinside>-0.937</td><td class=cellinside> 0.565</td></tr>
- <tr><td class=firstcolumn>I1.c3</td><td class=cellinside> 0.330</td><td class=cellinside> 0.809</td></tr>
- <tr><td class=firstcolumn>I2.c1</td><td class=cellinside>-0.251</td><td class=cellinside> 0.525</td></tr>
- <tr><td class=firstcolumn>I2.c2</td><td class=cellinside>-0.065</td><td class=cellinside> 0.580</td></tr>
- <tr><td class=firstcolumn>I2.c3</td><td class=cellinside> 1.432</td><td class=cellinside> 1.024</td></tr>
- <tr><td class=firstcolumn>I3.c1</td><td class=cellinside> 0.224</td><td class=cellinside> 0.561</td></tr>
- <tr><td class=firstcolumn>I3.c2</td><td class=cellinside> 0.074</td><td class=cellinside> 0.568</td></tr>
- <tr><td class=firstcolumn>I3.c3</td><td class=cellinside> 0.878</td><td class=cellinside> 0.784</td></tr>
- <tr><td class=firstcolumn>I4.c1</td><td class=cellinside>-0.423</td><td class=cellinside> 0.547</td></tr>
- <tr><td class=firstcolumn>I4.c2</td><td class=cellinside> 0.263</td><td class=cellinside> 0.689</td></tr>
- <tr><td class=firstcolumn>I4.c3</td><td class=cellinside>-0.140</td><td class=cellinside> 0.650</td></tr>
- <tr><td class=firstcolumn>I5.c1</td><td class=cellinside> 0.980</td><td class=cellinside> 0.788</td></tr>
- <tr><td class=firstcolumn>I5.c2</td><td class=cellinside>-0.096</td><td class=cellinside> 0.585</td></tr>
- <tr><td class=firstcolumn>I5.c3</td><td class=cellinside>-0.403</td><td class=cellinside> 0.579</td></tr>
- <tr><td class=firstcolumn>I6.c1</td><td class=cellinside> 0.586</td><td class=cellinside> 0.826</td></tr>
- <tr><td class=firstcolumn>I6.c2</td><td class=cellinside>-0.504</td><td class=cellinside> 0.634</td></tr>
- <tr><td class=firstcolumn>I6.c3</td><td class=cellinside>-1.100</td><td class=cellinside> 0.587</td></tr>
- <tr><td class=firstcolumn>I7.c1</td><td class=cellinside>-0.849</td><td class=cellinside> 0.504</td></tr>
- <tr><td class=firstcolumn>I7.c2</td><td class=cellinside> 0.298</td><td class=cellinside> 0.701</td></tr>
+<tr><td class=firstcolumn>I1.c2</td>
+<td class=cellinside>-0.937</td>
+<td class=cellinside> 0.565</td></tr>
+
+<tr><td class=firstcolumn>I1.c3</td>
+<td class=cellinside> 0.330</td>
+<td class=cellinside> 0.809</td></tr>
+
+<tr><td class=firstcolumn>I2.c1</td>
+<td class=cellinside>-0.251</td>
+<td class=cellinside> 0.525</td></tr>
+
+<tr><td class=firstcolumn>I2.c2</td>
+<td class=cellinside>-0.065</td>
+<td class=cellinside> 0.580</td></tr>
+
+<tr><td class=firstcolumn>I2.c3</td>
+<td class=cellinside> 1.432</td>
+<td class=cellinside> 1.024</td></tr>
+
+<tr><td class=firstcolumn>I3.c1</td>
+<td class=cellinside> 0.224</td>
+<td class=cellinside> 0.561</td></tr>
+
+<tr><td class=firstcolumn>I3.c2</td>
+<td class=cellinside> 0.074</td>
+<td class=cellinside> 0.568</td></tr>
+
+<tr><td class=firstcolumn>I3.c3</td>
+<td class=cellinside> 0.878</td>
+<td class=cellinside> 0.784</td></tr>
+
+<tr><td class=firstcolumn>I4.c1</td>
+<td class=cellinside>-0.423</td>
+<td class=cellinside> 0.547</td></tr>
+
+<tr><td class=firstcolumn>I4.c2</td>
+<td class=cellinside> 0.263</td>
+<td class=cellinside> 0.689</td></tr>
+
+<tr><td class=firstcolumn>I4.c3</td>
+<td class=cellinside>-0.140</td>
+<td class=cellinside> 0.650</td></tr>
+
+<tr><td class=firstcolumn>I5.c1</td>
+<td class=cellinside> 0.980</td>
+<td class=cellinside> 0.788</td></tr>
+
+<tr><td class=firstcolumn>I5.c2</td>
+<td class=cellinside>-0.096</td>
+<td class=cellinside> 0.585</td></tr>
+
+<tr><td class=firstcolumn>I5.c3</td>
+<td class=cellinside>-0.403</td>
+<td class=cellinside> 0.579</td></tr>
+
+<tr><td class=firstcolumn>I6.c1</td>
+<td class=cellinside> 0.586</td>
+<td class=cellinside> 0.826</td></tr>
+
+<tr><td class=firstcolumn>I6.c2</td>
+<td class=cellinside>-0.504</td>
+<td class=cellinside> 0.634</td></tr>
+
+<tr><td class=firstcolumn>I6.c3</td>
+<td class=cellinside>-1.100</td>
+<td class=cellinside> 0.587</td></tr>
+
+<tr><td class=firstcolumn>I7.c1</td>
+<td class=cellinside>-0.849</td>
+<td class=cellinside> 0.504</td></tr>
+
+<tr><td class=firstcolumn>I7.c2</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
+
+<p class='character'>Conditional log-likelihood: -121.1 <br />Number of iterations: 16 <br />Number of parameters: 19</p>
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.RSM_parameter_estimation.rkout
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.RSM_parameter_estimation.rkout 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.RSM_parameter_estimation.rkout 2010-08-06 15:06:45 UTC (rev 2959)
@@ -1,23 +1,49 @@
<h1>RSM parameter estimation</h1>
DATE<br>
-<p class='character'>Call:</p><pre>RSM(X = rsmdat)</pre>
+
+<p class='character'>Call:</p>
+<pre>RSM(X = rsmdat)</pre>
<h4>Coefficients:</h4>
<br>
+
<p align= center >
-<table cellspacing=0 border=1><caption align=bottom class=captiondataframe></caption><tr><td>
+<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>I2</td><td class=cellinside> 0.305</td><td class=cellinside> 0.205</td></tr>
- <tr><td class=firstcolumn>I3</td><td class=cellinside> 0.256</td><td class=cellinside> 0.203</td></tr>
- <tr><td class=firstcolumn>I4</td><td class=cellinside> 0.067</td><td class=cellinside> 0.197</td></tr>
- <tr><td class=firstcolumn>I5</td><td class=cellinside>-0.160</td><td class=cellinside> 0.197</td></tr>
- <tr><td class=firstcolumn>I6</td><td class=cellinside>-0.444</td><td class=cellinside> 0.208</td></tr>
- <tr><td class=firstcolumn>Cat 2</td><td class=cellinside>-0.167</td><td class=cellinside> 0.460</td></tr>
- <tr><td class=firstcolumn>Cat 3</td><td class=cellinside> 0.138</td><td class=cellinside> 0.732</td></tr>
+<tr><td class=firstcolumn>I2</td>
+<td class=cellinside> 0.305</td>
+<td class=cellinside> 0.205</td></tr>
+
+<tr><td class=firstcolumn>I3</td>
+<td class=cellinside> 0.256</td>
+<td class=cellinside> 0.203</td></tr>
+
+<tr><td class=firstcolumn>I4</td>
+<td class=cellinside> 0.067</td>
+<td class=cellinside> 0.197</td></tr>
+
+<tr><td class=firstcolumn>I5</td>
+<td class=cellinside>-0.160</td>
+<td class=cellinside> 0.197</td></tr>
+
+<tr><td class=firstcolumn>I6</td>
+<td class=cellinside>-0.444</td>
+<td class=cellinside> 0.208</td></tr>
+
+<tr><td class=firstcolumn>Cat 2</td>
+<td class=cellinside>-0.167</td>
+<td class=cellinside> 0.460</td></tr>
+
+<tr><td class=firstcolumn>Cat 3</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
+
+<p class='character'>Conditional log-likelihood: -107.6 <br />Number of iterations: 11 <br />Number of parameters: 7</p>
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.Rasch_parameter_estimation.rkout
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.Rasch_parameter_estimation.rkout 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.Rasch_parameter_estimation.rkout 2010-08-06 15:06:45 UTC (rev 2959)
@@ -1,21 +1,41 @@
<h1>Rasch parameter estimation</h1>
DATE<br>
-<p class='character'>Call:</p><pre>rasch(data = LSAT)</pre>
+
+<p class='character'>Call:</p>
+<pre>rasch(data = LSAT)</pre>
<h4>Coefficients:</h4>
<br>
+
<p align= center >
-<table cellspacing=0 border=1><caption align=bottom class=captiondataframe></caption><tr><td>
+<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>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>
+<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
+
+<p class='character'>Log-likelihood value at convergence: -2466.9</p>
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.Wald_test.messages.txt
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.Wald_test.messages.txt 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.Wald_test.messages.txt 2010-08-06 15:06:45 UTC (rev 2959)
@@ -1,7 +1,3 @@
Warning in Waldtest.Rm(estimates.rsm) :
The following items were excluded due to inappropriate response patterns within subgroups:
-Warning messages:
-1: In sqrt(diag(solve(parest$hessian))) : NaNs produced
-2: In sqrt(diag(lres$W %*% solve(parest$hessian) %*% t(lres$W))) :
- NaNs produced
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.Wald_test.rkcommands.R
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.Wald_test.rkcommands.R 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.Wald_test.rkcommands.R 2010-08-06 15:06:45 UTC (rev 2959)
@@ -2,10 +2,15 @@
## Prepare
require(eRm)
## Compute
+.GlobalEnv$fitctrl <- "optim"
waldtest.res <- Waldtest(estimates.rsm)
+rm(fitctrl, envir=.GlobalEnv)
## Print result
rk.header ("Wald test (estimates.rsm)")
-rk.print (waldtest.res)
+rk.print ("Call:")
+rk.print.literal (deparse(waldtest.res$call, width.cutoff=500))
+rk.header ("Wald test on item level (z-values):", level=4)
+rk.print(waldtest.res$coef.table)
})
-.rk.rerun.plugin.link(plugin="rkward::eRm_waldtest", settings="rad_splitcr.string=median\nx.available=estimates.rsm", label="Run again")
+.rk.rerun.plugin.link(plugin="rkward::eRm_waldtest", settings="drop_optimizer.string=optim\nrad_splitcr.string=median\nx.available=estimates.rsm", label="Run again")
.rk.make.hr()
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.Wald_test.rkout
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.Wald_test.rkout 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.Wald_test.rkout 2010-08-06 15:06:45 UTC (rev 2959)
@@ -1,6 +1,47 @@
<h1>Wald test (estimates.rsm)</h1>
DATE<br>
-<p class='character'><br>Wald test on item level (z-values):<br><br> z-statistic p-value<br>beta I1.c1 0.388 0.698<br>beta I1.c2 NaN NaN<br>beta I1.c3 NaN NaN<br>beta I6.c1 -0.388 0.698<br>beta I6.c2 NaN NaN<br>beta I6.c3 NaN NaN<br></p><h2>Messages, warnings, or errors:</h2>
+
+<p class='character'>Call:</p>
+<pre>Waldtest.Rm(object = estimates.rsm)</pre>
+<h4>Wald test on item level (z-values):</h4>
+<br>
+
+
+<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>z-statistic</th><th>p-value</th> </tr>
+
+<tr><td class=firstcolumn>beta I1.c1</td>
+<td class=cellinside> 0.39</td>
+<td class=cellinside> 0.70</td></tr>
+
+<tr><td class=firstcolumn>beta I1.c2</td>
+<td class=cellinside> 0.15</td>
+<td class=cellinside> 0.88</td></tr>
+
+<tr><td class=firstcolumn>beta I1.c3</td>
+<td class=cellinside> 0.16</td>
+<td class=cellinside> 0.88</td></tr>
+
+<tr><td class=firstcolumn>beta I6.c1</td>
+<td class=cellinside>-0.39</td>
+<td class=cellinside> 0.70</td></tr>
+
+<tr><td class=firstcolumn>beta I6.c2</td>
+<td class=cellinside> 0.14</td>
+<td class=cellinside> 0.89</td></tr>
+
+<tr><td class=firstcolumn>beta I6.c3</td>
+<td class=cellinside> 0.15</td>
+<td class=cellinside> 0.88</td></tr>
+
+ </tbody>
+</table>
+ </td></table>
+ <br>
+<h2>Messages, warnings, or errors:</h2>
<pre>I2 I3 I4 I5
Subgroup models are estimated without these items!
</pre>
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.goodnes-of-fit_Rasch.rkcommands.R
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.goodnes-of-fit_Rasch.rkcommands.R 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.goodnes-of-fit_Rasch.rkcommands.R 2010-08-06 15:06:45 UTC (rev 2959)
@@ -5,7 +5,12 @@
GoFRasch.res <- GoF.rasch(estimates.rasch)
## Print result
rk.header ("Goodness of Fit for Rasch Models (estimates.rasch)")
-rk.print (GoFRasch.res)
+rk.print ("Call:")
+rk.print.literal (deparse(GoFRasch.res$call, width.cutoff=500))
+rk.header ("Parametric Bootstrap test:", level=4)
+rk.print (paste("Chi-squared statistic:", round(GoFRasch.res$Tobs, digits=3)))
+rk.print (paste("Bootstrap samples:", GoFRasch.res$B))
+rk.print (paste("p-value:", GoFRasch.res$p.value))
})
.rk.rerun.plugin.link(plugin="rkward::ltm_gof_rasch", settings="spin_samples.real=49.00\nx.available=estimates.rasch", label="Run again")
.rk.make.hr()
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 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.goodnes-of-fit_Rasch.rkout 2010-08-06 15:06:45 UTC (rev 2959)
@@ -1,3 +1,13 @@
<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.94 <br></p>
\ No newline at end of file
+
+<p class='character'>Call:</p>
+<pre>rasch(data = LSAT)</pre>
+<h4>Parametric Bootstrap test:</h4>
+<br>
+
+<p class='character'>Chi-squared statistic: 18.345</p>
+
+<p class='character'>Bootstrap samples: 49</p>
+
+<p class='character'>p-value: 0.78</p>
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.item_fit_statistics.rkcommands.R
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.item_fit_statistics.rkcommands.R 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.item_fit_statistics.rkcommands.R 2010-08-06 15:06:45 UTC (rev 2959)
@@ -5,7 +5,12 @@
itemfit.res <- item.fit(estimates.rasch)
## Print result
rk.header ("Item-fit statistics (estimates.rasch)")
-rk.print (itemfit.res)
+rk.print ("Call:")
+rk.print.literal (deparse(itemfit.res$call, width.cutoff=500))
+rk.print ("Alternative: Items do not fit the model")
+rk.print (paste("Ability Categories:", itemfit.res$G))
+rk.header ("Item-Fit Statistics and P-values:", level=4)
+rk.print(cbind("X^2"=round(itemfit.res$Tobs, digits=3), "Pr (>X^2)"=format(round(itemfit.res$p.values, digits=3), nsmall=3)))
})
.rk.rerun.plugin.link(plugin="rkward::ltm_item_fit", settings="drop_sumgroups.string=median\nrad_pvalue.string=chi2\nspin_groups.real=10.00\nx.available=estimates.rasch", label="Run again")
.rk.make.hr()
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.item_fit_statistics.rkout
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.item_fit_statistics.rkout 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.item_fit_statistics.rkout 2010-08-06 15:06:45 UTC (rev 2959)
@@ -1,3 +1,43 @@
<h1>Item-fit statistics (estimates.rasch)</h1>
DATE<br>
-<p class='character'><br>Item-Fit Statistics and P-values<br><br>Call:<br>rasch(data = LSAT)<br><br>Alternative: Items do not fit the model<br>Ability Categories: 10<br><br> X^2 Pr(>X^2)<br>Item 1 61.93 <0.0001<br>Item 2 159.01 <0.0001<br>Item 3 233.79 <0.0001<br>Item 4 132.47 <0.0001<br>Item 5 84.76 <0.0001<br><br></p>
\ No newline at end of file
+
+<p class='character'>Call:</p>
+<pre>rasch(data = LSAT)</pre>
+
+<p class='character'>Alternative: Items do not fit the model</p>
+
+<p class='character'>Ability Categories: 10</p>
+<h4>Item-Fit Statistics and P-values:</h4>
+<br>
+
+
+<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>X^2</th><th>Pr (>X^2)</th> </tr>
+
+<tr><td class=firstcolumn>Item 1</td>
+<td class=cellinside>61.929 </td>
+<td class=cellinside>0.000 </td></tr>
+
+<tr><td class=firstcolumn>Item 2</td>
+<td class=cellinside>159.01 </td>
+<td class=cellinside>0.000 </td></tr>
+
+<tr><td class=firstcolumn>Item 3</td>
+<td class=cellinside>233.787</td>
+<td class=cellinside>0.000 </td></tr>
+
+<tr><td class=firstcolumn>Item 4</td>
+<td class=cellinside>132.473</td>
+<td class=cellinside>0.000 </td></tr>
+
+<tr><td class=firstcolumn>Item 5</td>
+<td class=cellinside>84.76 </td>
+<td class=cellinside>0.000 </td></tr>
+
+ </tbody>
+</table>
+ </td></table>
+ <br>
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.person_fit_statistics.rkcommands.R
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.person_fit_statistics.rkcommands.R 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.person_fit_statistics.rkcommands.R 2010-08-06 15:06:45 UTC (rev 2959)
@@ -5,7 +5,10 @@
personfit.res <- person.fit(estimates.rasch)
## Print result
rk.header ("Person-fit statistics (estimates.rasch)")
-rk.print (personfit.res)
+rk.print ("Call:")
+rk.print.literal (deparse(personfit.res$call, width.cutoff=500))
+rk.header ("Response patterns, person-fit statistics (L0, Lz) and probabilities for each response pattern (Pr):", level=4)
+rk.print(cbind(format(personfit.res$resp.patterns, nsmall=0), round(personfit.res$Tobs, digits=3), "Pr (<Lz)"=round(c(personfit.res$p.values), digits=3)))
})
.rk.rerun.plugin.link(plugin="rkward::ltm_person_fit", settings="rad_hypot.string=less\nrad_pvalue.string=normal\nrad_resppat.string=observed\nx.available=estimates.rasch", label="Run again")
.rk.make.hr()
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.person_fit_statistics.rkout
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.person_fit_statistics.rkout 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.person_fit_statistics.rkout 2010-08-06 15:06:45 UTC (rev 2959)
@@ -1,3 +1,319 @@
<h1>Person-fit statistics (estimates.rasch)</h1>
DATE<br>
-<p class='character'><br>Person-Fit Statistics and P-values<br><br>Call:<br>rasch(data = LSAT)<br><br>Alternative: Inconsistent response pattern under the estimated model<br><br> Item 1 Item 2 Item 3 Item 4 Item 5 L0 Lz Pr(<Lz)<br>1 0 0 0 0 0 -3.9914 -1.1186 0.1317<br>2 0 0 0 0 1 -3.9467 -1.0917 0.1375<br>3 0 0 0 1 0 -4.7396 -2.0187 0.0218<br>4 0 0 0 1 1 -4.1206 -1.3376 0.0905<br>5 0 0 1 0 0 -5.8062 -3.2656 0.0005<br>6 0 0 1 0 1 -5.1872 -2.4767 0.0066<br>7 0 0 1 1 0 -5.9802 -3.3235 0.0004<br>8 0 0 1 1 1 -4.7717 -2.0310 0.0211<br>9 0 1 0 0 0 -5.0475 -2.3786 0.0087<br>10 0 1 0 0 1 -4.4284 -1.6664 0.0478<br>11 0 1 0 1 1 -4.0129 -1.3093 0.0952<br>12 0 1 1 0 1 -5.0795 -2.3238 0.0101<br>13 0 1 1 1 0 -5.8725 -3.0780 0.001<br>14 0 1 1 1 1 -4.0395 -1.4551 0.0728<br>15 1 0 0 0 0 -3.3161 -0.3545 0.3615<br>16 1 0 0 0 1 -2.6970 0.1827 0.5725<br>17 1 0 0 1 0 -3.4900 -0.6641 0.2533<br>18 1 0 0 1 1 -2.2815 0.3375 0.6321<br>19 1 0 1 0 0 -4.5566 -1.8032 0.0357<br>20 1 0 1 0 1 -3.3481 -0.6769 0.2492<br>21 1 0 1 1 0 -4.1411 -1.4312 0.0762<br>22 1 0 1 1 1 -2.3081 0.0190 0.5076<br>23 1 1 0 0 0 -3.7978 -0.9929 0.1604<br>24 1 1 0 0 1 -2.5894 0.0447 0.5178<br>25 1 1 0 1 0 -3.3823 -0.7095 0.239<br>26 1 1 0 1 1 -1.5494 0.6651 0.747<br>27 1 1 1 0 0 -4.4489 -1.7240 0.0424<br>28 1 1 1 0 1 -2.6160 -0.2431 0.404<br>29 1 1 1 1 0 -3.4089 -0.9182 0.1793<br>30 1 1 1 1 1 -0.8945 0.8407 0.7997<br><br></p>
\ No newline at end of file
+
+<p class='character'>Call:</p>
+<pre>rasch(data = LSAT)</pre>
+<h4>Response patterns, person-fit statistics (L0, Lz) and probabilities for each response pattern (Pr):</h4>
+<br>
+
+
+<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>Item 1</th><th>Item 2</th><th>Item 3</th><th>Item 4</th><th>Item 5</th><th>L0</th><th>Lz</th><th>Pr (<Lz)</th> </tr>
+
+<tr><td class=firstcolumn>1</td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>-3.991</td>
+<td class=cellinside>-1.119</td>
+<td class=cellinside>0.132 </td></tr>
+
+<tr><td class=firstcolumn>2</td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>-3.947</td>
+<td class=cellinside>-1.092</td>
+<td class=cellinside>0.137 </td></tr>
+
+<tr><td class=firstcolumn>3</td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>-4.74 </td>
+<td class=cellinside>-2.019</td>
+<td class=cellinside>0.022 </td></tr>
+
+<tr><td class=firstcolumn>4</td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>-4.121</td>
+<td class=cellinside>-1.338</td>
+<td class=cellinside>0.091 </td></tr>
+
+<tr><td class=firstcolumn>5</td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>-5.806</td>
+<td class=cellinside>-3.266</td>
+<td class=cellinside>0.001 </td></tr>
+
+<tr><td class=firstcolumn>6</td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>-5.187</td>
+<td class=cellinside>-2.477</td>
+<td class=cellinside>0.007 </td></tr>
+
+<tr><td class=firstcolumn>7</td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>-5.98 </td>
+<td class=cellinside>-3.324</td>
+<td class=cellinside>0 </td></tr>
+
+<tr><td class=firstcolumn>8</td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>-4.772</td>
+<td class=cellinside>-2.031</td>
+<td class=cellinside>0.021 </td></tr>
+
+<tr><td class=firstcolumn>9</td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>-5.047</td>
+<td class=cellinside>-2.379</td>
+<td class=cellinside>0.009 </td></tr>
+
+<tr><td class=firstcolumn>10</td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>-4.428</td>
+<td class=cellinside>-1.666</td>
+<td class=cellinside>0.048 </td></tr>
+
+<tr><td class=firstcolumn>11</td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>-4.013</td>
+<td class=cellinside>-1.309</td>
+<td class=cellinside>0.095 </td></tr>
+
+<tr><td class=firstcolumn>12</td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>-5.08 </td>
+<td class=cellinside>-2.324</td>
+<td class=cellinside>0.01 </td></tr>
+
+<tr><td class=firstcolumn>13</td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>-5.872</td>
+<td class=cellinside>-3.078</td>
+<td class=cellinside>0.001 </td></tr>
+
+<tr><td class=firstcolumn>14</td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>-4.04 </td>
+<td class=cellinside>-1.455</td>
+<td class=cellinside>0.073 </td></tr>
+
+<tr><td class=firstcolumn>15</td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>-3.316</td>
+<td class=cellinside>-0.354</td>
+<td class=cellinside>0.361 </td></tr>
+
+<tr><td class=firstcolumn>16</td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>-2.697</td>
+<td class=cellinside>0.183 </td>
+<td class=cellinside>0.572 </td></tr>
+
+<tr><td class=firstcolumn>17</td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>-3.49 </td>
+<td class=cellinside>-0.664</td>
+<td class=cellinside>0.253 </td></tr>
+
+<tr><td class=firstcolumn>18</td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>-2.282</td>
+<td class=cellinside>0.338 </td>
+<td class=cellinside>0.632 </td></tr>
+
+<tr><td class=firstcolumn>19</td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>-4.557</td>
+<td class=cellinside>-1.803</td>
+<td class=cellinside>0.036 </td></tr>
+
+<tr><td class=firstcolumn>20</td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>-3.348</td>
+<td class=cellinside>-0.677</td>
+<td class=cellinside>0.249 </td></tr>
+
+<tr><td class=firstcolumn>21</td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>-4.141</td>
+<td class=cellinside>-1.431</td>
+<td class=cellinside>0.076 </td></tr>
+
+<tr><td class=firstcolumn>22</td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>-2.308</td>
+<td class=cellinside>0.019 </td>
+<td class=cellinside>0.508 </td></tr>
+
+<tr><td class=firstcolumn>23</td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>-3.798</td>
+<td class=cellinside>-0.993</td>
+<td class=cellinside>0.16 </td></tr>
+
+<tr><td class=firstcolumn>24</td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>-2.589</td>
+<td class=cellinside>0.045 </td>
+<td class=cellinside>0.518 </td></tr>
+
+<tr><td class=firstcolumn>25</td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>-3.382</td>
+<td class=cellinside>-0.709</td>
+<td class=cellinside>0.239 </td></tr>
+
+<tr><td class=firstcolumn>26</td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>-1.549</td>
+<td class=cellinside>0.665 </td>
+<td class=cellinside>0.747 </td></tr>
+
+<tr><td class=firstcolumn>27</td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>-4.449</td>
+<td class=cellinside>-1.724</td>
+<td class=cellinside>0.042 </td></tr>
+
+<tr><td class=firstcolumn>28</td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>-2.616</td>
+<td class=cellinside>-0.243</td>
+<td class=cellinside>0.404 </td></tr>
+
+<tr><td class=firstcolumn>29</td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>0 </td>
+<td class=cellinside>-3.409</td>
+<td class=cellinside>-0.918</td>
+<td class=cellinside>0.179 </td></tr>
+
+<tr><td class=firstcolumn>30</td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>1 </td>
+<td class=cellinside>-0.895</td>
+<td class=cellinside>0.841 </td>
+<td class=cellinside>0.8 </td></tr>
+
+ </tbody>
+</table>
+ </td></table>
+ <br>
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.unidimensional.rkcommands.R
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.unidimensional.rkcommands.R 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.unidimensional.rkcommands.R 2010-08-06 15:06:45 UTC (rev 2959)
@@ -5,7 +5,16 @@
unidim.res <- unidimTest(estimates.rasch)
## Print result
rk.header ("Unidimensionality check (estimates.rasch)")
-rk.print (unidim.res)
+rk.print ("Call:")
+rk.print.literal (deparse(unidim.res$call, width.cutoff=500))
+rk.header ("Matrix of tertachoric correlations:", level=4)
+rk.print (unidim.res$Rho)
+rk.header ("Unidimensionality Check using Modified Parallel Analysis:", level=4)
+rk.print ("Alternative hypothesis: <em>The second eigenvalue of the observed data is substantially larger than the second eigenvalue of data under the assumed IRT model</em>")
+rk.print (paste("Second eigenvalue in the observed data:", round(unidim.res$Tobs[2], digits=3)))
+rk.print (paste("Average of second eigenvalues in Monte Carlo samples:", round(mean(unidim.res$T.boot[,2]), digits=3)))
+rk.print (paste("Monte Carlo samples:", dim(unidim.res$T.boot)[1]))
+rk.print (paste("p-value:", round(unidim.res$p.value, digits=3)))
# keep results in current workspace
.GlobalEnv$unidim.res <- unidim.res
})
Modified: trunk/rkward/tests/item_response_theory/RKTestStandard.unidimensional.rkout
===================================================================
--- trunk/rkward/tests/item_response_theory/RKTestStandard.unidimensional.rkout 2010-08-02 09:58:43 UTC (rev 2958)
+++ trunk/rkward/tests/item_response_theory/RKTestStandard.unidimensional.rkout 2010-08-06 15:06:45 UTC (rev 2959)
@@ -1,3 +1,66 @@
<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.254<br>Monte Carlo samples: 100<br>p-value: 0.673<br></p>
\ No newline at end of file
+
+<p class='character'>Call:</p>
+<pre>rasch(data = LSAT)</pre>
+<h4>Matrix of tertachoric correlations:</h4>
+<br>
+
+
+<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>Item 1</th><th>Item 2</th><th>Item 3</th><th>Item 4</th><th>Item 5</th> </tr>
+
+<tr><td class=firstcolumn>Item 1</td>
+<td class=cellinside>1.000</td>
+<td class=cellinside>0.170</td>
+<td class=cellinside>0.228</td>
+<td class=cellinside>0.107</td>
+<td class=cellinside>0.067</td></tr>
+
+<tr><td class=firstcolumn>Item 2</td>
+<td class=cellinside>0.170</td>
+<td class=cellinside>1.000</td>
+<td class=cellinside>0.189</td>
+<td class=cellinside>0.111</td>
+<td class=cellinside>0.172</td></tr>
+
+<tr><td class=firstcolumn>Item 3</td>
+<td class=cellinside>0.228</td>
+<td class=cellinside>0.189</td>
+<td class=cellinside>1.000</td>
+<td class=cellinside>0.187</td>
+<td class=cellinside>0.105</td></tr>
+
+<tr><td class=firstcolumn>Item 4</td>
+<td class=cellinside>0.107</td>
+<td class=cellinside>0.111</td>
+<td class=cellinside>0.187</td>
+<td class=cellinside>1.000</td>
+<td class=cellinside>0.201</td></tr>
+
+<tr><td class=firstcolumn>Item 5</td>
+<td class=cellinside>0.067</td>
+<td class=cellinside>0.172</td>
+<td class=cellinside>0.105</td>
+<td class=cellinside>0.201</td>
+<td class=cellinside>1.000</td></tr>
+
+ </tbody>
+</table>
+ </td></table>
+ <br>
+<h4>Unidimensionality Check using Modified Parallel Analysis:</h4>
+<br>
+
+<p class='character'>Alternative hypothesis: <em>The second eigenvalue of the observed data is substantially larger than the second eigenvalue of data under the assumed IRT model</em></p>
+
+<p class='character'>Second eigenvalue in the observed data: 0.225</p>
+
+<p class='character'>Average of second eigenvalues in Monte Carlo samples: 0.255</p>
+
+<p class='character'>Monte Carlo samples: 100</p>
+
+<p class='character'>p-value: 0.634</p>
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