similar to: Compare a group of line slopes

Hi I have (8 measures * 96 groups) = 768 datasets for which I did linear regressions using lm(). Now I want to compare the slopes for each of the 8 measures in each of the 96 groups. As I understand , I can not use > anova(lm1, ..., lm8) as the lm1 ... lm8 are based on different datasets. I also read in previous discussions in this list, that I can see if the slope +- stddev(slope)

Hi, I've made a research about how to compare two regression line slopes (of y versus x for 2 groups, "group" being a factor ) using R. I knew the method based on the following statement : t = (b1 - b2) / sb1,b2 where b1 and b2 are the two slope coefficients and sb1,b2 the pooled standard error of the slope (b) which can be calculated in R this way: > df1 <-

Hello, Is there any formula or way to compare slopes of different functions? If we fit 2 functions in our data, and we have 2 slope parameters, how can we compare these slopes? Plotting y=5x and y=exp(5x) in which slope is equal to 5 in both of them.. doesn't seem that it makes sense to compare them. Maybe what I ask is basic statistics.. but you may be aware of some formula that could

Hi, I read almost all I found in prior R-Help list about How to compare two regression line slopes. So, I made a didactic example to illustrate a solution cited by Ben Bolker: =============================================== Subject: Re: [R] How to compare two regression line slopes From: Ben Bolker (bol... at ufl.edu) Date: Jan 27, 2009 1:52:20 pm List: org.r-project.r-help

Hi All, I am confused with SMATR's test for common slope. My null hypothesis here is that all slopes are parallel (common slopes?), right? So if I get a p value < 0.05 means that we can have confidence to reject it? That slopes are different? Or the other way around? it means that we have statistical confidence that the slopes are parallel? thanks -- Eugenio Larios PhD Student University

Hello, We've got a dataset with several variables, one of which we're using to split the data into 3 smaller subsets. (as the variable takes 1 of 3 possible values). There are several more variables too, many of which we're using to fit regression models using lm. So I have 3 models fitted (one for each subset of course), each having slope estimates for the predictor variables.

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

Hi, I'm trying to get confidence intervals to slopes from a linear model and I can't figure out how to get at them. As a cut 'n' paste example: ################# # dummy dataset - regression data for 3 treatments, each treatment with different (normal) variance x <- rep(1:10, length=30) y <- 10 - (rep(c(0.2,0.5,0.8), each=10)*x)+c(rnorm(10, sd=0.1), rnorm(10,

Hello everyone, I would like to test two regression slopes:do they differ significantly?The data and commands I've used so far: x<-8.5:32.5 #Vektor x y<-c( NA , NA , 5.67 , 6.53 , 6.83, 7.41 , 7.93 , 8.5 , 8.86, 9.46 , 9.82 , 10 ,10.35 , 10.7 ,11.03 ,11.37 ,11.61 ,11.84, 12.12, 12.39 ,12.67 ,12.96, 13.28 ,13.47, 13.65) #Vektor y (regression<-lm(y~x)) summary(regression)

Not really r-specific: Z = (b1 - b2) / SQRT ( SEb1^2 + SEb2^2) -------Original Message------- From: Gijsbert Stoet <stoet at volition.wustl.edu> Sent: 07/20/03 09:51 PM To: r-help at stat.math.ethz.ch Subject: [R] how to test whether two slopes are sign. different? > > Hi, suppose I do want to test whether the slopes (e.g. determined with lsfit) of two different population are

Hello, I have a linear model with different slopes for different treatment groups. I need to pairwise compare the different slope estimates for the different treatment groups. Is there a package that does pairwise comparisons of slope coefficients, making the appropriate adjustments in the P values? Thanks, John. -- ========================================== John Janmaat Department of

Hello, I've written a simple (although probably overly roundabout) function to test whether two regression slope coefficients from two linear models on independent data sets are significantly different. I'm a bit concerned, because when I test it on simulated data with different sample sizes and variances, the function seems to be extremely sensitive both of these. I am wondering if

Would like to try my luck to see if I can catch your eyes. I was trying to do some contrasts within ANOVA. I searched the archive and found a clue posted by Steffen Katzner ( http://tolstoy.newcastle.edu.au/R/help/06/01/19385.html) I have three levels for a factor names "StdLot" and want to make three comparisons, 1 vs 2, 1 vs 3 and 2 vs 3. First,

I have data from a design in which items are completely nested within subjects. Subject is the only second-level factor, but I have multiple first-level factors (IVs). Say there are 2 such independent variables that I am interested in. What is the proper syntax to fit a mixed-effects model with a by-subject random intercept, and by-subject random slopes for both the 2 IVs? I can

Hi I want to extract the random slopes from a lmer (I am doing a random regression), but are the answers obtained from ranef or coef? My model is: mod1<-lmer(B~ A +(A|bird), family=quasibinomial) And I want to obtain a slope for each individual bird but am not sure which output I need and can't find the answer anywhere. Thanks Sam Dr Samantha Patrick EU INTERREG Post Doc Davy 618

Dear list, I want to show all the objects starting with "d" and ending with a digit. How do I specify these conditions in the pattern argument I can do one condition but not two ls(pattern='^d') ls(pattern='[[:digit:]]$') are working. But, ls(pattern='^d'&'[[:digit:]]$') is not working. Appreciate any comment. Jun Shen [[alternative HTML

Can you provide the summary(m2) results? > -----Original Message----- > From: Simon Pickett [mailto:S.Pickett at exeter.ac.uk] > Sent: Wednesday, August 16, 2006 7:14 AM > To: Doran, Harold > Cc: r-help at stat.math.ethz.ch > Subject: [SPAM] - RE: [R] REML with random slopes and random > intercepts giving strange results - Bayesian Filter detected spam > > Hi again,

Hi, everyone, I was trying to sort a data frame by two columns, one increasing, the other decreasing and got an error. "Error in FUN(left) : invalid argument to unary operator", The command is "BA[order(BA[1],-BA[2]),]". BA is the data frame. It was working if I used increasing on both columns. Why the decreasing symbol "-" is not working here? Thanks. -- Jun

Hi, How would one go about determining if the slope terms from an analysis of covariance model are different from eachother? Based on the example from MASS: library(MASS) # parallel slope model l.para <- lm(Temp ~ Gas + Insul, data=whiteside) # multiple slope model l.mult <- lm(Temp ~ Insul/Gas -1, data=whiteside) # compare nested models: anova(l.para, l.mult) Analysis of Variance

Dear all, I am trying to estimate the prediction from a fixed effects model and their confidence intervals as well. Though I do not want to include in the prediction and at the confidence intervals the intercept. For that reason I used the argument incl.non.slopes=FALSE. But either if it is TRUE or FALSE it does not have any difference and also the system does not provide any warning. I really