• R/O
  • SSH

标签
No Tags

Frequently used words (click to add to your profile)

javac++androidlinuxc#windowsobjective-ccocoa誰得qtpythonphprubygameguibathyscaphec計画中(planning stage)翻訳omegatframeworktwitterdomtestvb.netdirectxゲームエンジンbtronarduinopreviewer

File Info

Rev. 932c039d9e58e0a14ec72293bf4015ac5c7bf39e
大小 83,434 字节
时间 2011-06-23 07:08:10
作者 lorenzo
Log Message

I am now generating new panel plots.

Content

rm(list=ls())
## library(igraph)
library(ggplot2)
library(tikzDevice)


vplayout <- function(x, y) viewport(layout.pos.row=x, layout.pos.col=y)


draw4eps <- function(pdfname, a,b,c,d,w,h) {

   postscript(pdfname, width=w,height=h,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )


    ## pdf(pdfname,width=w, height=h)
    grid.newpage()
    pushViewport(viewport(layout=grid.layout(2,2) ) )

        print(a, vp=vplayout(1,1))
        print(b, vp=vplayout(1,2))
        print(c, vp=vplayout(2,1))
        print(d, vp=vplayout(2,2))


    dev.off()
} 




N <- 42+36
nq <- 13


#size of panel plots
width <- 10
height <- 10


plotall <- 1


data <- read.csv("GRIGLIA2.csv", header=FALSE)

data$V1 <- levels(data$V1)[data$V1]

names(data) <- c("Domanda", "M","F","20-30","31-40","41-50","51-60",">60",
                 "<1","1-5",">5","moltissimo","molto","abbastanza","poco",
                 "per niente", "Italia", "Svizzera")


## new_D <- c()

## for (i in seq(length(data$Domanda))){

## temp <- unlist(strsplit(data$Domanda[i], split=""))[1]
  
## new_D <- c(new_D,temp)

## }

my_seq <- rep(seq(nq), N)

data$Domanda <- my_seq


temp <- seq(N*nq)



MF <- seq(N*nq)

swissita <- seq(N*nq)


experience <- seq(N*nq)
age <- seq(N*nq)


motivation_mf <- seq(N*6)
motivation_age<- seq(N*6)
motivation_nat <- seq(N*6)
motivation_exp <- seq(N*6)

temp_mot <- seq(N*6)

training_mf <- seq(N*6)
training_age<- seq(N*6)
training_nat <- seq(N*6)
training_exp <- seq(N*6)

temp_tra<- seq(N*6)




temp[] <- -10

for (i in seq(nq)){


sel <- which(data$Domanda==i)

data_restr <- data[sel, ]

for (j in seq(N)){

## print("i and j are")
## print(i)
## print(j)

## print("data_restr[j,12:16] is, ")
## print(data_restr[j,12:16])

sel2 <- which(data_restr[j,12:16]==1)

if (length(sel2)!=1){

print("length(sel2) is, ")

print(length(sel2))

print("(i-1)*N+j is")
print((i-1)*N+j)

print("data_restr[j,12:16] is, ")

print(data_restr[j,12:16])

}
temp[(i-1)*N+j] <- sel2

if (data_restr$M[j]==1){
MF[(i-1)*N+j] <- "Uomini"
}else{

MF[(i-1)*N+j] <- "Donne"

}




if (data_restr$Italia[j]==1){
swissita[(i-1)*N+j] <- "Italiani"
}else{

swissita[(i-1)*N+j] <- "Svizzeri"

}




if (data_restr$"<1"[j]==1){
experience[(i-1)*N+j] <- "<1 anno"
} else if (data_restr$"1-5"[j]==1){
experience[(i-1)*N+j] <- "1-5 anni"
} else {
experience[(i-1)*N+j] <- ">5 anni"
  
}

if (data_restr$"20-30"[j]==1){
age[(i-1)*N+j] <- "20-30"
} else if (data_restr$"31-40"[j]==1){
age[(i-1)*N+j] <- "31-40"
} else if (data_restr$"41-50"[j]==1){
age[(i-1)*N+j] <- "41-50"
} else if (data_restr$"51-60"[j]==1){
age[(i-1)*N+j] <- "51-60"
} else {
age[(i-1)*N+j] <- ">60"

}




}
}



for (i in seq(N)){

sel <- (i-1)*nq
  
motivation_mf[((i-1)*6+1):((i-1)*6+6)] <- MF[((i-1)*nq+1):((i-1)*nq+6)]
motivation_exp[((i-1)*6+1):((i-1)*6+6)] <- experience[((i-1)*nq+1):((i-1)*nq+6)]
motivation_nat[((i-1)*6+1):((i-1)*6+6)] <- swissita[((i-1)*nq+1):((i-1)*nq+6)]
motivation_age[((i-1)*6+1):((i-1)*6+6)] <- age[((i-1)*nq+1):((i-1)*nq+6)]

temp_mot[((i-1)*6+1):((i-1)*6+6)] <- temp[((i-1)*nq+1):((i-1)*nq+6)] 

training_mf[((i-1)*6+1):((i-1)*6+6)] <- MF[((i-1)*nq+7):((i-1)*nq+12)]
training_exp[((i-1)*6+1):((i-1)*6+6)] <- experience[((i-1)*nq+7):((i-1)*nq+12)]
training_nat[((i-1)*6+1):((i-1)*6+6)] <- swissita[((i-1)*nq+7):((i-1)*nq+12)]
training_age[((i-1)*6+1):((i-1)*6+6)] <- age[((i-1)*nq+7):((i-1)*nq+12)]

temp_tra[((i-1)*6+1):((i-1)*6+6)] <- temp[((i-1)*nq+7):((i-1)*nq+12)]

}


## data_MF <- as.data.frame(cbind(temp,MF))
## names(data_MF) <- c("voto", "sesso")




print("Done with the data")

if (plotall==1){

for (i in seq(nq)){

  plist <- list()
  plist2 <- list()

  
  ini <- (i-1)*N+1
  fin <- i*N
  ans <- temp[ini:fin] 

  MF_ans <- MF[ini:fin] 
  
  ans <- as.data.frame(ans)
  names(ans) <- c("val")

  exp_ans <- experience[ini:fin]
  age_ans <- age[ini:fin]

  swissita_ans <- swissita[ini:fin]
  
  ## sel_M <- which(MF_ans=="Uomini")

  ## x <- as.data.frame(ans[sel_M,])
  ## names(x) <- "val"

  ## x$set <- "Uomini"

  ## x$set <- factor(x$set)

  
  ## sel_F <- which(MF_ans=="Donne")

  ## y <- as.data.frame(ans[sel_F,])
  ## names(y) <- "val"

  ## y$set <- "Donne"

  ## y$set <- factor(y$set)

  ## ans <- ans## -0.5

  ## xyz <- rbind(x, y)

  xyz <- cbind(as.data.frame(ans), as.data.frame(MF_ans),
               as.data.frame(exp_ans), as.data.frame(age_ans), as.data.frame(swissita_ans))

  names(xyz) <- c("val", "set","experience", "age", "nationality")

  xyz$set <-  factor(xyz$set,c("Uomini", "Donne") )
  xyz$experience <-  factor(xyz$experience, c("<1 anno", "1-5 anni", ">5 anni"))

  xyz$nationality <-  factor(xyz$nationality, c("Italiani", "Svizzeri"))

  
  ## xyz$age <- factor(xyz$age, c(""))

  
  lbls <- levels(xyz$set)

  fnq <- paste("panel_question_",i,"_.eps", sep="")

  
  fn <- paste("aggregated_",i,"_.tex", sep="")
  fnps <- paste("aggregated_",i,"_.eps", sep="")

  fnmf <- paste("aggregated_MF_",i,"_.tex", sep="")
  fnpsmf <- paste("aggregated_MF_",i,"_.eps", sep="")

  fnpsexp <- paste("aggregated_experience_",i,"_.eps", sep="")
  fnpsage <- paste("aggregated_age_",i,"_.eps", sep="")

  fnpsnat <- paste("aggregated_nationality_",i,"_.eps", sep="")



  fnpsmfnorm <- paste("aggregated_MF_norm_",i,"_.eps", sep="")

  fnpsexpnorm <- paste("aggregated_experience_norm_",i,"_.eps", sep="")
  fnpsagenorm <- paste("aggregated_age_norm_",i,"_.eps", sep="")

  fnpsnatnorm <- paste("aggregated_nationality_norm_",i,"_.eps", sep="")


  fnpsmfnormd <- paste("aggregated_MF_norm_dodge_",i,"_.eps", sep="")

  fnpsexpnormd <- paste("aggregated_experience_norm_dodge_",i,"_.eps", sep="")
  fnpsagenormd <- paste("aggregated_age_norm_dodge_",i,"_.eps", sep="")

  fnpsnatnormd <- paste("aggregated_nationality_norm_dodge_",i,"_.eps", sep="")


  
  gpl2 <- ggplot(ans,aes(x=val))+
    geom_histogram(aes(y=..count..),binwidth=1,
                   colour="black", fill="blue",origin=0.5)+
    ## opts( panel.background=theme_rect(fill='white', size=1.5))+
    opts(panel.grid.minor = theme_blank())+
    ## opts(panel.grid.major = theme_blank())+
    scale_x_continuous(limits=c(0.5,5.5), breaks=seq(5))+
  ## scale_y_continuous(limits=c(0,.4), breaks=seq(0,0.4, by=0.1))+
scale_y_continuous(breaks=seq(0,300, by=2))+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
  xlab("Risposta")+
  ylab("Numero di risposte")+
  opts(axis.title.x = theme_text(size = 20))+
  opts(axis.title.y = theme_text(size = 20, angle=90,hjust=0.5))+
  opts(plot.title = theme_text(size = 25))+
  opts(axis.text.x = theme_text(size=18, colour="black",vjust=1))+
  opts(axis.text.y = theme_text(size=18, colour="black", hjust=1))





## tikz(fn, standAlone = TRUE, width=5,height=5)
## print(gpl2)
## dev.off() 

## tools::texi2dvi(fn,pdf=T)

postscript(fnps, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl2)
dev.off() 





  
##############################################################àà

gpl <- ggplot(xyz, aes(x=val,y=..count.., fill=set, group=set)) +
  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+

  
 

## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +
 scale_fill_manual("Sesso", breaks=lbls, labels=lbls, values=c("blue","red"))+

  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Numero di risposte")+

scale_y_continuous(breaks=seq(0,300, by=2))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


  
 


## tikz(fnmf, standAlone = TRUE, width=5,height=5)
## print(gpl)
## dev.off() 

## tools::texi2dvi(fnmf,pdf=T)

plist[[1]] <- gpl

postscript(fnpsmf, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 




if (i==13){

  gpl <- ggplot(xyz, aes(x=val,y=..count../sum(..count..)*100, fill=set, group=set)) +
  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+

  
 

## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +
 scale_fill_manual("Sesso", breaks=lbls, labels=lbls, values=c("blue","red"))+

  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Percentuale di risposte")+

scale_y_continuous(breaks=seq(0,100, by=5))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


  
 


## tikz(fnmf, standAlone = TRUE, width=5,height=5)
## print(gpl)
## dev.off() 

## tools::texi2dvi(fnmf,pdf=T)


postscript(fnpsmfnorm, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 

plist2[[1]] <- gpl

    gpl <- ggplot(xyz, aes(x=val,y=..count../sum(..count..)*100, fill=set, group=set)) +
  geom_histogram(alpha=1,
                 binwidth=1,  position="dodge",
                 origin=0.5 ##,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1,  position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+

  
 

## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +
 scale_fill_manual("Sesso", breaks=lbls, labels=lbls, values=c("blue","red"))+

  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Percentuale di risposte")+

scale_y_continuous(breaks=seq(0,100, by=5))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


  
 


## tikz(fnmf, standAlone = TRUE, width=5,height=5)
## print(gpl)
## dev.off() 

## tools::texi2dvi(fnmf,pdf=T)


postscript(fnpsmfnormd, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 




}


  
#############################################################à

lbls <- levels(xyz$experience)

  
gpl <- ggplot(xyz, aes(x=val,y=..count.., fill=experience, group=experience)) +
  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +
 scale_fill_manual("Anzianità", breaks=lbls,
                   labels=lbls, values=c("blue","red",
                                                   "darkgreen"))+

  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Numero di risposte")+

scale_y_continuous(breaks=seq(0,300, by=2))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsexp, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off()

plist[[2]] <- gpl


if (i==13){

gpl <- ggplot(xyz, aes(x=val,y=..count../sum(..count..)*100, fill=experience, group=experience)) +
  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +
 scale_fill_manual("Anzianità", breaks=lbls,
                   labels=lbls, values=c("blue","red",
                                                   "darkgreen"))+

  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Percentuale di risposte")+

scale_y_continuous(breaks=seq(0,300, by=5))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsexpnorm, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off()


plist2[[2]] <- gpl

gpl <- ggplot(xyz, aes(x=val,y=..count../sum(..count..)*100, fill=experience, group=experience)) +
  geom_histogram(alpha=1,
                 binwidth=1,  position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1,  position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +
 scale_fill_manual("Anzianità", breaks=lbls,
                   labels=lbls, values=c("blue","red",
                                                   "darkgreen"))+

  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Percentuale di risposte")+

scale_y_continuous(breaks=seq(0,300, by=5))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsexpnormd, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off()



}
  

#######################################################################

lbls <- levels(xyz$nationality)

  
gpl <- ggplot(xyz, aes(x=val,y=..count.., fill=nationality, group=nationality)) +
  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +
 scale_fill_manual("Nazionalità", breaks=lbls,
                   labels=lbls, values=c("blue","red",
                                                   "darkgreen"))+

  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Numero di risposte")+

scale_y_continuous(breaks=seq(0,300, by=2))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsnat, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 

plist[[3]] <- gpl


  
if (i==13){


gpl <- ggplot(xyz, aes(x=val,y=..count../sum(..count..)*100, fill=nationality, group=nationality)) +
  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +
 scale_fill_manual("Nazionalità", breaks=lbls,
                   labels=lbls, values=c("blue","red",
                                                   "darkgreen"))+

  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Percentuale di risposte")+

scale_y_continuous(breaks=seq(0,300, by=5))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsnatnorm, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 


plist2[[3]] <- gpl


gpl <- ggplot(xyz, aes(x=val,y=..count../sum(..count..)*100, fill=nationality, group=nationality)) +
  geom_histogram(alpha=1,
                 binwidth=1,  position="dodge",
                 origin=0.5 ##,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1,  position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +
 scale_fill_manual("Nazionalità", breaks=lbls,
                   labels=lbls, values=c("blue","red",
                                                   "darkgreen"))+

  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Percentuale di risposte")+

scale_y_continuous(breaks=seq(0,300, by=5))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsnatnormd, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 



}

  
##################################################################
lbls <- levels(xyz$age)

  
gpl <- ggplot(xyz, aes(x=val,y=..count.., fill=age, group=age)) +
  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +
 scale_fill_manual("Età", breaks=lbls, labels=lbls,
                   values=c("blue","red","darkgreen", "brown", "purple"))+

  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Numero di risposte")+

scale_y_continuous(breaks=seq(0,300, by=2))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsage, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 

plist[[4]] <- gpl



draw4eps(fnq, plist[[1]],plist[[2]],plist[[3]],plist[[4]], width, height) 


  
  if (i==13){


gpl <- ggplot(xyz, aes(x=val,y=..count../sum(..count..)*100, fill=age, group=age)) +
  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +
 scale_fill_manual("Età", breaks=lbls, labels=lbls,
                   values=c("blue","red","darkgreen", "brown", "purple"))+

  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Percentuale di risposte")+

scale_y_continuous(breaks=seq(0,300, by=5))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsagenorm, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 


plist2[[4]] <- gpl

draw4eps("cambiamento_panel.eps", plist2[[1]],plist2[[2]],plist2[[3]],plist2[[4]], width, height) 


gpl <- ggplot(xyz, aes(x=val,y=..count../sum(..count..)*100, fill=age, group=age)) +
  geom_histogram(alpha=1,
                 binwidth=1,  position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1,  position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +
 scale_fill_manual("Età", breaks=lbls, labels=lbls,
                   values=c("blue","red","darkgreen", "brown", "purple"))+

  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Percentuale di risposte")+

scale_y_continuous(breaks=seq(0,300, by=5))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsagenormd, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 



  }
  
  
}

}

########################################################################################
#aggregation by macroarea

###############macroarea about motivation

  plist <- list()


#for age

  xyz <- cbind(as.data.frame(temp_mot), as.data.frame(motivation_mf),
               as.data.frame(motivation_exp), as.data.frame(motivation_age), as.data.frame(motivation_nat))

  names(xyz) <- c("val", "set","experience", "age", "nationality")

  xyz$set <-  factor(xyz$set,c("Uomini", "Donne") )
  xyz$experience <-  factor(xyz$experience, c("<1 anno", "1-5 anni", ">5 anni"))

  xyz$nationality <-  factor(xyz$nationality, c("Italiani", "Svizzeri"))

  
  ## xyz$age <- factor(xyz$age, c(""))

  
  lbls <- levels(xyz$age)


fnpsage <- paste("motivation_age.eps")


gpl <- ggplot(xyz, aes(x=val,y=..count.., fill=age, group=age)) +
  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +
 scale_fill_manual("Età", breaks=lbls, labels=lbls,
                   values=c("blue","red","darkgreen", "brown", "purple"))+

  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Numero di risposte")+

scale_y_continuous(breaks=seq(0,400, by=10))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsage, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 



fnpsage <- paste("motivation_age_norm.eps")


gpl <- ggplot(xyz, aes(x=val,y=..count../sum(..count..)*100, fill=age, group=age)) +
  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +
 scale_fill_manual("Età", breaks=lbls, labels=lbls,
                   values=c("blue","red","darkgreen", "brown", "purple"))+

  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Percentuale di risposte")+

## scale_y_continuous(breaks=seq(0,400, by=10))+

    scale_y_continuous(breaks=seq(0,40, by=5), limits=c(0,42))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsage, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 

plist[[1]] <- gpl



fnpsage <- paste("motivation_age_norm_dodge.eps")


gpl <- ggplot(xyz, aes(x=val,y=..count../sum(..count..)*100, fill=age, group=age)) +
  geom_histogram(alpha=1,
                 binwidth=1,  position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +
 scale_fill_manual("Età", breaks=lbls, labels=lbls,
                   values=c("blue","red","darkgreen", "brown", "purple"))+

  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Percentuale di risposte")+

## scale_y_continuous(breaks=seq(0,400, by=10))+

    scale_y_continuous(breaks=seq(0,40, by=5))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsage, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 






########## for experience



  
  lbls <- levels(xyz$experience)


fnpsage <- paste("motivation_experience.eps")


gpl <- ggplot(xyz, aes(x=val,y=..count.., fill=experience, group=experience)) +
  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +


 scale_fill_manual("Anzianità", breaks=lbls,
                   labels=lbls, values=c("blue","red",
                                                   "darkgreen"))+

  
  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Numero di risposte")+

scale_y_continuous(breaks=seq(0,400, by=10))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsage, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 




fnpsage <- paste("motivation_experience_norm.eps")


gpl <- ggplot(xyz, aes(x=val,y=..count../sum(..count..)*100, fill=experience, group=experience)) +
  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +


 scale_fill_manual("Anzianità", breaks=lbls,
                   labels=lbls, values=c("blue","red",
                                                   "darkgreen"))+

  
  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Percentuale di risposte")+

## scale_y_continuous(breaks=seq(0,400, by=10))+
    scale_y_continuous(breaks=seq(0,40, by=5), limits=c(0,42))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsage, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 


plist[[2]] <- gpl

fnpsage <- paste("motivation_experience_norm_dodge.eps")


gpl <- ggplot(xyz, aes(x=val,y=..count../sum(..count..)*100, fill=experience, group=experience)) +
  geom_histogram(alpha=1,
                 binwidth=1,  position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1,  position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +


 scale_fill_manual("Anzianità", breaks=lbls,
                   labels=lbls, values=c("blue","red",
                                                   "darkgreen"))+

  
  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Percentuale di risposte")+

## scale_y_continuous(breaks=seq(0,400, by=10))+
    scale_y_continuous(breaks=seq(0,40, by=5))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsage, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 






########males/females


  lbls <- levels(xyz$set)


fnpsage <- paste("motivation_male-female.eps")


gpl <- ggplot(xyz, aes(x=val,y=..count.., fill=set, group=set)) +
  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +


 scale_fill_manual("Sesso", breaks=lbls,
                   labels=lbls, values=c("blue","red",
                                                   "darkgreen"))+

  
  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Numero di risposte")+

scale_y_continuous(breaks=seq(0,400, by=10))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsage, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 




fnpsage <- paste("motivation_male-female_norm.eps")


gpl <- ggplot(xyz, aes(x=val,y=..count../sum(..count..)*100, fill=set, group=set)) +
  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +


 scale_fill_manual("Sesso", breaks=lbls,
                   labels=lbls, values=c("blue","red",
                                                   "darkgreen"))+

  
  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Percentuale di risposte")+

## scale_y_continuous(breaks=seq(0,400, by=10))+
    scale_y_continuous(breaks=seq(0,40, by=5), limits=c(0,42))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsage, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 



plist[[3]] <- gpl

fnpsage <- paste("motivation_male-female_norm_dodge.eps")


gpl <- ggplot(xyz, aes(x=val,y=..count../sum(..count..)*100, fill=set, group=set)) +
  geom_histogram(alpha=1,
                 binwidth=1,  position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1,  position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +


 scale_fill_manual("Sesso", breaks=lbls,
                   labels=lbls, values=c("blue","red",
                                                   "darkgreen"))+

  
  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Percentuale di risposte")+

## scale_y_continuous(breaks=seq(0,400, by=10))+
    scale_y_continuous(breaks=seq(0,40, by=5))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsage, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 






#############nationality


  lbls <- levels(xyz$nationality)


fnpsage <- paste("motivation_nationality.eps")


gpl <- ggplot(xyz, aes(x=val,y=..count.., fill=nationality, group=nationality)) +
  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +


 scale_fill_manual("Nazionalità", breaks=lbls,
                   labels=lbls, values=c("blue","red",
                                                   "darkgreen"))+

  
  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Numero di risposte")+

scale_y_continuous(breaks=seq(0,400, by=10))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsage, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 


fnpsage <- paste("motivation_nationality_norm.eps")


gpl <- ggplot(xyz, aes(x=val,y=..count../sum(..count..)*100, fill=nationality, group=nationality)) +
  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+

  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +


 scale_fill_manual("Nazionalità", breaks=lbls,
                   labels=lbls, values=c("blue","red",
                                                   "darkgreen"))+

  
  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Percentuale di risposte")+

## scale_y_continuous(breaks=seq(0,400, by=10))+
    scale_y_continuous(breaks=seq(0,40, by=5), limits=c(0,42))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsage, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 




plist[[4]] <- gpl

draw4eps("motivation_panel.eps", plist[[1]],plist[[2]],plist[[3]],plist[[4]], width, height) 


fnpsage <- paste("motivation_nationality_norm_dodge.eps")


gpl <- ggplot(xyz, aes(x=val,y=..count../sum(..count..)*100, fill=nationality, group=nationality)) +
  geom_histogram(alpha=1,
                 binwidth=1,  position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1,  position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+

  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +


 scale_fill_manual("Nazionalità", breaks=lbls,
                   labels=lbls, values=c("blue","red",
                                                   "darkgreen"))+

  
  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Percentuale di risposte")+

## scale_y_continuous(breaks=seq(0,400, by=10))+
    scale_y_continuous(breaks=seq(0,40, by=5))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsage, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 





#########################################################################################
#########################################################################################

# macroarea about training



#for age

  xyz <- cbind(as.data.frame(temp_tra), as.data.frame(training_mf),
               as.data.frame(training_exp), as.data.frame(training_age), as.data.frame(training_nat))

  names(xyz) <- c("val", "set","experience", "age", "nationality")

  xyz$set <-  factor(xyz$set,c("Uomini", "Donne") )
  xyz$experience <-  factor(xyz$experience, c("<1 anno", "1-5 anni", ">5 anni"))

  xyz$nationality <-  factor(xyz$nationality, c("Italiani", "Svizzeri"))

  
  ## xyz$age <- factor(xyz$age, c(""))

  
  lbls <- levels(xyz$age)


fnpsage <- paste("training_age.eps")


gpl <- ggplot(xyz, aes(x=val,y=..count.., fill=age, group=age)) +
  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +
 scale_fill_manual("Età", breaks=lbls, labels=lbls,
                   values=c("blue","red","darkgreen", "brown", "purple"))+

  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Numero di risposte")+

scale_y_continuous(breaks=seq(0,400, by=10))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsage, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 




fnpsage <- paste("training_age_norm.eps")


gpl <- ggplot(xyz, aes(x=val,y=..count../sum(..count..)*100, fill=age, group=age)) +
  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+

  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +
 scale_fill_manual("Età", breaks=lbls, labels=lbls,
                   values=c("blue","red","darkgreen", "brown", "purple"))+

  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Percentuale di risposte")+

scale_y_continuous(breaks=seq(0,40, by=5), limits=c(0,40))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsage, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 


plist[[1]] <- gpl



fnpsage <- paste("training_age_norm_dodge.eps")


gpl <- ggplot(xyz, aes(x=val,y=..count../sum(..count..)*100, fill=age, group=age)) +
  geom_histogram(alpha=1,
                 binwidth=1,  position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+

  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +
 scale_fill_manual("Età", breaks=lbls, labels=lbls,
                   values=c("blue","red","darkgreen", "brown", "purple"))+

  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Percentuale di risposte")+

scale_y_continuous(breaks=seq(0,40, by=5))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsage, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 











########## for experience



  
  lbls <- levels(xyz$experience)


fnpsage <- paste("training_experience.eps")


gpl <- ggplot(xyz, aes(x=val,y=..count.., fill=experience, group=experience)) +
  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +


 scale_fill_manual("Anzianità", breaks=lbls,
                   labels=lbls, values=c("blue","red",
                                                   "darkgreen"))+

  
  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Numero di risposte")+

scale_y_continuous(breaks=seq(0,400, by=10))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsage, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 





fnpsage <- paste("training_experience_norm.eps")


gpl <- ggplot(xyz, aes(x=val,y=..count../sum(..count..)*100, fill=experience, group=experience)) +
  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +


 scale_fill_manual("Anzianità", breaks=lbls,
                   labels=lbls, values=c("blue","red",
                                                   "darkgreen"))+

  scale_y_continuous(breaks=seq(0,40, by=5), limits=c(0,40))+

  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Percentuale di risposte")+

## scale_y_continuous(breaks=seq(0,400, by=10))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsage, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 


plist[[2]] <- gpl 



fnpsage <- paste("training_experience_norm_dodge.eps")


gpl <- ggplot(xyz, aes(x=val,y=..count../sum(..count..)*100, fill=experience, group=experience)) +
  geom_histogram(alpha=1,
                 binwidth=1, position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +


 scale_fill_manual("Anzianità", breaks=lbls,
                   labels=lbls, values=c("blue","red",
                                                   "darkgreen"))+

  scale_y_continuous(breaks=seq(0,40, by=5))+

  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Percentuale di risposte")+

## scale_y_continuous(breaks=seq(0,400, by=10))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsage, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 






########males/females


  lbls <- levels(xyz$set)


fnpsage <- paste("training_male-female.eps")


gpl <- ggplot(xyz, aes(x=val,y=..count.., fill=set, group=set)) +
  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +


 scale_fill_manual("Sesso", breaks=lbls,
                   labels=lbls, values=c("blue","red",
                                                   "darkgreen"))+

  
  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Numero di risposte")+

scale_y_continuous(breaks=seq(0,400, by=10))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsage, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 




fnpsage <- paste("training_male-female_norm.eps")


gpl <- ggplot(xyz, aes(x=val,y=..count../sum(..count..)*100, fill=set, group=set)) +
  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +


 scale_fill_manual("Sesso", breaks=lbls,
                   labels=lbls, values=c("blue","red",
                                                   "darkgreen"))+

    scale_y_continuous(breaks=seq(0,40, by=5), limits=c(0,40))+

  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Percentuale di risposte")+

## scale_y_continuous(breaks=seq(0,400, by=10))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsage, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 



plist[[3]] <- gpl


fnpsage <- paste("training_male-female_norm_dodge.eps")


gpl <- ggplot(xyz, aes(x=val,y=..count../sum(..count..)*100, fill=set, group=set)) +
  geom_histogram(alpha=1,
                 binwidth=1,  position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1,  position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +


 scale_fill_manual("Sesso", breaks=lbls,
                   labels=lbls, values=c("blue","red",
                                                   "darkgreen"))+

    scale_y_continuous(breaks=seq(0,40, by=5))+

  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Percentuale di risposte")+

## scale_y_continuous(breaks=seq(0,400, by=10))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsage, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 






#############nationality


  lbls <- levels(xyz$nationality)


fnpsage <- paste("training_nationality.eps")


gpl <- ggplot(xyz, aes(x=val,y=..count.., fill=nationality, group=nationality)) +
  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +


 scale_fill_manual("Nazionalità", breaks=lbls,
                   labels=lbls, values=c("blue","red",
                                                   "darkgreen"))+

  
  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Numero di risposte")+

scale_y_continuous(breaks=seq(0,400, by=10))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsage, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 




fnpsage <- paste("training_nationality_norm.eps")


gpl <- ggplot(xyz, aes(x=val,y=..count../sum(..count..)*100, fill=nationality, group=nationality)) +
  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1, ## position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +


 scale_fill_manual("Nazionalità", breaks=lbls,
                   labels=lbls, values=c("blue","red",
                                                   "darkgreen"))+

    scale_y_continuous(breaks=seq(0,40, by=5), limits=c(0,40))+

  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Percentuale di risposte")+

## scale_y_continuous(breaks=seq(0,400, by=10))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsage, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 


plist[[4]] <- gpl


draw4eps("training_panel.eps", plist[[1]],plist[[2]],plist[[3]],plist[[4]], width, height) 


fnpsage <- paste("training_nationality_norm_dodge.eps")


gpl <- ggplot(xyz, aes(x=val,y=..count../sum(..count..)*100, fill=nationality, group=nationality)) +
  geom_histogram(alpha=1,
                 binwidth=1, position="dodge",
                 origin=0.5## ,colour="black"
                 )+

  geom_histogram(alpha=1,
                 binwidth=1,  position="dodge",
                 origin=0.5 ,colour="black",
                 legend=FALSE
                 )+

              ## opts( panel.background=theme_rect(fill='white', size=1.5))+
scale_x_continuous()+
scale_x_continuous()+
  opts(legend.text = theme_text(size = 14, vjust=0.4))+
  opts(legend.title = theme_text(size = 18,hjust=0))+
  opts(legend.position = c(0.35, 0.75), legend.background=theme_rect(col=0),
       legend.key = theme_rect(colour = NA))+


## scale_colour_manual("", breaks=lbls, labels=lbls, values="black") +


 scale_fill_manual("Nazionalità", breaks=lbls,
                   labels=lbls, values=c("blue","red",
                                                   "darkgreen"))+

    scale_y_continuous(breaks=seq(0,40, by=5))+

  
  opts(panel.grid.minor = theme_blank())+
    opts(axis.ticks = theme_segment(colour = "black", size=1),
         axis.ticks.length = unit(0.15, "cm"))+
## opts(panel.grid.major = theme_blank())+
## xlab("$\\Pi^{max}$")+
## ylab("$P(\\Pi^{max})$")+
  xlab("Risposta")+
  ylab("Percentuale di risposte")+

## scale_y_continuous(breaks=seq(0,400, by=10))+

  opts(axis.title.x = theme_text(size = 20))+
opts(axis.title.y = theme_text(size = 20, angle=90))+
#    opts(title = expression(paste("Dublin")))+
opts(axis.text.x = theme_text(size=15, colour="black",vjust=1))+
    opts(plot.title = theme_text(size = 25))+
     ## opts(plot.title = theme_blank())+
 ## opts(title = expression(paste("SG+elsewhere")))+
 
  
opts(axis.text.y = theme_text(size=15, colour="black", hjust=1))


postscript(fnpsage, width=5,height=5,horizontal = FALSE,
           onefile = FALSE, paper = "special"## ,
           ## family = "ComputerModern", encoding = "TeXtext.enc"
           )
print(gpl)
dev.off() 




print("So far so good")