"Yule2poly.matrix" <- function(x,v) { sizex <- dim(x)[2] if (!is.vector(v)) v <- as.vector(v) nv <- length(v) sizey <- dim(x)[1] for (i in 1:sizex) { xc <- v[i] for (j in 1:i) { yc <- v[j] if(x[i,j] > .999) {x[i,j] <- 1} else {x[i,j] <- Yule2poly(x[i,j],c(xc,yc))} x[j,i] <- x[i,j] } } return(x) } "phi2poly.matrix" <- function(x,v) { sizex <- dim(x)[2] if (!is.vector(v)) v <- as.vector(v) nv <- length(v) sizey <- dim(x)[1] for (i in 1:sizex) { xc <- v[i] for (j in 1:i) { yc <- v[j] if(x[i,j] > .999) {x[i,j] <- 1} else {x[i,j] <- phi2poly(x[i,j],xc,yc)} x[j,i] <- x[i,j] } } return(x) } "Yule2phi.matrix" <- function(x,v) { sizex <- dim(x)[2] if (!is.vector(v)) v <- as.vector(v) nv <- length(v) sizey <- dim(x)[1] for (i in 1:sizex) { xc <- v[i] for (j in 1:i) { yc <- v[j] if(x[i,j] > .999) {x[i,j] <- 1} else {x[i,j] <- Yule2phi(x[i,j],c(xc,yc))} x[j,i] <- x[i,j] } } return(x) } #revised August 29, 2010 "poly.mat" <- function(x,short=TRUE,std.err=FALSE,ML=FALSE) { .Deprecated("polychoric",msg="poly.mat is deprecated. Please use the polychoric function instead.") return(polychoric(x)) } "phi2polychor.matrix" <- function(x,v) { sizex <- dim(x)[2] if (!is.vector(v)) v <- as.vector(v) nv <- length(v) sizey <- dim(x)[1] for (i in 1:sizex) { xc <- v[i] for (j in 1:i) { yc <- v[j] if(x[i,j] > .999) {x[i,j] <- 1} else {x[i,j] <- phi2poly(x[i,j],xc,yc)} x[j,i] <- x[i,j] } } return(x) }