self.pathout = analyse.pathout
self.analyse = analyse
self.parametres = analyse.parametres
+ #self.scriptout = self.pathout['lastRscript.R']
self.scriptout = self.pathout['temp']
self.script = u"#Script genere par IRaMuTeQ - %s\n" % datetime.now().ctime()
-
+
def add(self, txt) :
self.script = '\n'.join([self.script, txt])
-
+
def defvar(self, name, value) :
self.add(' <- '.join([name, value]))
""" % ffr(DicoPath['TableUc2'])
txt += """
log1 <- "%s"
- chd1<-CHD(data1, x = nbt, mode.patate = mode.patate, svd.method = svd.method, libsvdc.path = libsvdc.path, log.file = log1)
+ chd1<-CHD(data1, x = nbt, mode.patate = mode.patate, svd.method =
+ svd.method, libsvdc.path = libsvdc.path, log.file = log1)
""" % ffr(DicoPath['log-chd1.txt'])
if classif_mode == 0:
txt += """
log2 <- "%s"
- chd2<-CHD(data2, x = nbt, mode.patate = mode.patate, svd.method = svd.method, libsvdc.path = libsvdc.path, log.file = log2)
+ chd2<-CHD(data2, x = nbt, mode.patate = mode.patate, svd.method =
+ svd.method, libsvdc.path = libsvdc.path) log.file = log2)
""" % ffr(DicoPath['log-chd2.txt'])
txt += """
seuil <- %s
if (!is.null(seuil)) {
if (method!='cooc') {
- seuil <- seuil/100
+ seuil <- seuil/1000
}
}
""" % seuil
self.pathout = PathOut()
self.parametres = parametres
self.scriptout = self.pathout['temp']
-
+
def make_script(self) :
#FIXME
-
+
txt = """
library(igraph)
library(Matrix)
RData = os.path.join(path,'RData.RData')
txt += load % (ffr(RData), gname)
self.add(txt)
- self.sources(['/home/pierre/workspace/iramuteq/Rscripts/simi.R'])
+ self.sources([self.analyse.parent.RscriptsPath['simi']])
txt = """
ng <- merge.graph(graphs)
ngraph <- list(graph=ng, layout=layout.fruchterman.reingold(ng, dim=3), labex.cex=V(ng)$weight)
tgen <- read.csv2("%s", row.names = 1, sep = '\\t')
""" % ffr(self.parametres['tgeneff'])
txt += """
+ tgenlem <- read.csv2("%s", row.names = 1, sep = '\\t')
+ """ % ffr(self.parametres['tgenlemeff'])
+ txt += """
res <- build.prof.tgen(tgen)
write.table(res$chi2, file = "%s", sep='\\t', col.names = NA)
write.table(res$pchi2, file = "%s", sep='\\t', col.names = NA)
""" % (ffr(self.pathout['tgenchi2.csv']), ffr(self.pathout['tgenpchi2.csv']))
+ txt += """
+ reslem <- build.prof.tgen(tgenlem)
+ write.table(reslem$chi2, file = "%s", sep='\\t', col.names = NA)
+ write.table(reslem$pchi2, file = "%s", sep='\\t', col.names = NA)
+ """ % (ffr(self.pathout['tgenlemchi2.csv']), ffr(self.pathout['tgenlempchi2.csv']))
self.add(txt)
class FreqMultiScript(PrintRScript):
dev.off()
""" % ffr(self.pathout['barplotrow.png'])
self.add(txt)
- self.write()
\ No newline at end of file
+ self.write()