parallelpam (version 1.0)

CalcAndWriteDissimilarityMatrix: CalcAndWriteDissimilarityMatrix

Description

Writes a binary symmetric matrix with the dissimilarities between ROWS of the data stored in a binary matrix in the scellpam package format.
Notice that, differently from the common practice in single cell, the rows represent cells. This is for efficiency reasons and it is transparent to the user, as long as he/she has generated the binary matrix (with CsvToBinMat, dgCMatToBinMat or SceToBinMat) using the option transpose=TRUE.
The input matrix of vectors can be a full or a sparse matrix. Output matrix type can be float or double type (but look at the comments in 'Details').

Usage

CalcAndWriteDissimilarityMatrix(
  ifname,
  ofname,
  distype = "L2",
  restype = "float",
  comment = "",
  nthreads = 0L
)

Value

No return value, called for side effects (creates a file)

Arguments

ifname

A string with the name of the file containing the counts as a binary matrix, as written by CsvToBinMat, dgCMatToBinMat or SceToBinMat

ofname

A string with the name of the binary output file to contain the symmetric dissimilarity matrix.

distype

The dissimilarity to be calculated. It must be one of these strings: 'L1', 'L2' or 'Pearson'.
Default: 'L2'.

restype

The data type of the result. It can be one of the strings 'float' or 'double'. Default: float (and don't change it unless you REALLY need to...).

comment

Comment to be added to the dissimilary matrix. Default: "" (no comment)

nthreads

Number of threads to be used for the parallel calculations with this meaning:
-1: don't use threads.
0: let the function choose according to the number of individuals (cells) and to the number of available cores.
Any possitive number > 1: use that number of threads. You can use even more than cores, but this is discouraged and raises a warning.
Default: 0.

Details

The parameter restype forces the output to be a matrix of either floats or doubles. Precision of float if normally good enough; but if you need double precision (may be because you expect your results to be in a large range, two to three orders of magnitude), change it.
Nevertheless, notice that this at the expense of double memory usage, which is QUADRATIC with the number of individuals (rows) in your input matrix.

Examples

Rf <- matrix(runif(50000),nrow=100)
JWriteBin(Rf,"Rfullfloat.bin",dtype="float",dmtype="full",
          comment="Full matrix of floats, 100 rows, 500 columns")
JMatInfo("Rfullfloat.bin")
CalcAndWriteDissimilarityMatrix("Rfullfloat.bin","RfullfloatDis.bin",distype="L2",
                        restype="float",comment="L2 distance matrix from full",nthreads=0)
JMatInfo("RfullfloatDis.bin")
JWriteBin(Rf,"Rsparsefloat.bin",dtype="float",dmtype="sparse",
                        comment="Sparse matrix of floats, 100 rows, 500 columns")
JMatInfo("Rsparsefloat.bin")
CalcAndWriteDissimilarityMatrix("Rsparsefloat.bin","RsparsefloatDis.bin",distype="L2",
                        restype="float",comment="L2 distance matrix from sparse",nthreads=0)
JMatInfo("RsparsefloatDis.bin")
Dfu<-GetJManyRows("RfullfloatDis.bin",c(1:nrow(Rf)))
Dsp<-GetJManyRows("RsparsefloatDis.bin",c(1:nrow(Rf)))
max(Dfu-Dsp)
file.remove("Rfullfloat.bin")
file.remove("Rsparsefloat.bin")
file.remove("RfullfloatDis.bin")
file.remove("RsparsefloatDis.bin")