Displaying 1 result from an estimated 1 matches for "rplink".
Did you mean:
uplink
2011 Aug 30
2
Error in evalauating a function
...2, 0, 1, 1, 2, 0, 0, 1, 0, 1, -1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 1, 0, 1, 0, 1, 0, 2, 0, 1, 1, 0, 0, 1, 1, 1, 0, 0, 0, 0, 1, 1, 2, 0, 0, 0, 1, 1, 1, 0, 0 )
GENO <- matrix( g , nrow = n ,byrow=T)
GENO[GENO == -1 ] <- NA?
# Function
library(ordinal)
???? Rplink <- function(PHENO,GENO,CLUSTER,COVAR)
???? {
????? f1 <- function(s)
????? {
??????? m <- summary(clm(PHENO ~ s))
??????? r <- c(m$coef)
??????? c( length(r) , r )
????? }
???? apply( GENO , 2 , f1 )
???? }
Thanks,
?Debs