Hola a todos, Estoy intentando utilizar este script para hacer un plot con valores x ,y,z; representando los valores TP(y) y FP(x) y en funcion de la velocidad que seria el factor alpha. Y me da este en el ultimo punto de hacer el plot, alguien sabe que significa? library(ROCR) data(ROCR.simple) pred <- prediction( ROCR.simple$predictions, ROCR.simple$labels ) perf <- performance( pred, "tpr", "fpr" ) filename<-"rocdata csv.csv" Anthias<- read.table(filename, sep=",",header=TRUE) windows(width = 9.5, height = 10.50) par(mfrow = c(1,1)) cc<-Anthias[1] ee<-Anthias[2] aa<-Anthias[3] perf en x.name<-"FP" perf en y.name<-"TP" perf en alpha.name<-"Speed" bb<-data.frame(aa) perf en x.values<-bb dd<-data.frame(cc) perf en y.values<-dd ff<-data.frame(ee) perf en alpha.values<-ff plot(perf, main="Teste", colorize=TRUE, coloraxis.at=c(0,2,4,6,8,10), coloraxis.cex.axis=0.8, colorize.palette=(rainbow(256,start=0, end=0.7)), lwd=7, colorkey.relwidth=0.4, yaxis.las=1) Si alguien puede darme alguna idea, agradeceria.Muchas gracias¡ ------------ próxima parte ------------ Se ha borrado un adjunto en formato HTML... URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20140627/8faabbad/attachment-0001.html> ------------ próxima parte ------------ A non-text attachment was scrubbed... Name: rocdata csv.csv Type: text/csv Size: 20763 bytes Desc: no disponible URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20140627/8faabbad/attachment-0001.bin>
Hola, A la hora de leer los datos, tienes dos errores: el separador de campos es ";" y los decimales son con "," en vez de ".". Modificando los parámetros de read.table en la forma que te indico, desaparece el problema con el gráfico. Aunque el gráfico que obtienes ya te adelanto que no tiene mucho sentido... #------------------- Anthias<- read.table(filename, sep="*;*",header=TRUE, *dec=","*) #------------------- Saludos, Carlos Ortega www.qualityexcellence.es El 27 de junio de 2014, 21:09, Marta valdes lopez <martavaldes85@gmail.com> escribió:> Hola a todos, > > Estoy intentando utilizar este script para hacer un plot con valores x > ,y,z; representando los valores TP(y) y FP(x) y en funcion de la velocidad > que seria el factor alpha. > Y me da este en el ultimo punto de hacer el plot, alguien sabe que > significa? > > library(ROCR) > data(ROCR.simple) > pred <- prediction( ROCR.simple$predictions, ROCR.simple$labels ) > perf <- performance( pred, "tpr", "fpr" ) > filename<-"rocdata csv.csv" > > Anthias<- read.table(filename, sep=",",header=TRUE) > > windows(width = 9.5, height = 10.50) > par(mfrow = c(1,1)) > > cc<-Anthias[1] > ee<-Anthias[2] > aa<-Anthias[3] > > perf@x.name<-"FP" > perf@y.name<-"TP" > perf@alpha.name<-"Speed" > > bb<-data.frame(aa) > perf@x.values<-bb > > dd<-data.frame(cc) > perf@y.values<-dd > > ff<-data.frame(ee) > perf@alpha.values<-ff > > > plot(perf, main="Teste", colorize=TRUE, coloraxis.at=c(0,2,4,6,8,10), > coloraxis.cex.axis=0.8, colorize.palette=(rainbow(256,start=0, end=0.7)), > lwd=7, colorkey.relwidth=0.4, yaxis.las=1) > > > Si alguien puede darme alguna idea, agradeceria.Muchas gracias¡ > > _______________________________________________ > R-help-es mailing list > R-help-es@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es > >-- Saludos, Carlos Ortega www.qualityexcellence.es [[alternative HTML version deleted]]
Hola, Pues muchas gracias Carlos, no me acabo de aclarar con lo de las comas o puntos y comas para leer los archivos. Ya probe y funciono pero si queda bastante raro.. tambien probe con el programa lattice puede estar bien tambien. Muchas gracias, un saludo El 27 de junio de 2014, 22:58, Carlos Ortega <cof@qualityexcellence.es> escribió:> Hola, > > No sé si lo siguiente te resultará útil, pero sí que podrás ver de forma > fácil las diferencias de lo que has obtenido ("x" e "y") para cada nivel de > alpha. Para poder utilizar el siguiente código tienes que descargarte el > paquete "lattice". > > #----------------------------- > library(lattice) > names(Anthias) <- c('y', 'alpha','x') > alpha.tmp <- unique(Anthias$alpha) > alpha.fac <- factor(alpha.tmp, levels=alpha.tmp, ordered=T) > xyplot( > y ~ x | alpha.fac > ,data=Anthias > ,layout=c(5,2) > ,as.table=T > ) > #----------------------------- > > Saludos, > Carlos Ortega > www.qualityexcellence.es > > > El 28 de junio de 2014, 0:36, Carlos Ortega <cof@qualityexcellence.es> > escribió: > > Hola, >> >> A la hora de leer los datos, tienes dos errores: el separador de campos >> es ";" y los decimales son con "," en vez de ".". >> Modificando los parámetros de read.table en la forma que te indico, >> desaparece el problema con el gráfico. Aunque el gráfico que obtienes ya te >> adelanto que no tiene mucho sentido... >> >> #------------------- >> Anthias<- read.table(filename, sep="*;*",header=TRUE, *dec=","*) >> #------------------- >> >> Saludos, >> Carlos Ortega >> www.qualityexcellence.es >> >> >> El 27 de junio de 2014, 21:09, Marta valdes lopez < >> martavaldes85@gmail.com> escribió: >> >>> Hola a todos, >>> >>> Estoy intentando utilizar este script para hacer un plot con valores x >>> ,y,z; representando los valores TP(y) y FP(x) y en funcion de la velocidad >>> que seria el factor alpha. >>> Y me da este en el ultimo punto de hacer el plot, alguien sabe que >>> significa? >>> >>> library(ROCR) >>> data(ROCR.simple) >>> pred <- prediction( ROCR.simple$predictions, ROCR.simple$labels ) >>> perf <- performance( pred, "tpr", "fpr" ) >>> filename<-"rocdata csv.csv" >>> >>> Anthias<- read.table(filename, sep=",",header=TRUE) >>> >>> windows(width = 9.5, height = 10.50) >>> par(mfrow = c(1,1)) >>> >>> cc<-Anthias[1] >>> ee<-Anthias[2] >>> aa<-Anthias[3] >>> >>> perf@x.name<-"FP" >>> perf@y.name<-"TP" >>> perf@alpha.name<-"Speed" >>> >>> bb<-data.frame(aa) >>> perf@x.values<-bb >>> >>> dd<-data.frame(cc) >>> perf@y.values<-dd >>> >>> ff<-data.frame(ee) >>> perf@alpha.values<-ff >>> >>> >>> plot(perf, main="Teste", colorize=TRUE, coloraxis.at=c(0,2,4,6,8,10), >>> coloraxis.cex.axis=0.8, colorize.palette=(rainbow(256,start=0, end=0.7)), >>> lwd=7, colorkey.relwidth=0.4, yaxis.las=1) >>> >>> >>> Si alguien puede darme alguna idea, agradeceria.Muchas gracias¡ >>> >>> _______________________________________________ >>> R-help-es mailing list >>> R-help-es@r-project.org >>> https://stat.ethz.ch/mailman/listinfo/r-help-es >>> >>> >> >> >> -- >> Saludos, >> Carlos Ortega >> www.qualityexcellence.es >> > > > > -- > Saludos, > Carlos Ortega > www.qualityexcellence.es >[[alternative HTML version deleted]]