similar to: multiple selection syntax

Dear R-help: Does anyone know of a package in R that will do Welch's modified F-test for heterogeneous error variances? Are there other statistical techniques available in R that test the equality of means when homoscedastisity is violated? 't.test' does this in the pairwise sense when var.equal = TRUE. With best wishes and kind regards I am Sincerely, Corey A. Moffet Support

Dear R-Help, Does anyone know of a way to use CGIwithR in Windows? David Firth the author of CGIwithR states "It ought to be possible to adapt it for use also in conjunction with web servers running under Windows NT ... but the author has no plans for that." I have developed some applications in R that I would like to make available through a web interface using CGIwithR but at this

On Thu, Jan 7, 2016 at 1:35 PM, Krzysztof Parzyszek <kparzysz at codeaurora.org > wrote: > On 1/7/2016 3:25 PM, Phil Tomson wrote: > >> >> That's better, but now I get: >> >> XSTGInstrInfo.td:902:3: error: In RelAddr: XSTGRELADDR node requires >> exactly 2 operands! >> >> Which makes some sense as XSTGRELADDR is defined as: >> def

The funciton c() works differently for strings and for factors: For strings: > l = c('a', 'b') > l [1] "a" "b" For factors: > l = c(factor('a'), factor('b')) > l [1] 1 1 What should be the right technique for merging factors? -- Svetlana Eden Biostatistician II School of Medicine

I'm sure this is a total noob question, but half an hour of searching bore no fruit: How do you select a subset of a vector by negative value? If I try > hist(sp$p[sp$r>0.01]) # all is well > hist(sp$p[sp$r<-0.01]) # this obviously causes issues Also, putting -0.01 in parentheses didn't help. Thanks in advance -- b

Is there a function in R, which would return index of maximum value in a vector ? e.g. > v <- round(10*rnorm(8)) > v [1] 6 -3 -6 15 7 9 0 -19 > max(v) [1] 15 ??? index.max(v) ??? 4

Hi sorta piggybacking on the other thread. I am looking for some feedback on how to implement the following idea in llvm. The really short version of the idea is this: * I want to alloca a field (record/struct), so that its size is an even multiple of 64 bytes. [^1] * This allocaed field will be exclusively used as an argument to functions * llvm should be aware of the extra bytes and should

Hi I wrote a function implementing an algoritm for simulation of spatial sampling plans (using geoR). I'm using a range of parameters to perform several simulation and sometimes this parameters create function errors wich cause the algoritm to stop. How can I avoid this ? Below there's a small example that I made to study the problem but I was not able to solve it. I've tried

Hi, Does anyone know how to manually configure the number of digits printed out from xtable()? For example, I'm exporting a data frame through xtable() into a LaTeX table, I only have two columns in the data frame so by default I only get two decimal places. But I'd like at least 5 decimal places. I had a look at ?xtable() but can't seem to find an example. Cheers, Kevin

Probably a stupid question, but I don't seem to be able to find the answer I'm looking for from any of the R literature. Basically I have a matrix with several thousand rows and 20 columns(weather stations) of wind direction data. I am wanting to extract a matrix which contains data for all columns conditional on column 20 having a value of _either_ less than 45 or greater than 315. (ie I

Hi R experts, I'm trying to use nls() for a piecewise linear regression with the first slope constrained to 0. There are 10 data points and when it does converge the second slope is almost always over estimated for some reason. I have many sets of these 10-point datasets that I need to do. The following segment of code is an example, and sorry for the overly precise numbers, they are just

We'd like to be able to vary the latency of our load instructions based on what address space is being loaded from. I was thinking I could do this by overriding getOperandLatency in our target, but I'm wondering if the addrspace info is available when instructions are scheduled? For example, I have this in our llvm IR: %0 = load i32 addrspace(4)* @answer, align 4 store i32 %0, i32*

Hi everyone, How might I go about configuring R to keep add-on packages loaded from session to session? Is this undesireable for some reason? At present, I keep a file called "pckgs.txt" in my working directory with, e.g., library(gnlm) library(Hmisc) ... and then type source("add-ons.txt") every time I start a new session. I suspect there's a more elegant way to hold

Another newbie question.... I want to create a factor (say cT)from a numerical variable (sy temp) by regrouping the values in classes (say cT <390, [390,400[,[400,409[...>=550) Is there a simple way of doing that using the factor function? AND I do not find how to manipulate strings (I want to concatenete characters strings ("abkdas","chjw") into something like

I am having difficulty with plot.formula. For example: ------------------- dat <- data.frame(y = c(1,1,1,2,2,2,3,3,3,4,4,4), x1 = as.ordered(c(1,1,1,1,1,1,2,2,2,2,2,2)), x2 = as.ordered(c(1,1,1,2,2,2,1,1,1,2,2,2))) plot(y ~ x1 + x2, data = dat) # Works fine "interactively" asking: # "Hit <Return> to see next plot:" #

Hello, I'm a newcomer to R so please forgive me if this is a silly question. It's that I have a linear regression: fm <- lm (x ~ y) and I want to test whether the slope of the regression is significantly less than 1. How can I do this in R? I'm also interested in comparing the slopes of two regressions: fm1 <- lm (x ~ y) fm2 <- lm (a ~ b) and asking if the slope of fm1 is

I know that there must be a cool way of doing this, but I can't think of it. Let's say I have an dataframe with NA's. > x <- data.frame(a = c(0,1,2,NA), b = c(0,NA,1,2), c = c(NA, 0, 1, 2)) > x a b c 1 0 0 NA 2 1 NA 0 3 2 1 1 4 NA 2 2 > I know it is easy to replace all the NA's with zeroes. > x[is.na(x)] <- 0 > x a b c 1 0 0 0 2 1 0 0 3 2 1

Dear R-helpers, I'd like to understand how to test the statistical significance of a random effect in gamm. I am using gamm because I want to test a model with an AR(1) error structure, and it is my understanding neither gam nor gamm4 will do the latter. The data set includes nine short interrupted time series (single case designs in education, sometimes called N-of-1 trials in medicine)

Hi, I have a small dataset that I'm fitting a segmented regression using nls on. I get a step below minimum factor error, which I presume is because residual sum of square is still "not small enough" when steps in the parameter space is already below specified/default value. However, when I look at the trace, the convergence seems to have been reached. I initially thought I might

Hello R-helpers I want to put a color scale in a plot: I've got an xy plot where the values of the response (z=f(x,y)) is symbolically given by colors (like heat or rainbow color scale) I would like to put such a scale with apprpriates labels in the plot, so as to facilitate the interpretation (like in a finite elements result plot) How is taht possible? Thanks Anne [[alternative