similar to: R 1.80 version

Hi All, I have already used gcc to creat a library, for example, it is called "test.a", which includes many c objects. How to link my main "a.c" function to this library when I used Rcmd SHLIB to complier my main c program? I tried "Rcmd SHLIB a.c test.a". It does not work this way. Does anyone have idea about this? Thanks, Xiaoping

I have written a wrapper for heatmap.2 called heatmap.w.row.and.col.clust which auto-generates breaks using breaks<-round((c(seq(from=(-20 * stddev), to=(20 * stddev))))/20, digits = 2) #(stddev in this case = 2.5) This has always worked well in the past but now I am getting an error that non-finite breaks are being generated. Drilling down, it seems that my wrapper is generating finite

I'm a programmer in a biology lab who is starting to use R to automate some of our statistical analysis of growth rate determination. But I'm running into some problems as I re-code. 1) Hypotheses concerning Slope similarity/difference: I'm using R's anova(lm()) methods to analyse a model which looks like this: growth.metric ~ time * test.tube I understand that

Dear useRs, `multcomp' version 0.991-1 will be shortly available from CRAN near you. Nearly all functionality contained in the package has been re-implemented from scratch. The focus of the package has been extended to general linear hypotheses in arbitrary parametric models and the most important function to check out is `glht()'. Multiple comparison of means procedures (for example

Dear useRs, `multcomp' version 0.991-1 will be shortly available from CRAN near you. Nearly all functionality contained in the package has been re-implemented from scratch. The focus of the package has been extended to general linear hypotheses in arbitrary parametric models and the most important function to check out is `glht()'. Multiple comparison of means procedures (for example

Hello, I am using R2.5 under Windows. Looks like the following statement vars <- (obj$sigma^2)*vw in getVarCov.gls method (nlme package) needs to be replaced with: vars <- (obj$sigma*vw)^2 With best regards Andrzej Galecki Douglas Bates wrote: >I'm not sure when the getVarCov.gls method was written or by whom. To >tell the truth I'm not really sure what

Hi, I am experimenting with using glht() from multcomp package together with coxph(), and glad to find that glht() can work on coph object, for example: > (fit<-coxph(Surv(stop, status>0)~treatment,bladder1)) coxph(formula = Surv(stop, status > 0) ~ treatment, data = bladder1) coef exp(coef) se(coef) z p treatmentpyridoxine -0.063 0.939 0.161

Dear all, For analysis of a longitudinal data set with fixed measurement in time I built a gls model (nlme). For testing hypotheses in this model I used the linear.hypothesis function from the car package. A check with the results obtained in SAS proc MIXED with a repeated statement revealed an inconsistency in the results. The problem can be that the linear.hypothesis function (1) only gives the

Hi all, I am fitting a linear mixed model with lme4 in R. The model has a single factor (des_days) with 4 levels (-1,1,14,48), and I am using random intercept and slopes. Fixed effects: data ~ des_days Value Std.Error DF t-value p-value (Intercept) 0.8274313 0.007937938 962 104.23757 0.0000 des_days1 -0.0026322 0.007443294 962 -0.35363 0.7237 des_days14 -0.0011319

Hello all! I would like to start off by saying that I am still really new to the vast world of R so please excuse my very limited vocabulary in the program. I have collected data from monkey videos and would like to setup some model(s) in R that can help with my hypotheses. I am having trouble figuring out which statistical tests/models to use for my two hypotheses. #1: Comparing the presence

[Env: Tested under Win Xp, R 2.9.2 and R 2.10.1; sessionInfo() at end] I've run into a problem with the heplot3d() function in my heplots package which causes the R GUI to crash ('R for Windows GUI encountered a problem and needs to close...'). I think the problem comes from an rgl call, but, I can't get a traceback or other information because my R session crashes. I've

Hi there, as part of a 2 x 2 contingency table analysis I would like to estimate conditional probabilities (success rates) in a Bernoulli experiment. In particular I want to test a null hypothesis p <= p0 versus the alternative hypothesis p > p0. As far as I understand the subject, there are UMPU tests for these types of hypotheses. Now I know about R's "binom.test" but the

No, since I'd like to test null: p <= p0 alternative: p > p0. and my understanding is that binom.test tests null: p = p0 (can only be a "simple" null hypothesis according to help(binom.test)) alternative: p > p0 (or p < p0 or p != p0). Thanks, Mirko. > -----Urspr?ngliche Nachricht----- > Von: Douglas Bates [mailto:bates at stat.wisc.edu] >

Dear all, I was wondering if you could give me any suggestions/help on the following issue. So I carried out the analysis of my data using generalized linear model (glm). After that, to check for multiple comparisons, I applied the glht function from the multcomp package in R. The output, however, gave me a warning (please see below). So my question is whether this warning is smth that I should

Dear All I have run the following GLM binominal model on a dataset composed by the following variables: TRAN_DURING_CAMP_FLG enviados bono_recibido 0 1 benchmark 0 1 benchmark 0 1 benchmark 0 1 benchmark 0 1 benchmark 0 1

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

Hello, I have problems interpreting the results of a Friedman test. It seems to me that the p-value resulting from a Friedman test and with it the "significance" has to be interpreted in another way than the p-value resulting from e.g. ANOVA? Let me describe the problem with some detail: I'm testing a lot of different hypotheses in my observer study and only for some the premises

Hello again R help. This is a simple question perhaps(laughing as I type) with a simple answer. Is the multcomp package appropriate for using with a lme built under the nlme package. I know i can get it to work (i.e report p-values for a tukey test) but i am not unsure if this is appropriate for a linear mixed effects model. Below is the code i used lmm1 <- lme(asinh(density_recruit) ~

Hi, I have data with binomial response variable (survival) and 2 categorical independent variables (site and treatment) (see below).? I have run a binomial GLM and found that both IVs and the interaction are significant.? Now I want to do a post-hoc test for all pairwise comparisons to see which treatment groups differ.? I tried the glht function in the multcomp package, but I get surprising