similar to: power of 2 way ANOVA with interaction

I've tried without success to calculate interaction and main effects contrasts using R. I've found the functions C(), contrasts(), se.contrasts() and fit.contrasts() in package gmodels. Given the url for a small dataset and the two-way anova model below, I'd like to reproduce the results from appended SAS code. Thanks. --Dale. ## the dataset (from Montgomery) twoway <-

Hi, all There are 62 samples in my data and I tested 3 times for each one, then I want to use ICC(intraclass correlation) from irr package to test the consistency among the tests. *combatexpdata_p[1:62] is the first text results and combatexpdata_p[63:124] * is the second one and *combatexpdata_p[125:186]* is the third. Here is the result:

I sent a message a couple days ago about doing calculations for power of the ANOVA. Several people got back to me very quickly which I really appreciated. I'm working now on a similar problem, but instead of a balanced ANOVA, I have an unbalanced one. The first part of the question was: You assume that the within-population standard deviations all equal 9. You set the Type 1 error rate at á

Hi all, I've been running loglinear models for three-way tables: one of the variables having three levels, and the other two having two levels each. An example looks like below: > yes.no <- c("Yes","No") > switch <- c("On","Off") > att <- c("BB","AA","CC") > L <- gl(2,1,12,yes.no) > T <-

I have been attempting some basic power calculations using R and I am not getting the results I expect. I had a homework assignment in SAS, but I want to learn R as well, so I was attempting to reproduce my result. (No one else in the class is doing R, so there's no need to obsfucate the answer, the SAS code is what I get my grade for.) The code I am using is: # You assume that the

Hi... In Stata, there is the ability to display scatter plots with data points at the same (x,y) location, using the 'jitter' command of the twoway scatter stata command. Anyone know of a way that I can do the equivalent thing in R? For non-Stata readers, if jitter is enabled in stata, and n-points would be at the same (x,y) location, the points are offset a bit (according to

I am a new R user but a long time SAS user. I searched for a response to this question but no luck, so forgive me if this topic has been covered before. I am running a TukeyHSD post hoc test after running an ANOVA. I get the results of all pairwise comparisons, no problem. However, the output table is a little "busy", and I'd like to make the output easier to read. Specifically, I

Hi R Users, I have a dataframe like this: id x rho A 1 0.1 B 20 0.5 C 2 0.9 ... I want to do a scatter plot of "x" versus "rho" but for each point on the scatter plot I want the corresponding entry for "id" instead of points. In STATA I can do so by twoway (scatter x rho, mlabel(id)) How can I do the same in R? I am sure there is some simple way to do

Hello, I want to know if there are some functions or packages to solve differential and integral equation using R. Thanks. Shao chunxuan. [[alternative HTML version deleted]]

Has anyone ever setup Asterisk to pass Feature Group D access while using a CIC code for outbound calls? If so can you please email the configuration you have done? I have tried to get this up and running but with no luck. I have also contacted support and I cant seem to get this going. Thanks in Advance, Jason Miller

Hi, In SPlus there is a function called "twoway" for median polishing gridded data. Is there an equivalent function in R? I have been searching for it in R help without much success. Your help is much appreciated. Cheers Sorn [[alternative HTML version deleted]]

Hi, everyone: I ask for help about translating a stata program into R. The program perform cross validation as it stated. #1. Randomly divide the data set into 10 sets of equal size, ensuring equal numbers of events in each set #2. Fit the model leaving out the 1st set #3. Apply the fitted model in (2) to the 1st set to obtain the predicted probability of a prostate cancer diagnosis. #4. Repeat

For various reasons, I spent part of my time today looking at sample size and power calculation tools (don't ask, don't tell...). This seems to be one area that R is incredibly weak in (well, nearly all stat packages, except perhaps specialized tools and SAS); sure, there are a number of functions in various packages: base, statmod, Hmisc Have I missed something? (I

I am interested in using "multcompLetters" after running "kruskalmc" but I'm a newbie and I'm not having luck figuring it out. I can run "kruskalmc" just fine, but after studying the documentation for "multcompletters" for a long time, I cannot figure out how to make it work. Any ideas? R input is below > >> #here is my data set >>

Hi helpeRs. When I do two-ways tables dimnames in columns are almost left-adjusted: <- factor(rbinom(30,1,.3),labels=c("first","second")) > b <- a > table(a,b) b a first second first 21 0 second 0 9 If I do a oneway table however, dimnames are right-adjusted: table(a) a first second 21 9 Actually I must admit I

Hi I'd really like to get a bar plot showing the means of my anova data. I have looked everywhere and can only seem to find instructions for 2 way anova's. I basically want to look at the mean condition of my subjects spilt by age, sex and year (as a factor rather than a continuous variable, hence Anova and not Ancova). and want to show it firstly as a bar graph with standard error. I

On Thu, 16 Feb 2023 17:30:37 +0000, Ulrich Sibiller <ulrich.sibiller at atos.net> wrote: > Martin Schwenke schrieb am 15.02.2023 23:23: > > OK, this part looks kind-of good. It would be interesting to know how > > long the entire failover process is taking. > > What exactly would you define as the begin and end of the failover? From "Takeover run

Hi, In Splus if I call the function crosstabs() the output is a contigency table; in each cell of the table is printed: N, N/RowTotal, N/ColTotal, N/Total. N is the number of observations in each cell. The same call to xtabs() in R will produce the contigency table but the only entry in each cell is N. How can I get the same relative frequencies that crosstabs() gives? Thanks, mike --

Hello everyone, I have a peoblem to create the twoways effect in the plm package. when i try to create the following dsn1<-plm(lnQ~lnC+lnL+lnM+lnE+eco+RD,data=newdata,effect="twoways",model="within") i have this error: Error in rep.int(c(1, numeric(n)), n - 1L) : negative length vectors are not allowed and to be honest i have no idea what does it mean!! can someone

Dear R users I get different F-statistic results for a "within" model, when using "time" or "twoways" effects in plm() [1] and when manually specifying the time control dummies [2]. [1] vignette("plm") [2] http://cran.r-project.org/doc/contrib/Farnsworth-EconometricsInR.pdf Two examples below: library("AER") data("Grunfeld", package =