similar to: assign group letters to T/F matrix

Hello, following problem: i have a written my own function to draw some sophisticated graphic, which after manipulating data somewhere contains a plot comand function <- function(x) { ... plot(xy) } to tidy up my final graphs, it would be very handy to be able forward all low-level graphic parameters set in function to plot i.e. function(x, cex=0.5) calls plot (xy, cex=0.5) is there a way

dear r community... it loks like i won't be able to reach homogenity of variance for my dataset, so i end up with welch anova instead of regular anova. documentation on this test is rather scarce, so maybe someone here can enlighten me a bit: - do i understand that no two-way implementation of the welch anova has been developed yet? - is there a post-hoc test for welch anovas implemented in

Is there code in R to automatically assign letters to different groups that are found to be statistically significant using a post-hoc test following an ANOVA? For example, let's say I found that relative rates of tree growth were statistically significant between tree species using an ANOVA and then I ran a TukeyHSD post-hoc test to determine which species were statistically

Dear all, we have an open position for a research assistant for Biostatistics in a project focusing on patient (instance) specific statistical learning methods. Application deadline is September 27th 2013. Details can be found at http://www.hs-furtwangen.de/willkommen/aktuelles/termine-fristen/aktuelles-einzelansicht/1426-academic-employer-m-f-area-biostatistics.html Regards, Matthias --

Hi, All How does one remove relevant information from a regression output besides just the coefficients? I've been able to modify the example given under "help(by)" to give me some additional information, but not everything I need. If you adjust the call statement from what is listed by adding the summary statement like so: tmp <- by(warpbreaks, tension, function(x)

Dear R users, I'm trying to create images of 2 symmetrical matrices using heatmap.2 of the 'gplots' library. Both matrices have the same row and column names For the first matrix I'm using the default clustering and column denrogram options: heatmap.2(Matrix1, symm=T, breaks=20, col=rich.colors(19), main="matrix1", density.info="histogram",

Hello, I have some conflicting output from an aov summary and tukey contrasts with a mixed effects model I was hoping someone could clarify. I am comparing the abundance of a species across three willow stand types. Since I have 2 or 3 sites within a habitat I have included site as a random effect in the lme model. My confusion is that the F test given by aov(model) indicates there is no

Dear R-help, I am trying to express myself as best as I can here. If you also use Latex to edit math reports or other languages with similar editing method, you'll see what I'm talking about. My sincere appologies if my question is not clear enough to some extend, as also I'm not able to provide my code here because I don`t know which one I can use... When editing the title in R

Hello, I am interested in generating a compact letter display for the results of Tukey HSD tests that contain interaction effects. The 'cld' method in the 'multcomp' package seems only to work for main effects. Does such a thing exist already? Thank you for any thoughts, Josh -- Joshua Caplan, PhD Postdoctoral Associate Department of Ecology, Evolution, and Natural Resources

After some reorganization, we would like to announce the availability on CRAN of a new version (1.6) of our package "distr" as well as the availabilty on CRAN of the new packages "distrEx", "distrSim", "distrTEst", "RandVar". ----------------------------------------------------------------------------------------- Changes from 1.5 to 1.6 -[former

After some reorganization, we would like to announce the availability on CRAN of a new version (1.6) of our package "distr" as well as the availabilty on CRAN of the new packages "distrEx", "distrSim", "distrTEst", "RandVar". ----------------------------------------------------------------------------------------- Changes from 1.5 to 1.6 -[former

Hello, JIT compiler interferes with internal R completions: compiler::enableJIT(2) utils:::functionArgs("density", '') gives: utils:::functionArgs("density", '') Note: no visible global function definition for 'bw.nrd0' Note: no visible global function definition for 'bw.nrd' Note: no visible global function definition for 'bw.ucv'

It's annoying, for instance when I am typing anything in google, it will up something like this: h3l5 do02loadin2 and so on, I had recently downloaded "word processor" so I can write a paper for my class and I think I downloaded the wrong one and ended up downloading a virus... So I tried downloading AVG or anything that will get rid of the virus, which is not working I tried opening

This is my current work.Now i am trying to use a function to do the normal distribution simulation. rm(list=ls()) t <- u<- mann<- rep(0, 45) Nsimulation<-function(S1,S2,Sds,nSims) { set.seed(1) for (sim in 1:nSims) { matrix_t <-matrix(0,nrow=nSims,ncol=3) matrix_u<-matrix(0,nrow=nSims,ncol=3)

Hi. I was having trouble getting Greek letters and mathematical symbols in maple under wine, and a search turned up a few messages suggesting others have the same problem. In particular, Thomas Koenig wrote... http://www.winehq.com/hypermail/wine-users/2001/04/0028.html > I wrote: > > >I have tried several applications which use Greek letters, such as > >Maple V R5 and

Hi, thanks for this great work. Just a minor point - maybe you should mention the rm -rf $FAKE_ROOT command within your README file: It caught me unexpected and I was looking a couple of minutes for the package artifact before I actually read the last line. Regards, Burkhard Kohl Lead Engineer Content Management burkhard.kohl at intershop.de INTERSHOP? http://www.intershop.de INTERSHOP

Hi, Can anyone help me with the following. I have been using R for Monte Carlo simulations and got some results I couldn't explain. Therefor I performed following short test: -------------- mean.sds <- NULL sample.sizes <- 3:30 for(N in sample.sizes){ dum <- NULL for(I in 1:5000){ x <- rnorm(N,0,1) dum <- c(dum,sd(x)) } mean.sds<- c(mean.sds,mean(dum)) }

HI, I am trying to use switch () function to connect the three distribution (normal ,gamma with equal skewness and gamma with unequal skewness. But i am losing my ideas since i have sample sizes-(10,10),(10,25),(25,25),(25,50),(25,100),50,25),(50,100), (100,25),(100,100) standard deviation ratio- (1.00, 1.50, 2.00, 2.50, 3.00 and 3.50) distribution of gamma distribution with unequal skewness

Dear list, I have a vector of letter strings as follow: > LETTERS[c(1:7,9,15,18:25)]  [1] "A" "B" "C" "D" "E" "F" "G" "I" "O" "R" "S" "T" "U" "V" "W" "X" "Y"   I need to find sequences of length 11, which are made from  all

What is the easiest way to change within vector of strings each letter after a space into a capital letter? E.g.: c( "this is an element of the vector of strings", "second element" ) becomes: c( "This Is An Element Of The Vector Of Strings", "Second Element" ) My reason to try to do this is to get more readable abbreviations. (A suggestion would be to