similar to: lattice, different plotting symbols

Hi, I realize the answer is very likely in the section Generic functions and methods (or Adding new generics), but I'm not clear what to do with the following. Running R CMD check, I get the following warnings for my generic functions. Does this mean I need the argument * checking generic/method consistency ... WARNING leps: function(x, ...) leps.default: function(x, pred, titl,

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, 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, 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

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, 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

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

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

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() {

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

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)

For block effects with small variance, lmer will sometimes estimate the variance as being very close to zero and issue a warning. I don't have a problem with this -- I've explored things a bit with some simulations (see below) and conclude that this is probably inevitable when trying to incorporate random effects with not very much data (the means and medians of estimates are plausibly

Me funciono perfecto, muchas gracias Carlos ... si no es mucha la molestia, alguien podria explicarme la diferencia tecnica entre ambos signos, y porque == produce el resultado que produce ? en que casos deberia usar == y en que casos %in% ... cual es el significado de cada uno ? Disculpen que no me conforme solo con obtener el resultado adecuado, es solo que no me suena logico que == no

Estimados, tengo el siguiente problema: estoy trabajando con un data.frame que tiene 892 filas y 9 columnas, el cual adjunto para que puedan probar. Cuando consulto por las filas que cumplen una cierta condicion, por ejemplo, ... graph[graph$NCar==c("160s"),] obtengo el siguiente conjunto de datos: sol con dia NCar NIns isom area rep tipo 1 con 0.001 1 160s 0 s

Muchas gracias por la explicación Isidro! El 10 de septiembre de 2013 08:19, Isidro Hidalgo <ihidalgo@jccm.es>escribió: > Cuando utilizas "==" estás comparando el vector "graph$NCar" (longitud 892) > con otro de longitud 2: c("160s","180s"). > Es decir, comparas el primer valor de "graph$NCar" con el primer valor de >

I have a process where I am creating a effects plot similar to the cowles effect example. I would like to add the point estimates to the effects plot, can someone show me the correct syntax. I have included the "R" effects example, so you can show me the correct syntax. Thanks mod.cowles <- glm(volunteer ~ sex + neuroticism*extraversion, data=Cowles, family=binomial)