search for: ehler

Full_Name: Peter Ehlers Version: 1.6.0 OS: Windows 2000 Pro Submission from: (NULL) (136.159.61.178) In: p.adjust package:base R Documentation In the paragraph: Hochberg's and Hommel's methods are valid when the hypothesis tests are independent or when they are non-negatively a...

Pat Burns makes a good point. -Peter -------- Original Message -------- Subject: Re: [R] Using seq_len() vs 1:n Date: Fri, 12 Feb 2010 09:01:20 +0000 From: Patrick Burns <pburns at pburns.seanet.com> To: Peter Ehlers <ehlers at ucalgary.ca> References: <4B746AEF.10900 at ucalgary.ca> If you want your code to be compatible with S+, then 'seq_len' isn't going to work. On 11/02/2010 20:39, Peter Ehlers wrote: > R-people, > > Duncan Murdoch's response in > > https://sta...

How can I add an overall plot title to these four plots? I would like to have something that says, "Distribution Comparisons": par(mfrow = c(2, 2)) # Plot 1 plot(rnorm(10),type="l",col="red") title(main = list(paste("Normal"), ????????????????????????? col="black", cex = 1.0)) # Plot 2 plot(rpois(10,

...le values are non -ve), how about: > > x <- runif(20, -1, 100) > exp(median(log(pmax(0,x)))) > > It'll give -Inf if the two middle values are negative, when I guess we > should get NaN, but I can't see a 1-line way to handle that! > > Keith J > > "Peter Ehlers" <ehlers at ucalgary.ca> wrote in message > news:4D3468EF.5010601 at ucalgary.ca... >> I've been reminded by Prof. Brian Ripley that R's >> log() function will indeed handle zeros appropriately. >> >> Apologies to S Ellison and Hadley Wickham. >>...

I have seen some strange error log messages in our dovecot server (2.0.15, RHEL 5.7x64, mdbox+zlib, local storage, high load) Dec 14 22:01:27 buzon dovecot: imap(fmcarrasco): Error: /indices_dovecot/indices/07/47/fmcarrasco/mailboxes/INBOX/dovecot.index reset, view is now inconsistent Dec 14 22:00:27 buzon dovecot: imap(cfp): Error:

...g the wrong argument type; nevertheless, if the behaviour is general, it might be better to generate an error. For anyone with more know-how than I have, here are the details of the crash: AppName: rgui.exe AppVer: 2.1.41115.0 ModName: r.dll ModVer: 2.1.41115.0 Offset: 000ef1b3 Peter Ehlers U of Calgary ehlers@math.ucalgary.ca

Hello all, Is it possible to change the radii line type in radial plots? I wasn't able to find anything online. Thanks, Bern -- View this message in context: http://r.789695.n4.nabble.com/Changing-radii-line-type-in-radial-plots-tp4647238.html Sent from the R help mailing list archive at Nabble.com.

I have a new device that takes measurements anywhere from every second, to every 15 minutes (depending on changes). The matrix has a date, time and Y column (Y is the measurement). For three days it is 25,000 rows. How do I average the measurements by every 30 minutes so my matrix is 48 rows per day? I have been working on this and cannot figure out a simple method. Any ideas? Thank you. ----- In

#I wish to find the "next non-NA" value within each row of a data-frame. #e.g. I have a data frame mydata. Rows 1, 2 & 3 have soem NA values. mydata <- data.frame(matrix(seq(20*6), 20, 6)) mydata[1,3:5] <-  NA mydata[2,2:3] <-  NA mydata[2,5] <-  NA mydata[3,6] <-  NA mydata[1:3,] #this loop accomplishes the task; I am tryign toi learn a "better" way for(i

...levels(f))[f], the vector as.numeric(levels(f)) > > is indexed by as.integer(f). > > > > This appears to rely on the current implementation, as mentioned in > section > > 2.3.1 of the language definition. > > > > > > On Mon, Apr 1, 2013 at 1:49 PM, Peter Ehlers <ehlers@ucalgary.ca> wrote: > > > > > On 2013-04-01 10:48, Matthew Lundberg wrote: > > > > > >> These two seem to be at odds. Is this the case? > > >> > > >> From help(factor) - section Warning: > > >>> > > &g...

...ns) that the compiler can't tell you when the arguments are wrong. On a number of occasions it's taken far more time than it should to find errors of that sort. The convenience wasn't worth the wasted time. That's my rant for today. *-- Russ * On Thu, May 5, 2011 at 10:13 AM, P Ehlers <ehlers@ucalgary.ca> wrote: > Russ, > > All you have to do is replace > > > addTA(GSPC.EMA.3, on = 1, col = "#0000ff") > > with > > plot(addTA(GSPC.EMA.3, on = 1, col = "#0000ff")) > > etc. > > I can sympathize with the doc...

Hello, I have a simple question that I could not really figure out. I am plotting labels within a graph using the text function. I first plot the first label by specifying the x and y coordinates on the graph. Then to plot the second label next to it, I am using te strwidth function to get the width of the previous label in user coordinate then add the maximum width to the x value. However, this

Full_Name: Uffe H?gsbro Thygesen Version: 2.2.0 OS: linux Submission from: (NULL) (130.226.135.250) Hello all. pbinom(q=0,size=0,prob=0.5) returns the value NaN. I had expected the result 1. In fact any value for q seems to give an NaN. Note that dbinom(x=0,size=0,prob=0.5) returns the value 1. Cheers, Uffe

Hi list, One of my co-workers found this problem with 'cor' in his code and I confirm it too (see below). He's using R 2.2.1 under Win 2K and I'm using R 2.3.0 under Win XP. =========================================== > R.Version() $platform [1] "i386-pc-mingw32" $arch [1] "i386" $os [1] "mingw32" $system [1] "i386, mingw32" $status

....exe http://cran.us.r-project.org/R-2.3.0dev-win32.exe (Same for all other mirrors (about 6) that I've tried. Could someone please check if these files are accessible to them? I'm guessing there's something haywire in my setup. Note that I can ftp the files without problem. -- Peter Ehlers Department of Mathematics and Statistics University of Calgary, 2500 University Dr. NW ph: 403-220-3936

...ts which encourage the "single table with a gargantuan number of columns, most of which are missing data in most cases" approach to organization of longitudinal data is regrettable." http://n4.nabble.com/Hierarchical-data-sets-which-software-to-use-td1458477.html#a1470430 -- Peter Ehlers University of Calgary

Dear Community, I have a matrix with assigned colnames and rolnames as follows: A B NR 0.15 0,05 AL 0,05 0,05 . . . . . . . . . I want to extract the names of the rows for which A>0,1 and B<0,1. In the above example this would be observation NR only. Hence the output should write for

? stato filtrato un testo allegato il cui set di caratteri non era indicato... Nome: non disponibile Url: https://stat.ethz.ch/pipermail/r-help/attachments/20060525/ef972e9f/attachment.pl

Readers, I tried to the following commands: plot(y~x,ylab=expression(A[1]~B[2],xlab=expression(C~D)) mtext(expression(A[1]~B[2]),"additional text",side=3,line=1) I receive the text that I want, but the command terminal shows the following response: Warning message: NAs introduced by coercion in: mtext(expression(A[1]~B[2]),"additional text", side = 3, What is my mistake

Full_Name: foo ba baz Version: R2.2.0 OS: Mac OS X (10.4) Submission from: (NULL) (219.66.32.183) chisq.test(matrix(c(9,10,9,11),2,2)) Chi-square value must be 0, and, P value must be 0 R does over correction when | a d - b c | < n / 2 &#65292;chi-sq must be 0