\name{psych.misc} \alias{psych.misc} \alias{misc} \alias{tableF} \alias{lowerCor} \alias{lowerMat} \alias{matMult} \alias{progressBar} \alias{reflect} \alias{shannon} \alias{test.all} \alias{cor2} \alias{levels2numeric} \alias{char2numeric} \alias{nchar2numeric} \alias{isCorrelation} \alias{isCovariance} \alias{fromTo} \alias{cs} \alias{acs} \alias{SAPAfy} \title{Miscellaneous helper functions for the psych package} \description{This is a set of minor, if not trivial, helper functions. lowerCor finds the correlation of x variables and then prints them using lowerMat which is a trivial, but useful, function to round off and print the lower triangle of a matrix. reflect reflects the output of a factor analysis or principal components analysis so that one or more factors is reflected. (Requested by Alexander Weiss.) progressBar prints out ... as a calling routine (e.g., \code{\link{tetrachoric}}) works through a tedious calculation. shannon finds the Shannon index (H) of diversity or of information. test.all tests all the examples in a package. best.items sorts a factor matrix for absolute values and displays the expanded items names. fa.lookup returns sorted factor analysis output with item labels. \code{\link{cor2}} correlates two data.frames (of equal length). levels2numeric and char2numeric convert dataframe columns that are categorical/levels to numeric values. } \usage{ psych.misc() lowerCor(x,digits=2,use="pairwise",method="pearson",minlength=5,cor="cor",show=TRUE) cor2(x,y,digits=2,use="pairwise",method="pearson",cor="cor",show=TRUE) lowerMat(R, digits = 2,minlength=5) matMult(x,y) #multiply two matrices with missing data tableF(x,y) reflect(f,flip=NULL) progressBar(value,max,label=NULL) shannon(x,correct=FALSE,base=2) test.all(pl,package="psych",dependencies = c("Depends", "Imports", "LinkingTo"),find=FALSE,skip=NULL) levels2numeric(x) char2numeric(x,flag=TRUE) nchar2numeric(x,flag=TRUE) isCorrelation(x,na.rm=FALSE) #test if an object is a symmetric matrix # with diagonals of 1 and all values between -1 and 1 isCovariance(x) #test if an object is a symmetric matrix fromTo(data,from,to=NULL) #convert character names to locations as specified in colnames #of data cs(...) #convert a list of text words to character vector acs(...) #convert a list of text words to a single string SAPAfy(x,y) #sample y columns from x, replace other columns with NA } \arguments{ \item{R}{A rectangular matrix or data frame (probably a correlation matrix)} \item{x}{A data matrix or data frame or a vector depending upon the function.} \item{y}{A data matrix or data frame or a vector} \item{f}{The object returned from either a factor analysis (fa) or a principal components analysis (principal) } \item{digits}{round to digits} \item{minlength}{Abbreviate to minlength in lowerCor} \item{show}{Display the correlations from lowerCor or cor2 (default is TRUE)} \item{use}{Should pairwise deletion be done, or one of the other options to cor} \item{na.rm}{Should we check for NA on the diagonal of a correlation matices} \item{method}{"pearson", "kendall", "spearman"} \item{cor}{defaults to the normal cor function, but can also do tetrachoric, polychoric or covariances (cov)} \item{value}{the current value of some looping variable} \item{max}{The maximum value the loop will achieve} \item{label}{what function is looping} \item{flip}{The factor or components to be reversed keyed (by factor number)} \item{flag}{flag=TRUE in char2numeric will flag (with *) those variables that had been numeric. This changes the variable name. flag = FALSE does not mark those variables.} \item{correct}{Correct for the maximum possible information in this item} \item{base}{What is the base for the log function (default=2, e implies base = exp(1))} \item{pl}{The name of a package (or list of packages) to be activated and then have all the examples tested.} \item{package}{Find the dependencies for this package, e.g., psych} \item{dependencies}{Which type of dependency to examine?} \item{find}{Look up the dependencies, and then test all of their examples} \item{skip}{Do not test these dependencies} \item{data}{A dataframe or matrix to choose from} \item{from}{select from column with name from to column with name to} \item{to}{select from column from to column to} \item{...}{Any string of legitimate objects} } \value{ \code{\link{tableF}} is fast alternative to the table function for creating two way tables of numeric variables. It does not have any of the elegant checks of the table function and thus is much faster. Used in the \code{\link{tetrachoric}} and \code{\link{polychoric}} functions to maximize speed. \code{\link{lowerCor}} Finds and prints (using \code{\link{lowerMat}}) the lower diagonal correlation matrix but returns (invisibly) the full correlation matrix found with the use, method and cor parameters. The default values are for pairwise deletion of variables, "pearson" correlations. Specify cor="tetrachoric", "polychoric", or "mixed" for alternatives. By default, digits is set to print to 2 decimal places. By default, it will change character variables to numeric and flag them. \code{\link{lowerMat}}Shows the lower triangle of a matrix, rounded to digits with titles abbreviated to digits + 3 \code{\link{progressBar}} Display a series of dots as we progress through a slow loop (removed from anything using multicores). \code{\link{tableF}} (for tableFast) is a cut down version of table that does no error checking, nor returns pretty output, but is significantly faster than table. It will just work on two integer vectors. This is used in polychoric an tetrachoric for about a 50\% speed improvement for large problems. \code{\link{shannon}} finds Shannon's H index of information. Used for estimating the complexity or diversity of the distribution of responses in a vector or matrix. \deqn{H = -\sum{p_i log(p_i) }} \code{\link{test.all}} allows one to test all the examples in specified package. This allows us to make sure that those examples work when other packages (e.g., psych) are also loaded. This is used when developing revisions to the psych package to make sure the the other packages work. Some packages will not work and/or crash the system (e.g., DeducerPlugInScaling requires Java and even with Java, crashes when loaded, even if psych is not there!). Alternatively, if testing a long list of dependencies, you can skip the first part by specifying them by name. \code{\link{cor2}} will find and display the correlations between two sets of variables, rounded to digits, using the other options. If x is a list of multiple sets (two or more), then all sets are correlated. \code{\link{levels2numeric}}converts character data with levels to numeric data. Used in the SAPA analyses where we code some variables, (e.g., gender, education) with character codes to help in the documentation of files, but want to do analyses of correlations with other categorical variables. \code{\link{char2numeric}}converts character data with levels to numeric data. Used for cases when data from questionnaires include the response categories rathere than numeric data. Unless the levels of the data are in meaningful order, the numeric results are not useful. Most useful if doing polychoric analyses. Note this is not suitable for recoding numeric data stored as characters, for it will force them to levels first. See \code{\link{nchar2numeric}}. Note that for very small data sets, because the recoding done by char2numeric is column wise, it might result in different numerical results for the same characters in different columns. Problems with \code{\link{char2numeric}} can be solved by using the \code{psychTools::recode} function. \code{\link{nchar2numeric}}converts numbers coded as characters (quoted) to numeric without forcing them to factors first. \code{\link{fromTo}} selects the columns in data from to (see the examples) \code{\link{cs}} concatenates strings without the need to identify variables by " ". } \details{ \code{\link{lowerCor}} prints out the lower off diagonal matrix rounded to digits with column names abbreviated to digits + 3 characters, but also returns the full and unrounded matrix. By default, it uses pairwise deletion of variables. It in turn calls \code{\link{lowerMat}} which does the pretty printing. It is important to remember to not call \code{\link{lowerCor}} when all you need is \code{\link{lowerMat}}! \code{\link{cs}} is a direct copy of the Cs function in the Hmisc package by Frank Harrell. Added to psych to avoid the overhead of the Hmisc package. } \seealso{\code{\link{corr.test}} to find correlations, count the pairwise occurrences, and to give significance tests for each correlation. \code{\link{r.test}} for a number of tests of correlations, including tests of the difference between correlations. \code{\link{lowerUpper}} will display the differences between two matrices.} \examples{ lowerMat(Thurstone) lb <- lowerCor(psychTools::bfi[1:10]) #finds and prints the lower correlation matrix, # returns the square matrix. lowerCor(psychTools::ability[,1:5]) #the Pearson correlations, compare with lowerCor(psychTools::ability[,1:5],cor="tetra") #fiml <- corFiml(psychTools::bfi[1:10]) #FIML correlations require lavaan package #lowerMat(fiml) #to get pretty output f3 <- fa(Thurstone,3) f3r <- reflect(f3,2) #reflect the second factor #find the complexity of the response patterns of the iqitems. round(shannon(psychTools::iqitems),2) #test.all('BinNor') #Does the BinNor package work when we are using other packages bestItems(lb,"A3",cut=.1,dictionary=psychTools::bfi.dictionary[1:2]) #to make this a latex table #df2latex(bestItems(lb,2,cut=.2)) # data(psychTools::bfi.dictionary) f2 <- fa(psychTools::bfi[1:10],2) fa.lookup(f2,psychTools::bfi.dictionary) sa1 <-sat.act[1:2] sa2 <- sat.act[3:4] sa3 <- sat.act[5:6] cor2(sa1,sa2) cor2(list(sa1,sa2)) #show within set and between set cors cor2(list(sa1,sa2,sa3)) lowerCor(fromTo(sat.act,"ACT","SATQ")) #show some correlations vect <- cs(ACT,SATQ) #skip the quotes vect #they are in this vector #to combine longer terms vect <- cs("Here is a longish",vector, that, we ,"want to combine", into, several) vect temp <- acs("Here is a longish",vector, that, we ,"want to combine", into, one) temp lowerCor(fromTo(sat.act,cs(ACT,SATQ))) lowerCor(fromTo(bfi,cs(A3,C4)),cor="poly") #compare with lowerCor(fromTo(bfi,cs(A3,C4))) set.seed(42) temp <- SAPAfy(bfi[1:10],3) #30 % sample from bfi f2 <- fa(bfi[1:10],2) f2s <- fa(temp,2) fa.congruence(f2s,f2) #the two factor structure are almost identical #although the scores are not identical cor2(f2$scores, f2s$scores) } % Add one or more standard keywords, see file 'KEYWORDS' in the % R documentation directory. \keyword{multivariate}