search for: pallier

...comes to my mind: I could save all the graphics in png format (using R <script.R), and when it is finished, call ImageMagick's display to show all the png (or use any other diaporama system). However, I find it a dirty hack. Is there a simpler and cleaner way to achieve this? Christophe Pallier www.pallier.org

...ch the person managing the list at > r-help-owner at stat.math.ethz.ch > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of R-help digest..." > > > Today's Topics: > > 1. Re: "Statistiques avec R" (pallier) > 2. Re: graphic device MetaPost (Prof Brian Ripley) > 3. Re: graphic device MetaPost (Gabor Grothendieck) > 4. Re: graphic device MetaPost (Itay Furman) > 5. applying data generating function (Fred J.) > 6. missing values (Grace Conlon) > 7. drawing filled coun...

Hi, I am trying to write a function with the following codes and I would like it to return the values for "alpha beta para parab " seperately. Then I would like to use this funstion for "variable" with factor "a" and "b". But the result turns out to be a matrix with element like "Numeric,2" ... I guess they are just the values for

Dear All, I am doing some exercise in statistics textbook on comparison of two experimental means. Is it possible to use t.test() do t-test when I have only two means, sample size, two standard deviations ? (no raw data). Thanks. Pramote

...e is categorical. But in all the examples provided, the units of analysis contribute only one measurement. Should I include the subject and sentences factors in the formula? But then they would be treated as fixed-factors in the analysis, would they not? Any suggestion is welcome. Christophe Pallier www.pallier.org

Is there a conventional way to test for nested factors? I.e., if 'a' and 'b' are lists of same-length factors, does each level specified by 'a' correspond to exactly one level specified by 'b'? The function below seems to suffice, but I'd be happy to know of a more succinct solution, if it already exists. Thanks, Tim. --- "%nested.in%" <-

I'm trying to learn how to do a repeated measures ANOVA in R using lme(). A data set that comes from the book Design and Analysis has the following structure: Measurements (DV) were taken on 8 subjects (SUB) with two experimental levels (GROUP) at four times (TRIAL). In SAS, I use the code: PROC MIXED DATA=[data set below]; CLASS sub group trial; MODEL dv = group trial group*trial;

I have 2 files containing data analysed by 2 different methods. I would like to find out which genes appear in both analyses. Can someone show me how to do this? _________________________________________________________________ [[trailing spam removed]] [[alternative HTML version deleted]]

...o of aggregate.formula: my.aggregate.formula = function(formula,FUN=mean) { { d=model.frame(formula) factor.names=lapply(names(d)[sapply(d,is.factor)],as.name) factor.list=lapply(factor.names,eval) names(factor.list)=factor.names aggregate(d[1],factor.list,FUN) } Christophe Pallier http://www.pallier.org --------------- HEre is the code for aggregate.data.frame that recovers the name sof the factors: my.aggregate.data.frame <- function (x, by, FUN, ...) { if (!is.data.frame(x)) { x <- as.data.frame(x) } if (!is.list(by)) stop(...

...;. Currently, it does not mention the use of factors as indexes (It is true that factors are only defined later). At least, this manual could tell the reader to read the help about Extract, or type '?"["' (which did not occur to me before writing this message...). Christophe Pallier

Dear all, I would like to take out the values from one vector that are equal to the values in another vector. Example: a <- c(1,2,3,4,5,6,7,8,9) b <- c(3,10,20,5,6) b_noRepeats = c(10,20) So I would like to have the vector b without the same values as vector a. Kind regards, João Fadista [[alternative HTML version deleted]]

I am using R to make two-way ANOVA on a number of variables using g <- aov(var ~ fact1*fact2) where var is a matrix containing the variables. However the outcome seem to be dependent on the order of fact1 and fact2 (i.e. fact2*fact1) gives a slightly (factor of 1.5) different result. Any ideas why this is? Thanks for any help Anders

Hello All, I have a more statistical question, and how this is implemented in R. The problem is the following: We have 2 different solutions (samples), which are filtered and then the concentration of the filtrate is measured. We want to evaluate how the filter proces and the concentration measurement influences the detection of the difference of the two solutions and which step has which

Dear R-help, I've encounter what seems to me a strange problem with "names<-". Suppose I define the function: fun <- function(x, f) { m <- tapply(x, f, mean) ans <- x - m[match(f, unique(f))] names(ans) <- names(x) ans } which subtract out the means of `x' grouped by `f' (which is the same as, e.g., resid(lm(x~f)) if `f' is a factor).

Hello Coming from matlab background, I could use some help on this one. I need to generate a data set based on this equation X(t) = 3.8x(t-1) (1-x(t-1)) + e(t), where e(t) is a N(0,0,001) random variable I need say 100 values. How do I do this? Thanks

If I combine elements into a list b <- c(22.4, 12.2, 10.9, 8.5, 9.2) my.c <- sample.int(round(2*mean(b)), 5) my.list <- list(b, my.c) the names of the elements seems to get lost in the process: > str(my.list) List of 2 $ : num [1:5] 22.4 12.2 10.9 8.5 9.2 $ : int [1:5] 11 8 6 9 20 If I explicitly name the elements at list-creation, I get what I want: my.list <- list(b=b,

Hi, I have around 300 files with the x, y and z coordinates of a DEM that I would like to join in a single data frame object. I know how to automate the import of the 300 data frames. in Bash ls > names in R names<-scan(names...) With rep() and data.frame() construct a series of read.table() commands, which you write to a file and execute them with source(). I do not know however how

I'm new to R and am trying to perform a simple, yet problematic task. I have two variables for which I would like to measure the correlation and plot versus each other. However, I have ~30 million data points measurements of each variable. I can read this into R from file and produce a plot with plot(x0, x1) but as you would expect, its not pretty to look at and produces a postscript

Hi, I am using R to perform wilcox.test and wondering if you know how the p-value in wilcox.test is calculated? Thank you for your help:-) Regards, Latha

Hi all, So I think I have seen some similar questions to mine when I searched the archives, but have not seen any concrete answers and was wondering if any one could help. I have been trying to use R's aov() function to analyze my data. I have a 3 x 4 x 2 repeated measures design. All of the IVs are within subjects. I do also have missing values (unequal N), as I have to remove