# Phi coefficient matrix (package psych) Classic List Threaded 4 messages Open this post in threaded view
|

## Phi coefficient matrix (package psych)

 I want to calculate phi coefficient for every pair of the columns. Is there a way to generate a matrix like a correlation matrix? I know cor function in the case below gives same answer as phi coefficient. ​x <- sample(c(0,1), 10, replace=TRUE) y <- sample(c(0,1), 10, replace=TRUE) z <- sample(c(0,1), 10, replace=TRUE) df <- data.frame(x,y,z) cor(df) library(psych) phi(df)​ Thank you, Kumar Mainali Postdoctoral Associate Department of Biology University of Maryland ᐧ         [[alternative HTML version deleted]] ______________________________________________ [hidden email] mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.
Open this post in threaded view
|

## Re: Phi coefficient matrix (package psych)

 HI Kumar, A simple way is: phimat<-function(x) {  xcol<-dim(x)  newx<-matrix(NA,nrow=xcol,ncol=xcol)  for(i in 1:xcol) {   for(j in 1:xcol) newx[i,j]<-phi(table(x[,i],x[,j]))  }  rownames(newx)<-colnames(newx)<-colnames(x)  return(newx) } phimat(df) Jim On Wed, Apr 22, 2015 at 6:34 AM, Kumar Mainali <[hidden email]> wrote: > I want to calculate phi coefficient for every pair of the columns. Is there > a way to generate a matrix like a correlation matrix? I know cor function > in the case below gives same answer as phi coefficient. > > x <- sample(c(0,1), 10, replace=TRUE) > y <- sample(c(0,1), 10, replace=TRUE) > z <- sample(c(0,1), 10, replace=TRUE) > df <- data.frame(x,y,z) > cor(df) > library(psych) > phi(df) > > Thank you, > Kumar Mainali > Postdoctoral Associate > Department of Biology > University of Maryland > ᐧ > >         [[alternative HTML version deleted]] > > ______________________________________________ > [hidden email] mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html> and provide commented, minimal, self-contained, reproducible code. ______________________________________________ [hidden email] mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.