Attachment | Size |
---|---|

twinAceOrd.R.txt [6] | 7.36 KB |

Hello everyone,

I am trying to learn how to fit a univariate ACE model when the variable of interest is ordinal. For this I've started with the script called Univariate Twin Ordinal-Matrix/twinAceOrd.R that I found under the TC 2012 - OpenMx website. When I run the code i get the following four error messages:

1st,

> # Generate Descriptive Statistics

> #colMeans(mzDataOrd,na.rm=TRUE)

> #colMeans(dzDataOrd,na.rm=TRUE)

> cov(mzDataOrd,use="complete")

Error: is.numeric(x) || is.logical(x) is not TRUE

> cov(dzDataOrd,use="complete")

Error: is.numeric(x) || is.logical(x) is not TRUE

2nd,

> ACEpathMatrices <- c("a","c",path"e","iSD","iSD %*% a","iSD %*% c","iSD %*% e")

Error: unexpected string constant in "ACEpathMatrices <- c("a","c",path"e""

3rd,

> ACEpathLabels <- c("path_a","path_c","path_e","isd","spath_a","spath_c","spath_e")

> formatOutputMatrices(AceFit,ACEpathMatrices,ACEpathLabels,Vars,4)

[1] "Matrix ACE.a"

path_a1

bmi 0.8646

[1] "Matrix ACE.c"

path_c1

bmi 0.0000

[1] "Matrix ACE.e"

path_e1

bmi -0.5025

[1] "Matrix ACE.iSD"

Error in print(formatOutputMatrix(evalQuote(matricesList[[k]], fittedModel), :

error in evaluating the argument 'x' in selecting a method for function 'print': Error in eval(expr, envir, enclos) : object 'ACE.iSD' not found

4th,

> round(rbind(AceFit@output$estimate,AeFit@output$estimate,CeFit@output$estimate,eFit@output$estimate),4)

a11 c11 e11 t1Z t2Z t3Z t4Z

[1,] 0.8646 0.0000 -0.5025 -1.2212 0.8049 0.7102 0.5686

[2,] 0.8646 -0.5025 -1.2212 0.8049 0.7102 0.5687 0.8646

[3,] -0.7677 -0.6408 -1.2177 0.8078 0.7090 0.5719 -0.7677

[4,] -0.7677 -0.6408 -1.2177 0.8078 0.7090 0.5719 -0.7677

Warning message:

In rbind(AceFit@output$estimate, AeFit@output$estimate, CeFit@output$estimate, :

number of columns of result is not a multiple of vector length (arg 2)

I managed to fix the problem and get rid of the 2nd error message, but haven't really understood what the other errors are about.

I would be grateful for any kind of help. I am attching the script that I am using, where I haven't saved any changes.

/Natalja