search for: pocernich

...e 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 levels of factors? Thanks, Matt Pocernich Y <- c(40,38,30,47,25,32,13,35,27,24,20,13) X<-c(110,75,93,97,60,75,38,140,62,90,45,59) manuf <- factor(rep(LETTERS[1:3],c(4,4,4))) data3 <- data.frame(X,Y,manuf) attach(data3) plot(X, Y, type="n") points(X[manuf=="A"], Y[manuf=="A"], pch = "A")...

...eneric/method consistency ... WARNING leps: function(x, ...) leps.default: function(x, pred, titl, plot) reliability.plot: function(x, ...) reliability.plot.default: function(x, obar.i, prob.y, titl, mod.names) summary: function(object, ...) summary.prob.bin: function(object) Matt Pocernich NCAR - Research Applications Program 303-497-8312

...ory)" on how to increase the number of cons cells. >From reading the help file, I've learned a bit more about static memory, but I'm not sure where to place the R --vsize=15M --nsize=1000k line. Where should this line be placed? (Or in which file?) Thanks for your time. Matt Pocernich -------------- next part -------------- An HTML attachment was scrubbed... URL: https://stat.ethz.ch/pipermail/r-help/attachments/20000317/5ffd82eb/attachment.html

...he preceeding and succeeding day is below the threshold) I also wish to select that record. For example from the data set Day Flow 1 10 4 13 5 20 6 15 9 13 I would like the 1st, 3rd and 5th record filered. Any ideas on how I would write such and algorithm would be appreciated. Thanks, Matt Pocernich -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat...

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 http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "...

...thing like this to produce a large, appended file made from different large datasets. This doesn't work either. trellis.device(postscript, file = "/d1/pocernic/test.ps", onefile = TRUE) for(i in 1:5){ xyplot(y~x, data = dat) } dev.off() Any thoughts would be welcomed. Matt Pocernich NCAR - Research Applications Program 303-497-8312

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

...e'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 , panel = function(x,y){panel.abline(h = 0) panel.xyplot(x,y)}, strip = TRUE, ) Matt Pocernich NCAR - Research Applications Program 303-497-8312

...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 -------------- An HTML attachment was scrubbed... URL: https://stat.ethz.ch/pipermail/r-help/attachments/20000202/b943cd4c/attachment.html

...-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 as.POSIXct.default(ddd) : Don't know how to convert `ddd' to class "POSIXct" > Thanks, Matt Matt Pocernich NCAR - Research Applications Program 303-497-8312

...ng to predict) is not found. I receive an identical message using the command > step(x0.lm). I am using R version 1.1.1. I can run all the examples successfully. Also, I realize that there are differences in the algorithms, but I can run the code above on Splus version 5. Thanks Matt Pocernich -------------- next part -------------- An HTML attachment was scrubbed... URL: https://stat.ethz.ch/pipermail/r-help/attachments/20001129/2f219ad5/attachment.html

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

I am having a series of problems using date time data that has been converted into a POSIXt and POSIXlt classes. I have hourly time series data from 1900 that has been converted from text data. I assume most of my problems come from a mis-underdanding of the POSIX class. My matrix named (aa) for this year is approx 8700 by 4. When I try to calculate the length of posit column ( which is the

...ore. Suggestions and comments on the package can be directed to me at pocernic at ucar.edu. More information on verification (with a somewhat meteorolgical perspective) can be found at http://www.bom.gov.au/bmrc/wefor/staff/eee/verif/verif_web_page.html. Thanks, Matt ##################### Matt Pocernich NCAR - Research Applications Program 303-497-8312 _______________________________________________ R-packages mailing list R-packages at stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-packages

How do I change the dimensions and proportions of my graphics device? I've used the fin parameter to change the shape of the figure region, but the files are still square. I'm using the windows version R 0.90.0. Thanks, Matt -------------- next part -------------- An HTML attachment was scrubbed... URL:

I am producing a series of histograms and would like to control the format of the tick labels. I tried using the exp=1 option, but apparently this can't be set in a higher level plot function. I've unsuccessfully tried to modify the hist command. Any suggestions? Thanks Matt Pocernich -------------- next part -------------- An HTML attachment was scrubbed... URL: https://stat.ethz.ch/pipermail/r-help/attachments/20000216/e9165757/attachment.html

...ot; "single" "single" "single" [22] "single" "single" "single" "single" ... etc ## when I try to read it, I get the following error. > bbb<- read.netCDF(aaa) Error: dim: Invalid dimension vector Thanks Matt Matt Pocernich NCAR - Research Applications Group 303-497-8312

...) nn<- layout(matrix(c(0,1), nrow = 1), widths = c(1,4)) plot(0,0) mtext("Sample graph", outer = TRUE, line = 1, cex = 2) par(op) } practice() # works fine. par(mfrow = c(2,2)) practice() practice() practice() practice() ## produces 4 graphs on seperate pages. Thanks. Matt Matt Pocernich NCAR - Research Applications Program 303-497-8312

...ore. Suggestions and comments on the package can be directed to me at pocernic at ucar.edu. More information on verification (with a somewhat meteorolgical perspective) can be found at http://www.bom.gov.au/bmrc/wefor/staff/eee/verif/verif_web_page.html. Thanks, Matt ##################### Matt Pocernich NCAR - Research Applications Program 303-497-8312 _______________________________________________ R-packages mailing list R-packages at stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-packages

On Wed, 29 Nov 2000, Matt Pocernich wrote: > 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, dat...