[R] xyplot: adjusting the scale (min, max tick)

Hi, I'm terribly sorry but it seems it cannot figure this one out by myself so, please, if somebody could help I would be very grateful. So, when I plot with xyplot() I get an y-axis that is very ugly... starting from a random number and having so many ticks that it becomes unreadable. How do I

Re: [R] xyplot: problems with column names legend

converted, make.names('my variable') [1] my.variable HTH, baptiste 2010/1/3 Jay josip.2...@gmail.com: Hello! one more question about xyplot. If I have data which have space in the column names, say xyz 123. How do I create a working graph where this text is displayed

[R] xyplot: several plots in one creates y-scale problem

Hello, I've been looking for a solution to this problem for some time now but I seem unable to solve it.. so this is the case: I want to plot 4 time series in the same graph using xyplot(). When I do this with xyplot(mydata[,2]+mydata[,3]+mydata[,4]+mydata[,5] ~ mydata[,1], data = mydata,

[R] xyplot: problems with column names legend

Hello! one more question about xyplot. If I have data which have space in the column names, say xyz 123. How do I create a working graph where this text is displayed in the legend key? Now when I try something like xyplot(xyz 123 ~ variable1, data = mydata, ...) I get nothing. Also, is it

[R] Positioning plots on top of each other (aligment borders)

Hello, I want to place two plots on top of each other. However, the problem is that I can't figure out a simple way to align them correctly. Is there a way to specify this? Since the data is bunch of coordinates and the second layer is an outline of a map (a .ps file I import using the grImport

Re: [R] Starting with R and distributions graphics

the computer, pick up a pencil, and solve the problems by hand. Wait a minute, maybe that should be step 1? Have fun, Jay *** G. Jay Kerns, Ph.D. Associate Professor Department of Mathematics Statistics Youngstown State University Youngstown, OH

Re: [R] Need help to forecasting the data of the time series .

response, Stephen; he has seen (most of) this before. Unbelievable. Jay *** G. Jay Kerns, Ph.D. Associate Professor Department of Mathematics Statistics Youngstown State University Youngstown, OH 44555-0002 USA Office: 1035 Cushwa Hall Phone

Re: [R] Combinations and joint probabilities

– 8.50, 8.50 – 10.00), etc, in the first line instead of LETTERS. Welcome to R. Jay P.S. Please don't send HTML. Again, in the Posting Guide. [[alternative HTML version deleted]] *** G. Jay Kerns, Ph.D. Associate Professor Department

Re: [R] Multicore package: sharing/modifying variable accross processes

beautifully with multicore. Feel free to email us with questions, and we appreciate feedback. Jay Original message: Hi, I want to parallelize some computations when it's possible on multicore machines. Each computation produces a big objects that I don't want to store if not necessary: in the end only

Re: [R] Estimation in a changepoint regression with R

strucchange, among others). Jay Message: 4 Date: Thu, 15 Oct 2009 03:56:22 -0700 (PDT) From: FMH kagba2...@yahoo.com Subject: [R] Estimation in a changepoint regression with R To: r-help@r-project.org Message-ID: 365399.56401

[R] Help producing plot for assessing forecasting accuracy

)/(Y_tot) across all values of x. I would be grateful for any guidance you can offer, and I'm sorry if I've overlooked some really simple solution; I'm fairly new to R and learning by doing. Regards, Jay -- Jay Ulfelder, Ph.D. Research Director Political Instability Task Force Science Applications

Re: [R] reading web log file into R

this is cool, great! If not, well... hire a programmer, or if you're lucky Microsoft or Apache have tools to help you with this. There might be something in the Perl/Python world. Or maybe there's a package in R designed just for this, but I encourage students to develop the raw skills... Jay

Re: [R] generating unordered combinations

(prob) urnsamples(0:2, size = 3, replace = TRUE, ordered = FALSE) Regards, Jay *** G. Jay Kerns, Ph.D. Associate Professor Department of Mathematics Statistics Youngstown State University Youngstown, OH 44555-0002 USA Office: 1035 Cushwa Hall Phone

Re: [R] faulty formatting of toLatex(sessionInfo())

discussion last April: https://stat.ethz.ch/pipermail/r-devel/2009-April/053094.html and IIRC this was fixed for R version 2.10. Hope this helps, Jay -- *** G. Jay Kerns, Ph.D. Associate Professor Department of Mathematics Statistics Youngstown

Re: [R] Logistic regression and R

Chapter 5 (and also some of 4). and second, what could I do when I have some independent variables that are not only numerical but categorical too, i.e. mixed (categorical and numerical), can I still use a logistic regression? Easy peasy, lemon squeezy. See page 78. Hope this helps, Jay

Re: [R] kmeans.big.matrix

is right: big data packages (like bigmemory and ff) can't be used directly with R functions (like lm). And because of R's design you can't extract subsets with more than 2^31-1 elements, even though the big.matrix can be as large as you need (with filebacking). I hope that helps. Jay -- John W

Re: [R] Building a big.matrix using foreach

general iteration in parellel. Regards, Jay Original message: Hi there! I have become a big fan of the 'foreach' package allowing me to do a lot of stuff in parallel. For example, evaluating the function f on all elements in a vector x is easily accomplished: foreach(i=1:length(x),.combine=c) %dopar

[R] ' R ' - General Question (newbie)

/ online documentation for the beginning user of R. Thanks, Jay [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting

Re: [R] setdiff bizarre (was: odd behavior out of setdiff)

Jason, (moved back to R-help) On Sat, May 30, 2009 at 3:30 PM, Jason Rupert jasonkrup...@yahoo.com wrote: Jay, I really appreciate all your help help. I posted to Nabble an R file and input CSV files more accurately demonstrating what I am seeing and the output I desire to achieve when

Re: [R] Odd Behavior Out of setdiff(...) - addition of duplicate entries is not identified

posted solution, please specify exactly what you want and I'll wager an R Ninja could dispatch it in moments. Regards, Jay *** G. Jay Kerns, Ph.D. Associate Professor Department of Mathematics Statistics Youngstown State University Youngstown

Re: [R] Odd Behavior Out of setdiff(...) - addition of duplicate entries is not identified

Jason, On Fri, May 29, 2009 at 5:58 PM, Jason Rupert jasonkrup...@yahoo.com wrote: Jay, Thanks much for the reply.    I think you are right about the prob. Unfortunately, I was not able to find the old emails I had discussing the use of the more powerful setdiff that essentially inherits

Re: [R] R Books listing on R-Project

Jay: why not post your R-books how to on the wiki itself??? Because I thought that it would be better to write the instructions in R-wiki language that anybody could modify rather than post a PDF by me. Here is what I had in mind: http://wiki.r-project.org/rwiki/doku.php?id=links:books:howto

Re: [R] can you tell what .Random.seed *was*?

better than this? Conceivably, some of my students will be searching these archives in the future; please feel free to respond off-list if appropriate. Jay -- *** G. Jay Kerns, Ph.D. Associate Professor Department of Mathematics

[R] can you tell what .Random.seed *was*?

to figure it out in 2*10^68 years doesn't count, but within a couple months is acceptable.) 3) does the answer change if there is a remove(.Random.seed) command right before the save.image() command? Any thoughts would be appreciated. Cheers, Jay

Re: [R] can you tell what .Random.seed *was*?

Dear Duncan, On Thu, May 14, 2009 at 4:39 PM, Duncan Murdoch murd...@stats.uwo.ca wrote: On 5/14/2009 3:36 PM, G. Jay Kerns wrote: Question: 1) can you tell me what my original set.seed() value was?  (I wouldn't be able to figure it out, but maybe someone can) The only way I know

Re: [R] Kumaraswamy distribution

View on Probability distributions: http://cran.r-project.org/web/views/Distributions.html and in particular, check out package VGAM. HTH, Jay __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] problem with rgl package

a few days ago, and followed the recommendations of this message: http://tolstoy.newcastle.edu.au/R/e6/help/09/03/9250.html which fixed it for me. HTH, Jay *** G. Jay Kerns, Ph.D. Associate Professor Department of Mathematics Statistics

[R] Bivariate pdf and marginal pdf of non-standard distributions

a suggestion to it. Thanks, Jay Liu, University of Tennesse at Knoxville ## # f(x,y) = 4xy, range of x is (0,1), range of y is (0,1) # Checking if f (x,y) is a joint pdf Pxy-integrate(function(y) { sapply(y, function(y) { integrate(function(x) { sapply(x, function(x) (4*x*y)) }, 0, 1

[R] [R-pkgs] Major bigmemory revision released.

matrices for larger-than-RAM applications. We're working on updating the package vignette, and a draft is available upon request (just send me an email if you're interested). The user interface is largely unchanged. Feedback, bug reports, etc... are welcome. Jay Emerson Michael Kane -- John W

Re: [R] How to set up a function for Central Limit Theorem

- Continuous - Exponential - Sample... If you were wanting to normalize the row sums (or investigate the sampling distribution of some other statistic, for that matter) then check out the Sampling Distributions... menu in RcmdrPlugin.IPSUR. Good luck, Jay Basically I need to prove the Central Limit

Re: [R] setdiff for data.frames?

-devel/archive/26683.html and as a consequence of that discussion: library(prob) setdiff(A,B) Best, Jay *** G. Jay Kerns, Ph.D. Associate Professor Department of Mathematics Statistics Youngstown State University Youngstown, OH 44555-0002 USA

Re: [R] Using very large matrix

Steve et.al., The old version is still on CRAN, but I strongly encourage anyone interested to email me directly and I'll make the new version available. In fact, I wouldn't mind just pulling the old version off of CRAN, but of course that's not a great idea. !-) Jay On Mon, Mar 2, 2009 at 8

Re: [R] Using very large matrix

to email me for more information (and this invitation applies to anyone who is interested in this). Cheers, Jay #Dear friends, # #I have to use a very large matrix. Something of the sort of #matrix(8,8,n) where n is something numeric of the sort 0.xx # #I have not found a way

Re: [R] Bivarite Weibull Distribution

approach, one can generate bivariate samples that have any given marginal distribution. See Nelson for details. Best, Jay thanks A.S. Qureshi __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] generate combination multiset (set with repetition)

way: library(prob) urnsamples(1:3, size = 2, ordered = FALSE, replace = TRUE) You can convert to a matrix with as.matrix(), if desired. Regards, Jay -- *** G. Jay Kerns, Ph.D. Associate Professor Department of Mathematics Statistics Youngstown

Re: [R] generate combination multiset (set with repetition)

carefully about _why_ it is desired to have all possible combinations. In some circumstances, it is good enough to randomly generate combinations and draw inferences from a sampling distribution of some sort associated with the problem. Good luck. Jay

Re: [R] generate combination multiset (set with repetition)

? There isn't a definitive answer to that question. Assuming that it is desired to find an 'average risk' of some sort, it would be good to increase the number of samples until the risk estimate stabilizes to a value with which you are comfortable. Best, Jay

Re: [R] Create unique sets of 3 from a vector of IDs?

? expand.grid(letters[1:5], letters[1:5], letters[1:5]) D Have a look at urnsamples() in the prob package. ID - LETTERS[1:5] urnsamples(ID, size = 3, replace = FALSE, ordered = FALSE) Best, Jay *** G. Jay Kerns, Ph.D. Associate Professor

Re: [R] Create unique sets of 3 from a vector of IDs?

Dear Kingsford, You are quite right, my mistake: urnsamples(ID, size = 3, replace = FALSE, ordered = TRUE) Thanks. Jay On Wed, Dec 3, 2008 at 12:04 AM, Kingsford Jones [EMAIL PROTECTED] wrote: However, I believe Brandon was trying to get the permutations of size 3, rather than combinations

[R] how to test for the empty set

? and in particular, is there an elegant way to check in the case that the mode of the vector is not already known? Thanks in advance for any insight you may have. Best, Jay *** G. Jay Kerns, Ph.D. Associate Professor Department of Mathematics Statistics

Re: [R] how to test for the empty set

Dear Steven, length(x) Does this cover all your use cases? Yes, and thanks again to everybody else who later replied. I had falsely imagined something so much more complicated...!? Next time, I will wait 8*runif(1) before posting. :-) Best, Jay HTH Steven McKinney Statistician

Re: [R] how to test for the empty set

Dear Everybody, Thanks, length() is the answer. Best, Jay On Mon, Nov 24, 2008 at 1:15 PM, Ted Harding [EMAIL PROTECTED] wrote: On 24-Nov-08 17:41:25, G. Jay Kerns wrote: Dear R-help, I first thought that the empty set (for a vector) would be NULL. x - c() x However, the documentation

[R] split data, but ensure each level of the factor is represented

0.2 3 4.7 3.2 1.3 0.2 5 5.0 3.6 1.4 0.2 8 5.0 3.4 1.5 0.2 10 4.9 3.1 1.5 0.1 Thank you for your time, Jay __ R

Re: [R] split data, but ensure each level of the factor is represented

Thanks so much. On Oct 13, 1:14 pm, Henrique Dallazuanna [EMAIL PROTECTED] wrote: Try this: a-factor(c(3, 3, 3, 2, 3, 2, 2, 3, 2, 3), levels = 1:3) split(iris, a) lapply(split(iris, a), dim) On Mon, Oct 13, 2008 at 2:06 PM, Jay [EMAIL PROTECTED] wrote: Hello, I'll use part

Re: [R] Bivariate normal

simply be P( X x | Y=y ) * f(y), where f(y) is the marginal pdf of Y (a dnorm). Note that the above is assuming that y is a fixed constant; if not, then you may want to check out the Ryacas package. I hope that this helps, Jay *** G. Jay Kerns

Re: [R] Plot qnorm

(RcmdrPlugin.IPSUR) and after the Commander restarts take a look at the Distributions menu. You can plot any of the above. In addition, you will have code echoed to you for the console that you can use for other problems. Best wishes, Jay *** G. Jay

[R] call perl

Hi,   It may be the old question. can anyone tell me how to call perl in R? thanks   Y.     [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

Re: [R] Writing Rcmdr Plugins

...[snipped] The /inst directory: in here you will need to put a file menus.txt This should have been the /inst/etc directory, as you noted in your original post. Jay __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r

[R] Setting fixed size for segement plot using stars() (axes size vs print size)

I have been making some segment plots with five variables. They work great, especially when I used a different scale function, which scaled them by area of the circle rather than radius scale - function(x, Mr = 1 , Mx = 100) { ((x/Mx)^.5)*Mr} Where x is the the value, Mr is the Maximum

[R] Help with rpart.predict?

= sites1pabuild) OVEN.pa - predict(OVENrpart.pa, newdata = sites1patest) L3, Plantation, and Range are all character variables, the rest are numeric. Thanks! Jay [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https

Re: [R] SLLN for loggamma fails for 2 order moments

in the actuar package are correct. There may be other people who can identify numerical issues that slow down the convergence rate even more. Best, Jay P.S. the MGF of gamma is all you need for the moments of loggamma... E(X^k) = (1 - k / 2.6)^(-25) On Sat, Jun 14, 2008 at 4:54 AM, michaelhk82 [EMAIL

Re: [R] R package building

created more documentation or a 'white pages' on this, please do spread the word. Thanks to all who have -- and continue -- to work on the system! Jay Subject: [R] R package building In a few days I'll give a talk on R package development and my personal experience, at the 3rd Free / Libre / Open

Re: [R] Multiset Permutations

the columns collapsed, then you can use the apply(B, 1, paste, sep = , collapse = ) command that Johannes suggested. Details are in the prob package vignette, vignette(prob) I hope that this helps, Jay * fix: As it happens, your particular question helped to identify a bug in the current CRAN

Re: [R] R on a computer cluster

on the same computer). Both seem quite easy to use, essentially only needing one command to initiate the cluster and then one command to do something like apply() in parallel. It takes a little planning of your application, but the painfully obvious parallel problem should be painless to implement. Jay

Re: [R] (Small) problem with barchart

Dear Kimmo, It doesn't appear that anyone has yet mentioned pareto.chart() in the qcc package; it may serve your purposes. Please note that it does not require the data frame to be ordered beforehand. x - DATA[[2]] names(x) - DATA[[1]] library(qcc) pareto.chart(x) Best wishes, Jay On Feb 3

Re: [R] Memory problem?

to this is Thomas Lumley's biglm package, which processes the data in chunks. We need more tools like these. Ultimately, you'll need to find some method of analysis that is pretty smart memory-wise, and this may not be easy. Best of luck, Jay - Original message: I am

[R] [R-pkgs] prob package: elementary probability on finite sample spaces

your experiences with it in the classroom and elsewhere. The audience for this package would include teachers and students of elementary probability, or simply anyone wanting to dabble with probability on a finite sample space. Regards, Jay *** G

Re: [R] number of count of each unique row

), FUN = sum) The output will be a data frame with the unique rows, and a column at the end labeled x with the frequency of each unique row. Once you get this you can convert to a list, manipulate, etc. I am sure that there exist faster/better methods. Best, Jay On Dec 21, 2007 5:03 PM, Louis

Re: [R] analytical solution to Sum of binominal distributed random numbers?

Hi Rainer, The distr package can calculate the distribution for you: library(distr) X - Binom(size = 7, prob = 0.3) Y - Binom(size = 11, prob = 0.5) Z - X + Y d(Z)( 0:18 ) # the pmf r(Z)( n = 5 ) # random variates Please note, however, that size and prob must be of length 1. Best, Jay On 10

<    1   2