#Author: Pierre Ratinaud
#Copyright (c) 20010-2013 Pierre Ratinaud
-#Lisense: GNU/GPL
+#License: GNU/GPL
#fichier genere par IRaMuTeq
source('%s', encoding = 'utf8')
typegraph <- %i
+edgesfile <- "%s"
+nodesfile <- "%s"
what <- %i
x <- %i
y <- %i
cex.txt <- %s
txt.min <- %i
txt.max <- %i
-fileout <- '%s'
+fileout <- "%s"
width <- %i
height <- %i
taillecar <- %i
tchi <- %s
tchi.min <- %i
tchi.max <- %i
-dirout <- '%s'
+dirout <- "%s"
do.svg <- %s
xminmax <- NULL
yminmax <- NULL
if (exists('afctable')) {
eff <- rowSums(afctable)
} else {
- eff <- afctable$rowmass
+ eff <- afc$rowmass
}
if (!is.null(debsup)) {
cex.par <- eff[debsup:(debet-1)]
}
if ( qui == 2 ) {
+ fin <- nrow(table.in)
table.in <- table.in[debet:nrow(table.in),]
tablechi <- tablechi[debet:nrow(tablechi),]
- cex.par <- eff[debet:nrow(afctable)]
+ cex.par <- eff[debet:fin]
}
}
cex.par <- eff
}
} else {
+ fin <- nrow(table.in)
table.in <- table.in[debet:nrow(table.in),]
tablechi <- tablechi[debet:nrow(tablechi),]
- cex.par <- eff[debet:nrow(afctable)]
+ cex.par <- eff[debet:fin]
}
}
table.in <- table.in[row.keep,]
tablechi <- tablechi[row.keep,]
} else if (do.select.chi.classe) {
- row.keep <- select.chi.classe(tablechi, ptbycluster)
+ row.keep <- select.chi.classe(tablechi, ptbycluster, active=FALSE)
table.in <- table.in[row.keep,]
tablechi <- tablechi[row.keep,]
} else {
yminmax <- c(min(table.in[,2], na.rm = TRUE) + ((max(cex.par)/10) * min(table.in[,2], na.rm = TRUE)), max(table.in[,2], na.rm = TRUE) + ((max(cex.par)/10) * max(table.in[,2], na.rm = TRUE)))
}
-if (typegraph == 0) {
+if (typegraph == 0 || typegraph == 2) {
open_file_graph(fileout, width = width, height = height, svg = do.svg)
parcex <- taillecar/10
par(cex = parcex)
if (over) {
- table.in <- table.in[order(cex.par, decreasing = TRUE),]
- classes <- classes[order(cex.par, decreasing = TRUE)]
- cex.par <- cex.par[order(cex.par, decreasing = TRUE)]
- table.in <- stopoverlap(table.in, cex.par=cex.par, xlim = xminmax, ylim = yminmax)
- classes <- classes[table.in[,4]]
- cex.par <- cex.par[table.in[,4]]
+ table.in <- table.in[order(cex.par, decreasing = TRUE),]
+ classes <- classes[order(cex.par, decreasing = TRUE)]
+ cex.par <- cex.par[order(cex.par, decreasing = TRUE)]
+ table.out <- stopoverlap(table.in, cex.par=cex.par, xlim = xminmax, ylim = yminmax)
+ table.in <- table.out$toplot
+ notplot <- table.out$notplot
+ if (! is.null(notplot)) {
+ write.csv2(notplot, file = paste(fileout,'_notplotted.csv', sep=''))
+ }
+ classes <- classes[table.in[,4]]
+ cex.par <- cex.par[table.in[,4]]
}
- make_afc_graph(table.in, classes, clnb, xlab, ylab, cex.txt = cex.par, xminmax = xminmax, yminmax = yminmax)
-
-} else {
- library(rgl)
- rn <- vire.nonascii(rownames(table.in))
- rain = rainbow(clnb)
- colors = rain[classes]
- #rn <- rownames(table.in)
- rgl.open()
- text3d(table.in[,1], table.in[,2], table.in[,3], rn, col = colors , cex = cex.par)
- rgl.bg(col = c('white', "#99bb99"), front = "lines", box=FALSE, sphere = TRUE)
- par3d(cex=0.7)
- #par3d(windowRect = c(100,100,600,600))
- rgl.lines(c(rx), c(0, 0), c(0, 0), col = "#000000")
- rgl.lines(c(0,0),c(ry),c(0,0),col = "#000000")
- rgl.lines(c(0,0),c(0,0),c(rz),col = "#000000")
- text3d(rx[2]+1,0,0, xlab)
- text3d(0,ry[2]+1,0, ylab)
- text3d(0,0,rz[2]+1, zlab)
-
- if (tchi) {
- maxchi <- norm.vec(maxchi, tchi.min/100, tchi.max/100)
- } else if (!is.null(cex.par)) {
- maxchi <- norm.vec(cex.par, txt.min/100, txt.max/100)
+ if (typegraph == 0) {
+ make_afc_graph(table.in, classes, clnb, xlab, ylab, cex.txt = cex.par, xminmax = xminmax, yminmax = yminmax)
} else {
- maxchi <- 0.1
- }
-
-
- for (i in 1:clnb) {
- text3d(rx[2],(ry[2]+(0.2*i)),0,paste('classe',i),col=rain[i])
- }
- #if (tchi) {
- # rgl.spheres(table.in, col = colors, radius = maxchi, alpha = alpha)
- #}
-
- if (dofilm) {
- require(tcltk)
- ReturnVal <- tkmessageBox(title="RGL 3 D",message="Cliquez pour commencer le film",icon="info",type="ok")
-
- movie3d(spin3d(axis=c(0,1,0),rpm=6), movie = 'film', frames = "tmpfilm", duration=10, clean=TRUE, top = TRUE, dir = dirout)
- ReturnVal <- tkmessageBox(title="RGL 3 D",message="Fini !",icon="info",type="ok")
+ dev.off()
+ require(rgexf)
+ nodes.attr <- make.afc.attributes(rownames(table.in), afc_table, afctable, clnb)
+ if (qui != 3) {
+ tokeep <- rownames(chistabletot) %%in%% rownames(table.in)
+ chis <- chistabletot[tokeep,]
+ chis<-chis[rownames(table.in),]
+ nodes.attr$chiclasse <- chis
+ } else {
+ chis <- NULL
+ }
+ afctogexf(fileout, table.in, classes, clnb, cex.par, nodes.attr = nodes.attr)
}
- require(tcltk)
- ReturnVal <- tkmessageBox(title="RGL 3 D",message="Cliquez pour fermer",icon="info",type="ok")
- rgl.close()
+} else {
+ if (typegraph == 4) {
+ library(igraph)
+ rain = rainbow(clnb)
+ col = rain[classes]
+ g <- make_empty_graph()
+ vertex <- rownames(table.in)
+ g <- add.vertices(g, length(vertex), attr=list(weight=cex.par, names=vertex, color=col))
+ minx <- min(table.in[,1])
+ maxx <- max(table.in[,1])
+ miny <- min(table.in[,2])
+ maxy <- max(table.in[,2])
+ minz <- min(table.in[,3])
+ maxz <- max(table.in[,3])
+ table.in <- rbind(table.in, c(minx, 0, 0))
+ rminx <- nrow(table.in)
+ table.in <- rbind(table.in, c(maxx, 0, 0))
+ rmaxx <- nrow(table.in)
+ table.in <- rbind(table.in, c(0, miny, 0))
+ rminy <- nrow(table.in)
+ table.in <- rbind(table.in, c(0, maxy, 0))
+ rmaxy <- nrow(table.in)
+ table.in <- rbind(table.in, c(0, 0, minz))
+ rminz <- nrow(table.in)
+ table.in <- rbind(table.in, c(0, 0, maxz))
+ rmaxz <- nrow(table.in)
+ g <- add.vertices(g, 6, attr=list(weight=c(0.1,0.1,0.1,0.1,0.1,0.1), names=c(rminx,rmaxx,rminy,rmaxy,rminz,rmaxz), color=c('white','white','white','white','white','white')))
+ g <- g + edge(rminx, rmaxx, weight=0.1) + edge(rminy, rmaxy, weight=0.1) + edge(rminz, rmaxz, weight=0.1)
+ table.in <- layout.norm(table.in, -5,5,-5,5,-5,5)
+ graph.to.file2(g, table.in, nodesfile=nodesfile, edgesfile=edgesfile)
+ } else {
+ library(rgl)
+ rn <- vire.nonascii(rownames(table.in))
+ rain = rainbow(clnb)
+ colors = rain[classes]
+ #rn <- rownames(table.in)
+ #rgl.open()
+ bg3d('white')
+ #rgl.bg(col = c('white', "#99bb99"), front = "lines", box=FALSE, sphere = TRUE)
+
+ par3d('userMatrix' = matrix(c(1,0,0,0, 0,1,0,0,0,0,1,0,0,0,0,1), ncol=4, nrow = 4))
+ #par3d(cex=0.7)
+ #par3d(windowRect = c(100,100,600,600))
+ rgl.lines(c(rx), c(0, 0), c(0, 0), col = "#000000")
+ rgl.lines(c(0,0),c(ry),c(0,0),col = "#000000")
+ rgl.lines(c(0,0),c(0,0),c(rz),col = "#000000")
+ text3d(rx[2]+1,0,0, xlab)
+ text3d(0,ry[2]+1,0, ylab)
+ text3d(0,0,rz[2]+1, zlab)
+ splt <- split(seq_along(table.in[,1]), ceiling(seq_along(table.in[,1])/100))
+ #colsplt <- split(seq_along(colors), ceiling(seq_along(colors)/100))
+ #cexsplt <- split(seq_along(cex.par), ceiling(seq_along(cex.par)/100))
+ for (i in splt) {
+ rgl.texts(table.in[i,1], table.in[i,2], table.in[i,3], rn[i], col = colors[i] , cex = cex.par[i])
+ }
+ #rgl.texts(table.in[,1], table.in[,2], table.in[,3], rn, col = colors , cex = cex.par)
+
+ if (tchi) {
+ maxchi <- norm.vec(maxchi, tchi.min/100, tchi.max/100)
+ } else if (!is.null(cex.par)) {
+ maxchi <- norm.vec(cex.par, txt.min/100, txt.max/100)
+ } else {
+ maxchi <- 0.1
+ }
+
+
+ #for (i in 1:clnb) {
+ # text3d(rx[2],(ry[2]+(0.2*i)),0,paste('classe',i),col=rain[i])
+ #}
+ if (tchi) {
+ rgl.spheres(table.in, col = colors, radius = maxchi, alpha = alpha)
+ }
+ par3d(skipRedraw=FALSE)
+
+ if (dofilm) {
+ require(tcltk)
+ ReturnVal <- tkmessageBox(title="RGL 3 D",message="Cliquez pour commencer le film",icon="info",type="ok")
+
+ movie3d(spin3d(axis=c(0,1,0),rpm=6), movie = 'film', frames = "tmpfilm", duration=10, clean=TRUE, top = TRUE, dir = dirout)
+ ReturnVal <- tkmessageBox(title="RGL 3 D",message="Fini !",icon="info",type="ok")
+ }
+
+ if (typegraph == 1) {
+ require(tcltk)
+ ReturnVal <- tkmessageBox(title="RGL 3 D",message="Cliquez pour fermer",icon="info",type="ok")
+ } else {
+ writeWebGL(dir = fileout, width = width, height= height)
+ }
+ rgl.close()
+ }
}