similar to: Computing effect sizes based on mixed models

Dear All I would like to produce interaction boxplots and this seems to work: par(mfrow=c(2,2)) A=sample(rnorm(50,50,10)) B=sample(rnorm(50,100,10)) Test=merge(A,B,by=0)#by=0 where 0 is the row.names TreatA=(gl(2,50,100,labels=c("High","Low"))) TreatB=rep(gl(2,25,50,labels=c("High","Low")),2) Newdata=data.frame(TreatA,TreatB,Test)

Dear R-list, I have a data set from 20 pigs, each of which is tested at crossed 9 doses (logdose -4:4) and 3 skin treatment substances when exposed to a standard polluted environment. So there are 27 patches on each pig. The response is irritation=yes/no. I want to determine "equally effective 50% doses" (similar to old LD50), and to test the treatments against each other. I am looking

Gday, This is a repost since I only had one direct reply and I remain mystified- This may be stupidity on my part but it may not be so simple. In brief, my problem is I'm not sure how to extract parameter values/effect sizes from a nonlinear regression model with a significant interaction term. My data sets are dose response curves (force and dose) for muscle that also have two

G'day, My problem is I'm not sure how to extract effect sizes from a nonlinear regression model with a significant interaction term. My data sets are multiple measurements of force response to an agonist with two superimposed treatments each having two levels. This is very similar to the Ludbrook example in Venables and Ripley. The experiment is that a muscle is exposed to an agonist

Hi all. While working on something I got bitten by a mismatch between .o files due to a changed struct. This is easily prevented. Building proper dependency information into the Makefiles would be major surgery and require ongoing maintenance (although some of that could be automated by parsing the .depend files generated on OpenBSD). This patch is basically the brute-force approach: it will

?load Read this carefully. Pay attention to its instructions re: overwriting existing objects. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Tue, Jan 16, 2018 at 12:43 AM, Steven Yen <styen at ntu.edu.tw> wrote: >

Hello: I am a rookie in using R. When I used the unit root test in "CADFtest", I got the different t-test statistics between using criterion="BIC" and no using criterion. But when I checked the result with eviews, I find out that no using criterion is correct. Why after using criterion="BIC", I got the different result? Paul > data(Canada) > ADFt

I ran two separate hours-long projects. Results of each were saved to two separate .RData files. Content of each includes, among others, the following: ?????????????????? me??? se????? t???? p sig pc21.age??????? 0.640 0.219? 2.918 0.004 *** pc21.agesq????? 0.000 0.000??? NaN?? NaN pc21.inc??????? 0.903 0.103? 8.752 0.000 *** pc21.incsq????? 0.000 0.000??? NaN?? NaN pc21.sei10????? 0.451 0.145?

Having spent the last few weeks trying to decipher R, I feel I may finally be getting somewhere, but i'M still in need of some advice and all my tutors seem to be on holiday! Basically a bit of background, I have data collected on a population of Lizards which includes age,sex, and body condition. I collected data myself this year and I have data previously collected from 1999, 2002 and

Hi all, I have a question about linear model with interaction: I created a data frame df like this: >df V1 V2 V3 V4 V5 1 6.414094 c t a g 2 6.117286 t a g t 3 5.756922 a g t g 4 6.090402 g t g t ... which holds the response in the first column and letters (a,c,g,t) in the other columns. I am interested to see if there are interactions between the neigbouring letters

Hi all, I have a question about linear model with interaction: I created a data frame df like this: >df V1 V2 V3 V4 V5 1 6.414094 c t a g 2 6.117286 t a g t 3 5.756922 a g t g 4 6.090402 g t g t ... which holds the response in the first column and letters (a,c,g,t) in the other columns. I am interested to see if there are interactions between the neigbouring letters

Huh. I may by completely wrong but you cannot do such "merging". .RData files are AFAIK places where all objects from given session are stored. However you could load each .RData file and save/export result (one object). BTW, what do you mean exactly by "combine/consolidate"? And finally, post your questions in plain text not html, otherwise they can be mangled. Cheers

Hello. I'm a new user of R, and I have a question regarding the use of aov and lm-functions. I'm doing a fractional factorial experiment at our production site, and I need to familiarize myself with the analysis before I conduct the experiment. I've been working my way through the examples provided at http://www.itl.nist.gov/div898/handbook/pri/section4/pri472.htm

On 16/01/2018 6:33 AM, Steven Yen wrote: > Hi all, > This is great. Again, here is what I need. I run two separate jobs (a.R > and b.R) with results (say regression outputs) going to a.RData and > b.RData. I like to put all results in one place (where I can retrieve > them in one place, ab.RData). The following codes do it (I am not sure > if line 2 is needed but I am happy).

Understood. In my case, a.RData and b.RData contain identical variables/data, plus simulation outputs from separate runs. The codes deliver what I need. Good to know the three lines work. Thank you. On 1/16/2018 8:06 PM, Duncan Murdoch wrote: > On 16/01/2018 6:33 AM, Steven Yen wrote: >> Hi all, >> This is great. Again, here is what I need. I run two separate jobs (a.R >>

Hi all, This is great. Again, here is what I need. I run two separate jobs (a.R and b.R) with results (say regression outputs) going to a.RData and b.RData. I like to put all results in one place (where I can retrieve them in one place, ab.RData). The following codes do it (I am not sure if line 2 is needed but I am happy). Thank you all. load("a.RData")

Hello, While testing recent snapshots (20131023 and 20131024) I encountered a problem with ProxyCommand. The regression tests all passed, but the use of ProxyCommand's in my ~/.ssh/config resulted in name resolution errors; even if CanonicalizeHostname was explicitly set to "no." The patch included inline below fixed the issue: Index: ssh.c

Hello everyone I'm trying to use Scheffe's complete cubic model (mixture design). In the bibliographies, they indicate that the term is of the type: A * B * (A-B). But I see that trying to adjust the three cubic terms results in singularities. I know this implies not having the inverse matrix: solve (t (X)% *% X) does not exist. The bibliographies show all three cubic terms. So my

Since the fact that gfortran performance has improved over the major releases, I decided to benchmark the current releases on a MacPro with the Polyhedron 2005 benchmarks using -ffast-math -funroll-loops -msse3 -O3. The results are... gcc release gcc 4.2.4 gcc 4.3.3 gcc 4.4-pre gcc 4.3.3/ gcc 4.4-pre/

Hi, I just installed Sipp 3.3?on CentOS 6.3 and all of the calls Sipp is generating are failing. I am trying to run Sipp on the same machine as Asterisk PBX using the ./sipp -sn uac 192.168.1.115 command. SIpp output: ----------------------------- Statistics Screen ------- [1-9]: Change Screen -- ? Start Time???????????? | 2013-05-20?22:53:08:637?1369083188.637273??????????? ? Last Reset