1 ppam_exe = function(path=".", np=parallel::detectCores(), data=NULL,
2 args="DontLetMeEmptyPlease!")
4 command_line = paste("mpirun -np ",np," ",path,"/ppam.exe cluster",sep="")
6 #if data provided (as data.frame or matrix...): binarize it, and add it as first argument
9 if (!is.character(data))
11 #assuming matrix or data.frame, WITH row names
12 #( identifiers; could be line number... e.g. data <- cbind(1:nrow(data),data) )
13 write.table(data, "/tmp/data_csv", sep=",", row.names=FALSE, col.names=FALSE)
14 system(paste(path,"/ppam.exe serialize /tmp/data_csv /tmp/data_bin 0 0",sep=""))
17 system(paste(path,"/ppam.exe serialize ",data," /tmp/data_bin 0 0",sep=""))
19 command_line = paste(command_line," /tmp/data_bin",sep="")
22 command_line = paste(command_line," ",args,sep="")
26 #NOTE: identifiers in first column
27 getMedoids = function(path=".", xmlResult = "ppamResult.xml",
28 finalSeries = "ppamFinalSeries.bin")
30 system(paste(path,"/ppam.exe deserialize ",finalSeries," ppamFinalSeries.csv -1",sep=""))
31 curves = read.table("ppamFinalSeries.csv", sep=",")
33 ranks = as.integer( xmlToList( xmlParse(xmlResult) )$ranks )
34 return ( curves[ranks,] ) # == medoids
37 #TODO: check C function( is it correct?!)
38 getDistor = function(path=".", xmlResult = "ppamResult.xml",
39 finalSeries = "ppamFinalSeries.bin")
41 system(paste(path,"/ppam.exe classif ",finalSeries," ",xmlResult,sep=""))