search for: hypothesize

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

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

I'm coming from the AMOS world and am wondering if there is a simple way to do multiple hypothesis testing in the manner of BIC analyses in AMOS using the sem package in R. I've read the documentation, but don't see anything in there except for basic BIC scores. Perhaps someone has devised a simple way to compare the relative likelihood of all possible path-fittings within a

Hello everybody, I have some questions on ANOVA in general and on ANOVA in R particularly. I am not Statistician, therefore I would be very appreciated if you answer it in a simple way. 1. First of all, more general question. Standard anova() function for lm() or aov() models in R implements Type I sum of squares (sequential), which is not well suited for unbalanced ANOVA. Therefore it is better

Full_Name: Yves Rosseel Version: 2.2.1 OS: Submission from: (NULL) (157.193.116.152) The code below sketches a possible implementation of a function 'mlh.mlm' which I think would be a good complement to the 'anova.mlm' function in the stats package. It tests a single linear hypothesis of the form H_0: LBT'= 0 where B is the matrix of regression coefficients; L is a matrix

Dear All, I'm trying to model heteroscedasticity using a multilevel model. To do so, I make use of the nlme package and the weigths-parameter. Let's say that I hypothesize that the exam score of students (normexam) is influenced by their score on a standardized LR test (standLRT). Students are of course nested in "schools". These variables are contained in the Exam-data in the mlmRev package. library(nlme) library(mlmRev) lme(fixed = normexam ~ stand...

Dear Brian, As Duncan mentioned, the terms type-I, II, and III sums of squares originated in SAS. The type-II and III SSs computed by the Anova() function in the car package take a different computational approach than in SAS, but in almost all cases produce the same results. (I slightly regret using the "type-*" terminology for car::Anova() because of the lack of exact

I was wondering, whether there is a way to have statistical significance test for cluster agreement. I know that I can use classAgreement() function to get Rand index, which will give me some indication whether the clusters agree or not, but it would be interesting to have a formal test. Thanks.

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

Dear R Gurus: Is there a test for a single variance available in R, please? Thanks, Edna Bell

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 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

I have a three-way unbalanced ANOVA that I need to calculate (fixed effects plus interactions, no random effects). But word has it that aov() is good only for balanced designs. I have seen a number of different recommendations for working with unbalanced designs, but they seem to differ widely (car, nlme, lme4, etc.). So I would like to know what is the best or most usual way to go about working

Hi there, The following data is obtained from a long-term experiments. > mydata <- read.table(textConnection(" + y year Trt + 9.37 1993 A + 8.21 1995 A + 8.11 1999 A + 7.22 2007 A + 7.81 2010 A + 10.85 1993 B + 12.83 1995 B + 13.21 1999 B + 13.70 2007 B + 15.15 2010 B + 5.69 1993 C + 5.76 1995 C + 6.39 1999

I have a system with a binary response variable that was hypothesized to follow a simple logistic function. The relationship between the continuous independent variable and the logit is clearly not monotonic. I have two questions. 1) Can anyone recommend a reference that describes my modeling options in this case, and 2) what facilities does R have to deal with...

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

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

Could you please help me on the following ANCOVA issue? This is a part of my dataset: sampling dist h 1 wi 200 0.8687212 2 wi 200 0.8812909 3 wi 200 0.8267464 4 wi 0 0.8554508 5 wi 0 0.9506721 6 wi 0 0.8112781 7 wi 400 0.8687212 8 wi 400 0.8414646 9 wi 400 0.7601675 10 wi 900 0.6577048 11 wi 900