similar to: script works in Rgui, but failes in Rscript (coords, package 'pROC')

Hi I try to plot a ROC curve, using the ROC library from bioonductor (using the example in the help). Everything runs nicely except, that finally I get the error: "Error in plot.window(xlim, ylim, log, asp, ...) : need finite xlim values" --- > set.seed(123) > state <- c(0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 1, 1, + 1, 1, 1) > markers <- c(1, 2, 1, 1, 2,

Dear useRs: I use pROC package to compute the bootstrap C.I. of AUC. The command was as follows: roc1<-roc(all$D,all$pre,ci=TRUE,boot.n=200) However, the result was: Area under the curve: 0.5903 95% CI: 0.479-0.7016 (DeLong) Why the C.I. was computed by the Delong Method? Yao Zhu Department of Urology Fudan University Shanghai Cancer Center Shanghai, China [[alternative HTML version

I am currently trying to use R to construct a regression model to explain output based on temperature. I have combined my output and temp data into a notepad file. there is no problem with loading the data into R. > data.df output temp 1 850 17 2 849 17 3 905 17 4 925 17 5 1043 19 6 1104 20 7 1097 18 8 979 19 9 926 18 10 1133 18 ~~ ~~ 240 1124

Is there a way to teach xy.coords, when given log="xy", or just "x" or "y" separately, to do a decimal log10 instead of the natural log? Cheers, Alexy

I'm now trying to read.zoo in a rather long txt file with two columns: date/time and value in kW e.g. 432.2189 The read.zoo function finally ran without errors but not sure it is correct. I took the header off, and put in commas and added a " at the beginning and " at the end. z=read.zoo("Kevin-0-comma-ITPower.txt","%m/%d/%Y

Is there a function to convert from plot coordinates to device coordinates? It should be easy enough to work out from various par() parameters - namely usr, plt, and fig. 'usr' tells us where an (x,y) point is in the plot, then 'plt' tells us where it is in the figure, then 'fig' tells us where it is in the ((0,1),(0,1)) device space. Then there's the slight

In ?xy.coords it says: If 'y' is missing and 'x' is a formula: of the form 'yvar ~ xvar'. 'xvar' and 'yvar' are used as x and y variables. list: containing components 'x' and 'y', these are used to define plotting coordinates. time series: the x values are taken to be 'time(x)' and the y

Dear R users, I would like to draw map and import it in maptools/spatstat packages. The 'raw data' I have come from a web page (<map>...</map>) and are basically a list of coordinates of a polygon. I would like to know how to import them in R; I checked the maptools packages, but all the examples use existing .dbf files. I just have a (serie of) text file(s) looking like

Hi, I'm trying to plot a map with the location of the stations that I am using require(MASS)require(maps) mapaBrasil <- read.table("BASE_BR.txt",sep="") latlonRS <- read.csv2("coordenadas.txt",h=F,sep="") png(filename="mapa_grupos.png",width=600,height=800) # criating a map for RSmap(xlim=c(-80,-35),ylim=c(-45,10),

Full_Name: Richard Calaway Version: 2.9.2 OS: Mac OS X 10.5.3 Submission from: (NULL) (65.47.30.18) Here's my sessionInfo(): > sessionInfo() R version 2.9.2 (2009-08-24) i386-apple-darwin8.11.1 locale: en_US.UTF-8/en_US.UTF-8/C/C/en_US.UTF-8/en_US.UTF-8 attached base packages: [1] stats graphics grDevices utils datasets methods base Currently, if you pass a list

Hi there, When I run the code below I get the error Error in xy.coords(x, NULL, log = log) :(list) object cannot be coerced to type 'double' Any tips how I can resolve this? > > library("waveslim") > > vols=read.csv(file="C:/Users/ocuk/My Documents/Abs Vol.csv", header=TRUE, > sep=",") > x<-c(vols[,1]) > #x > #data(ibm) >

Hola! > rm(list=ls(all=TRUE)) > x <- 1:20 > y <- 1+x+rnorm(x) > xy.coords(y ~ x,NULL) ... expected output, correct, but when called from inside lowess: > lowess(y ~ x) Error in xy.coords(x, y) : x and y lengths differ > debug(xy.coords) > lowess(y ~ x) debugging in: xy.coords(x, y) ... long listing deleted if (is.language(x)) { if (inherits(x,

Dear all, I need to run a interactive procedure where the user will need click on the screen (over a XY plot) and I need to record the XY coordinate which the user clicked. Roughly I wrote a short code below. You see that I suppose that the user will choose four coordinates inside the region of intersection between three segmentes. On each click, I would like to record the position clicked in a

i am working on the project to analyze hedge fund performance, i would appreciate that if you guys could spare some time helping me out with the R code. Thanks. The senario is: i applied BOXPLOT() to plot the performance of all hedge funds with 7 strategies. And right now in this boxplot I need to plot the points of 30 individual hedge funds from my portfolio. And I applied POINTS() and

Hi, I noticed the following problem with xy.coords() in R 2.2.1-patched (version info at the foot of this email) and R 2.3.0 unstable (subversion no: r37123): > xy.coords(x = matrix(1:20, ncol = 2)) Error in xy.coords(x = matrix(1:20, ncol = 2)) : argument "y" is missing, with no default > xy.coords(x = matrix(1:20, ncol = 2), y = NULL) $x [1] 1 2 3 4 5 6 7 8 9

Hello, I want to add text annotation about correlation on "pairs" plots. I found that I could pass a function to the "panel" argument of pairs : panel.annot <- function(x, y, ...) { points(x, y, ...) c <- cor.test(x, y) legend("topleft", legend=substitute(rho == r, list(r=sprintf("%.2f", c$estimate))), bty="n") } And then :

Hi, people. xy.coords() does not behave like its documentation says, when given some matrices. ?xy.coords says: If 'y' is 'NULL' and 'x' is a [...] formula [...] list [...] time series [...] matrix with two columns [...] In any other case, the 'x' argument is coerced to a vector and returned as *y* component [...] Now, consider this short

Is there a way to estimate a nominal response model? To be more specific let's say I want to calibrate: \pi_{v}(\theta_j)=\frac{e^{\xi_{v}+\lambda_{v}\theta_j}}{\sum_{h=1}^m e^{\xi_{h}+\lambda_{h}\theta_j}} Where $\theta_j$ is a the dependent variable and I need to estimate $\xi_{h}$ and $\lambda_{h}$ for $h \in {1...,m}$. Thank you, Mauricio Romero Quantil S.A.S. Cel: 3112231150