similar to: generic/method consistency

Hi, I am having problems creating a pdf file of a lattice graph. Things work fine for a single image, but I am having trouble using the commands in a loop. To illustrate by example This works with both ps and pdf files. dat = list(x= 1:10, y = 1:10) trellis.device(postscript, file = "/d1/pocernic/test.ps") xyplot(y~x, data = dat) dev.off() This does not. It produces a very

Hi, I was lately debugging parts of my 'colAUC' function in caTools package, and in a process looked into other packages for calculating Areas Under ROC Curves (AUC). To my surprise I found at least 6 other functions: * wilcox.test * AUC from ROC package, * performance from ROCR package, * auROC from limma package, * ROC from Epi package, * roc.area from verification

Hi, I am trying to plot different variables from a data.frame using lattice's xyplot using code like that below. How do I specify a symbol and color for the variable 'prob' and different one's for 'll.prob'? Thanks, Matt xyplot( prob + ll.prob ~ time.eff |stat.id + time.out ,data = OUT, allow.multiple = TRUE, layout = c(6,3), as.table = TRUE ,

Hello, I have been having difficulties opening netcdf files using the netcdf library. I am able to successfully open the file and I can determine the names and size. However, when I try to read it I receive an error. Reading the documentation, I am wondering if this might be due to a difference between netcdf version 2 and more recent netcdf format? I've tried using different values for the

Hi, Using lattice's xyplot, is it possible to specify a variable to group plots by page. For example, if I have xyplot(y~x|A*B*C) could I get a page created for each unique value in variable C ? I am hoping to avoid having pages with the same strip above each plot on a page. Thanks, M Matt Pocernich NCAR - Research Applications Program 303-497-8312

Hi, How do I convert a numeric value indicating the time since 1970, back into a POSIXct class object? I have tried format.POSIXct and as.POSIXct without success. For example > ccc [1] "1945-01-01 15:00:00 MDT" > ddd<- as.numeric(ccc); > ddd [1] -788842800 > format.POSIXct(ddd) Error in format.POSIXct(ddd) : wrong class > as.POSIXct(ddd) Error in

The verification package has recently been posted to CRAN. This package was initially developed by people in the Verification Group at the National Center for Atmospheric Research to verfiy and study weather models and forecasts. It has been written in general terms to be applicable to other fields of study. Functions include receiver operating characteristic curves attribute diagrams

The verification package has recently been posted to CRAN. This package was initially developed by people in the Verification Group at the National Center for Atmospheric Research to verfiy and study weather models and forecasts. It has been written in general terms to be applicable to other fields of study. Functions include receiver operating characteristic curves attribute diagrams

I'm sure this question is answered in the help file, but likely I'm not reading it corrected. Running windows version 1.00.0, loading a table (35K rows by 10 columns) from Excel using the read.table command I receive the following message. Error: cons memory (350000 cells) exhausted See "help(Memory)" on how to increase the number of cons cells. >From reading the

Hi - I would like to create a scatterplot with the symbols coded by a factor in each record. Is their a better way than what I have used below? I attempted to use the plot.factor, but cannot force it to make a scatterplot instead of boxplots. I have considerd making a function with a 'for' statement for each level of the factor. Is there a better way for when I am dealing with more

Hello, We previously had Win2008 R2 DCs which we have begun to replace with Win2012 servers (forest is still at 2008 level however). I was able to query DFS via rpcclient but now I find that I am unable to. I know that SMB3 was introduced in Win2012 but I believe that it should auto-negotiate earlier versions. I did however test from v4.1.9 as well, unfortunately that still fails. From Samba3

Dear Geskus, I want to develop a prediction model. I followed your paper and analysed thro' weighted coxph approach. I can develop nomogram based on the final model also. But I do not know how to do internal validation of the model and subsequently obtain calibration plot. Is it possible to use Wolbers et al Epid 2009 approach 9 (R code for internal validation and calibration) . It is

Are there any R function or packages that generate wind roses? These are ``star''-like plots that summarize the frequency, direction and magnigude of windspeeds. Thanks, Matt -- Matt Pocernich University of Colorado - Denver Department of Mathematics -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hello, I am having difficulty filtering data. I am working with flow data collected at a stream gage. For each record, I have a date and flow value. I have filtered this data to only include days when flow values exceed a given threshold. Here is my problem. Within this subset of data, I often have several consecutive days above the threshold. From this group of days, I wish to select the

Dear Thomas, I want to use evaluate effect of Age using restricted cubic form in the FGR function as Fgr.crr <- FGR(Hist(time, event) ~ rcs(Age_years), data=dat) It provides error. " Error in parse(text = termtext, keep.source = FALSE): .... 1: response ~ rcs(Age_years Do I need to change any of the R code? Regards Amalraj Raja The University of Aberdeen is a charity

Dear All, I would like to use R function 'selectFGR' of fine gray model in competing risks model. I used the 'Melanoma' data in 'riskRegression' package. Some of the variables are factor. I get solution for full model but not in variable selection model. Any advice how to use factor variable in 'selectFGR' function. The following R code is produced below

Hello, I have a question about using the layout command within a function. I've written function that uses layout to create a figure from 2 plots. This works fine to create a figure. When I use par(mfrow = c(2,2)) to create multiple plots, it seems that the layout command resets the mfrow parmeter. Is there a way for me to avoid this problem? For example practice<- function() {

Hi , library(foreign) write.dta(data1, "data1.dta") should work. The file will be saved in the working directory. Use getwd() to know the working directory. Best wishes Amalraj Raja -----Original Message----- From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of rosario scandurra Sent: 22 March 2018 07:47 To: r-help at r-project.org Subject: [R] exporting data to stata

I have a simple question. On a plot, I would like to place text containing a calculated value. What is the format? For example. I would like the following text to appear on a plot. The average temperature is 23.1 degrees. I'm guessing the command is text( x, y, "The average temperature is " %variable name% "degrees") Thanks, Matt Pocernich -------------- next part

I am having problem using the step function for a linear regression model. I've created an initial model containing only the intercept. Then using the step function, I've selected three variables to be considered for the model. > x0.lm<- lm(MEDV~1, data = x) > > anova(x0.lm) Analysis of Variance Table Response: MEDV Df Sum Sq Mean Sq F value Pr(>F)