similar to: is this an ANOVA ?

Hi, there I'm trying to get the value of the Mean Square from the ANOVA model summary that comes from specifying the error term, and am wondering if one can actually do this ( I know it's possible when using anova(lm) objects and the like, but I'm having a tough time with it under this framework). There does appear to be some indexing in the output of this type, but perhaps not

I'm just getting to grips with using ode function and have used the examples and vignettes to produce a small model of a one-pool, michaelis-menten, enzyme kinetic reaction. The rate of flux of substrate into pool A is constant (fluxoa) however the rate of flux out of pool A is controlled by the HMM equation (v = Vmax/ ( 1 + (Km / Concentration A )) ). This function works fine and

Hello, In a book (David W. Stockburger, "Multivariate Statistics: Concepts, Models, and Applications", chapter 12 "Contrasts, Special and Otherwise", available online at http://www.psychstat.smsu.edu/multibook) I've found some examples of doing analysis of variance on a contrast basis. I attach my solution (in R, the book uses SPSS) to this problem. Am I computing the

I am a Professor of Statistics at Indira Gandhi Krishi Vishwavidyalaya, Raipur, India. While teaching in class about analysis of variance using R, I was doing a one-way analysis for the two data-sets given below in the R-class. I got a typical error in "emmeans" package, please help: Data-set-1: -------------- Medley and Clements (1998) investigated the impact of zinc contamination (and

Hi, I'm analyzing experimental results where two different events ("T1" and "T2") can occur or not during an experiment. I made my experiments with one factor ("Substrate") with two levels ("Sand" and "Clay"). I would like to know wether or not "Substrate" affects the occurrence probability of the two events. Moreover, for each

On 07/01/18 02:19, Akhilesh Singh wrote: > I am a Professor of Statistics at Indira Gandhi Krishi Vishwavidyalaya, > Raipur, India. While teaching in class about analysis of variance using R, > I was doing a one-way analysis for the two data-sets given below in the > R-class. I got a typical error in "emmeans" package, please help: > > Data-set-1: > --------------

Thanks for your kind reply. Problem is solved. However, it's "confidence interval / treatment comparison plot" is not taking main title. And the fonts of axes labels can not be changed using 'cex' parameter. I will appreciate if you could help in this matter too. Dr. A. K. Singh On 09-Jan-2018 8:18 PM, "Sal Mangiafico" <salvatore.s.mangiafico at gmail.com>

Hi all; I'm trying to display a 2 panel plot for the Puromycin data from R with 2 different non-linear models fitted to each group. The problem is that as far as I know panel.number doesn't work in the latest version of R. Can anyone give a hint how to solve this? Here is the code that I used before and now doesn't work xyplot(rate ~conc| state,Puromycin,

Hi, my question is not exactly about R... What I am looking for are hints and directions on suitable methods (available in R or elsewhere) to solve a grouping (or pattern recognition) problem of environmental features in an environmental gradient as described below. Given environmental sampling data set (Depth, Presence of sand, Presence of boulders, Presence of clay). 1 1 1 0 1 1 0 0 1 1

UseRs, I am working on a dataset (see small example below) where individuals were followed on a specific date-time combo and multiple repeated measurements were taken (e.g., height in meters, behavior class in 2 letter code). Observation numbers varied between individual (ranging from 1 observation for each date-time combo to >50) I am trying to summarize the data into 1 row per

1) deleting path/to/(Something) Word Word Anotherword--Word Word.pdf IO error encountered -- skipping file deletion I think either the -- is causing it to fail, or more likely, the next filename, which it didn't bother to print out, which has an ampersand in the name, causes it to fail. I am using SSH as a substrate, of course, and I think you guys failed to escape things, so the ampersand

I compiled R 2.1.0 under Windows XP SP2 as a preliminary to rebuilding a custom package for use with R 2.1.0. The compile completed successfully, and I was able to run demo(graphics) successfully. But make check and make check-recommended fail. > version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 2 minor 1.0 year 2005

Dear R users, Could somebody please help me to find a way of comparing nonlinear, non-nested models in R, where the number of parameters is not necessarily different? Here is a sample (growth rates, y, as a function of internal substrate concentration, x): x <- c(0.52, 1.21, 1.45, 1.64, 1.89, 2.14, 2.47, 3.20, 4.47, 5.31, 6.48) y <- c(0.00, 0.35, 0.41, 0.49, 0.58, 0.61, 0.71, 0.83, 0.98,

Hi! I wish to know if there is any R forum which is meant only for Statistics? I mean where we can clarify our statistics doubts and seek knowledge. I know there are lot many books and internet sites, but 'R forum' has altogether different standard and very high level and one can learn a lot from them. Regards Vincy [[alternative HTML version deleted]]

Hello, I am performing a BACI analysis with ANOVA using the following glm: fit1<-glm(log(Cucs_m+1)~(BA*Otter)+BA+Otter+ID+Primary, data=b1) The summary(aov(fit1)) shows significance in the interaction; however, now I would like to determine what combinations of BA and Otter are significantly different (each factor has two levels). ID and PRIMARY substrates are categorical and included in

Dear R_Users, I am confronted to a technical issue that I wouldn't even know how to start to handle and where to look at. It is about data manipulation in a dataframe and value substitution Imagine I have a matrix describing a set of chemical equations. Colnames are Components names and Rownames are Reaction names Each cell in that matrix corresponds to the stoechiometric coefficient to the

Dear all, I've taken a subset of data from a data frame using crb<-subset(all.raw, creek %in% c("CR") & year %in% c(2000,2001) & substrate %in% ("b")) this works fine, except that all of the original factor levels are maintained. This results in NA's for these empty levels when I try to do summaries based on factors using by(). Is there a simple way to drop

Hey All, I’m attempting to gather stats on usage of the “metal”, by which I mean the physical host’s hardware. I would like to know the CPU, IO, and network stats for the hardware. Normally I would get these from /proc/stat, /proc/diskstat, and ifconfig respectively. However, because of the segregation of VMs and Dom0’s special relationship as a paravirtualized I don’t know if I can trust this

Sorry this may well be defined as Off Topic. I apologize in advance. I am interested in performing what I think would be a probabilistic sensitivity simulation. I've done some crude ones before in excel but I'm wondering if R can help me do it more effectively? I have a set of theoretical variables for simplicity lets use (what I think) is an easier example: I have a peg and a hole

I'd like to model the relationship between m responses Y1, ..., Ym and a single set of predictor variables X1, ..., Xr. Each response is assumed to follow its own regression model, and the error terms in each model can be correlated. My understanding is that although lm() handles vector Y's on the left-hand side of the model formula, it really just fits m separate lm models. What should