search for: wid

Hi everyone, I would like to identify the case by groups that is just bigger that avg plus sd. For example, using species as group and petal.wid as my variable in the iris data. What's the better way to doit? creating a function? So,the question is to identify the single element of each species that is just larger than a cut-off point (i.e. larger than mean + sd) I made this, but I can not make a relation to to orginal data for identi...

...f step[i-1] then it is also independent of -step[i-1]): Cheers Joe > # your (unbounded) random walk > k <- cumsum(c(0,sample(c(-1,1), 1000, rep = T))) > > #shift it to positive, to allow following calculations k <- k - min(k) > # set your bound (remember min(k) is now zero) wid <- 10 > > # fold the walk back on itself with 'creases' at wid, 2*wid, 3*wid, ... > walls <- ifelse((k%/%wid)%%2==1, > # if k[i] = a*wid + b, with a odd and b < wid, > # invert it and map it to (0, wid) > (1+(k%/%wid))*wid...

...clear how what data object to put the standard errors for the thing to work... My code (the version that at least draws a plot): xyplot(inbreeding~generations|breeds, data = trellisplot, type = 'l', layout = c(2,5), col = 'black', lwd = 3, errbar = function(x, y, err, down = T, width = 0.005, lwd = 1, plot = T, colour="black") { up <- y + err dn <- y - err x.range <- max(x) - min(x) wid.lf <- x - (x.range * width) wid.rt <- x + (x.range * width) if(plot == F) { return(data.frame(x, up, dn)) } else { if(down == F) {...

Dear all, I am a total beginner in R, so sorry if this is the wrong place. I am using R 2.10.1 on a Mac (Mac OS 10.6.2). I have this small dataset : growth sugar 75 C 72 C 73 C 61 F 67 F 64 F 62 S 63 S I have no problem reading the table, or getting the summary, but if I try boxplot(growth~sugar, ylab="growth", xlab="sugar", data=Dataset), I have the following error :

...Stub; we set our timer resolution at minimum warn:dsound:DirectSoundDevice_SetCooperativeLevel level=DSSCL_PRIORITY not fully supported warn:dsalsa:SetFormat Could not set sound rate to 44100, but instead to 48000 warn:dsound:PrimaryBufferImpl_GetVolume control unavailable warn:mmtime:timeKillEvent wID=3 is not a valid timer ID warn:mmtime:timeKillEvent wID=3 is not a valid timer ID warn:heap:HEAP_Decommit Could not decommit 003d0000 bytes at 0xbf970000 for heap 0xbf940000 warn:heap:HEAP_Decommit Could not decommit 003d0000 bytes at 0xbf970000 for heap 0xbf940000 warn:file:wine_nt_to_unix_file_na...

Dear R-Help members, I have a question about how to save to the hard drive the one thousand datasets I generated in a simulation. The datasets are created in a "for" loop that repeatedly creates normally distributed datasets, such as the example below: Library(MASS) for (number in 1:1000) { dataset = mvrnorm(n = 400, mu = c(0,0,0), Sigma =

...hes have been truncated. ?? what should the behavior be ?? the diff is against the 11 Jan version of R 2.3.0 *** newboxplot.R 2006-01-23 14:32:12.000000000 -0500 --- oldboxplot.R 2006-01-23 14:29:29.000000000 -0500 *************** *** 84,98 **** bplt <- function(x, wid, stats, out, conf, notch, xlog, i) { ## Draw single box plot - conf.ok <- TRUE - if(!any(is.na(stats))) { - ## check for overlap of notches and hinges - if (notch && (stats[2]>conf[1] || stats[4]<conf[2])) { - conf.ok &lt...

...<- as.factor(data$Participant) summary(data) colwise(class)(data) colwise(mode)(data) colwise(mean)(data[3:4]) colwise(sd)(data[3:4]) ops <- options(contrasts = c('contr.sum', 'contr.poly')) ezANOVA(data, dv = .(Measure.2), wid = .(Participant), within = .(Factor), type = 3, detailed = T) ezANOVA(data, dv = .(Measure.1), wid = .(Participant), within = .(Factor), type = 3, detailed = T) options(ops) I get a model from the first call to ezANOVA, but the deficient rank error from the second (script...

Hi Everyone: I was having trouble getting the plot size correct when putting a plot in an sweave document. So, I searched on the internet and I found a clever solution where the person uses the cat function to write latex code so he can include the cat code in the chunk and place the plot using width and height. The solution is here: https://stat.ethz.ch/pipermail/r-help/2009-April/195094.html. But the solution uses a backslash before the ending quotes in the cat statements that might be windows relevant possibly ? I'm not sure about that part but I'm on linux and when I tried my wa...

...plotting more than one page on png()? 'Properly' means to me with leading 0s, so that the PNGs become easy to navigate in a file/image browser. Lacking a better solution I ended up using the code below, but would much prefer something like png("test_%d.png",bg="white",width=1000,height=700) where %d could be formatted like formatC(%d,digits=0,wid=3,flag="0",mode="integer") Thank you, b. #---works, but is rather complicated--- pngno <- 0 ; i <- 1 for (w in 1:53) { if (i %in% c(4*0:100+1)) { pngno <- pngno + 1 png(paste(&quo...

What is the likelyhood that major inexpensive webhosts like godaddy, bluehost, etc. will upgrade to RoR 1.1? Is this going to be like PHP 5 where it has to percolate for a year or more before it becomes widly available? Your thoughts? Along the same lines... is it possible to adopt some of the new improved Ajax / javascript capabilities without actually upgrading the ruby installation? Thanks in advance. -- Posted via http://www.ruby-forum.com/.

Prasad wrote: > I wrote a function in R which uses tcltk package .... essentially I wanted > to give within that function, a widget with 2 radiobuttons to choose > between plotting Precip and Temperature plots. After the user has chosen > one of the radiobuttons there is another widget that asking him to identify > outliers. However, I am having a lot of problems...what R does is evaluate > the whole function wit...

Hello, I wrote a function in R which uses tcltk package .... essentially I wanted to give within that function, a widget with 2 radiobuttons to choose between plotting Precip and Temperature plots. After the user has chosen one of the radiobuttons there is another widget that asking him to identify outliers. However, I am having a lot of problems...what R does is evaluate the whole function without pausing...if I...

...x file. I needed help from someone with more latex experience to progress Duncan At 10:12 20/08/2009, you wrote: >Ken > >Thanks for your inspiration: This combination of strwrap and writeLines >seems to "work" in one sense of the word: > > > ># Intro stuff: >wid <- 50 >option(width=wid) >hide <- capture.output(example(glm)) > ># The business end: >writeLines( strwrap( capture.output(summary(glm.D93)), wid)) > > > >Mind you, while it "works", it's an awful construct. I was actually >hoping to show this cod...

.... KevinW> I seem unable to left-justify exponential format KevinW> numbers. There appears to always be an extra space KevinW> inserted to the left. KevinW> Using the example from the formatC help page: R> xx <- pi * 10^(-5:4) R> cbind(formatC(xx, wid = 9, flag = "-")) KevinW> [,1] KevinW> [1,] " 3.142e-05" KevinW> [2,] "0.0003142" KevinW> [3,] "0.003142 " KevinW> [4,] "0.03142 " KevinW> [5,] "0.3142 " KevinW> [6,] &qu...

...le size is different(M:F=23:18) and one of condition(cond) and the interaction at each subset from 90 datasets So I perform anova 90 times using a command like below; for(i in 1:90) {results_ezANOVA = ezANOVA(data=subset(ast.ast_coef, ast.ast_coef$coef_thr==i), dv=.(ast.values), between=.(gender), wid=.(subj), within=.(cond))} But I got the last(90th) result, not all. Here are my questions. 1) Is my command correct? 2) If correct, please let me know if I can get all 90 results. 3) What kind of postHoc would be appropriate? Thank you, Jeong [[alternative HTML version deleted]]

...s something I haven't understood about the grid.layout function. The second time I use it, it does what I expect, but there's something wrong with the first one. require(grid) N <- 5 variates <- LETTERS[1:N] grid.newpage() plot.square <- 2 * length(variates)# cm text.wid <- 3 # cm vp <- grid.layout(nrow = 1, ncol = 2, widths = c(text.wid, plot.square), heights = plot.square, default.units = "cm") pushViewport(viewport(layout = vp, layout.pos.row = 1, layout.pos.col = 2)) ## vpp...

...;plummer@iarc.fr> >> To: r-testers@stat.math.ethz.ch >> #4) If you put a few lines of comments at the start of a function, >> then the comments always appear *after* the first line of code. This is REALLY painful; try, e.g., fcat <- function(..., f.dig= 4, f.wid = f.dig +5, f.flag = ' ', nl = TRUE, file = "", sep = " ", fill = FALSE, labels = NULL, append = FALSE) { ## Purpose: Formatted CAT -- for printing 'tables' ## -----------------------------------------------------------------...

...n large matrices, it's a familiar truism that as.vector(z %*% rep(1, dim(z)[2])) is MUCH faster than apply(z, 1, sum). The equivalent "inline" code for "max" is not so nice. So here's an alternative version (untested): z <- abs(outer(w, x, "-")) wid <- seq(along=w) top <- 1 + round(max(z), 0) # we know z is nonnegative by construction s <- z[ sort.list((z / top) + matrix(wid, nrow=dim(z)[1], ncol=dim(z)[2], byrow=F)[wid * dim(z)[2]] ] z <- (1 - (z / matrix(s, nrow=dim(z)[1], ncol=dim(z)[2], byrow=F))^3)^3 sums <- as.vect...

...ny device in that vendor/device range would bypass the iommu unless some new property "actually-works-like-a-real-pci-device" happens to exist in the device-tree, which we would then need to define somewhere and handle accross at least 3 different platforms who get their device-tree from widly different places. Also if tomorrow I create a PCI device that implements virtio-net and put it in a machine running IBM proprietary firmware (or Apple's or Sun's), it won't have that property... This is not hypothetical. People are using virtio to do point-to-point communication bet...