#par(family="Arial") tiff("gumbels.tiff",width=8.3,height=8.3,units="cm",res=300,pointsize=8) library(evd) #Score <- read.table("PDscore/gapless/data/PDNoSameFold",colClasses=c("numeric","numeric"),head=F) #Score <- Score[,2] #obs<- table(cut(Score,breaks=breaks)) breaks<-seq(0,1,0.01) breaks1000<-seq(0,1.01,0.001) X=read.table("PDObs100") obsPD<-X[,2] loc = 0.05773573 scale =0.02719071 teoricasPD <- pgumbel(breaks1000[11:1010],loc,scale )-pgumbel(breaks1000[1:1000],loc,scale ) # sol<-chisq.test(obs,p=teoricas,rescale.p=T) # print(sol) X=read.table("TMObs100") obsTM<-X[,2] loc = 0.16203446 scale =0.02108873 teoricasTM <- pgumbel(breaks1000[11:1010],loc,scale )-pgumbel(breaks1000[1:1000],loc,scale ) breaks=breaks[1:100] breaks1000=breaks1000[5:1004] #par(mar =c(3.8,3.8,0.5,0.5)) plot(breaks,obsTM/sum(obsTM),type="s",ylab="Relative Frequency",xlab="Score") lines(breaks1000,teoricasTM,col="red") lines(breaks1000,teoricasPD,col="red") lines(breaks,obsPD/sum(obsPD),type="s") dev.off() #index=[1:100] #esp=teoricas[index]*sum(obs[index]) #diff=(esp-obs[index])^2 #sdobs=sd(obs[index]) #R2=1-sum(diff)/(sdobs*sdobs*length(obs[index])) #print(R2)