similar to: contrasts (PR#7695)

Full_Name: David Clayton Version: 1.8.1 OS: Linux Submission from: (NULL) (131.111.126.242) qchisq behaves very strangely when ncp is passed as zero (forcing internal qnchisq to be called) when first argument is small. Eg > qchisq(1-1e-6, 1, ncp=0, lower.tail=TRUE) qchisq(1-1e-6, 1, ncp=0, lower.tail=TRUE) [1] 1024 while, if ncp is unspecified, > qchisq(1-1e-6, 1) qchisq(1-1e-6, 1)

Hello, I use a type III anova ("car" package) to analyse an unbalanced data design. I have two factors and I would have the effect of the interaction. I read that the result could be strongly influenced by the contrasts. I am really not an expert and I am not sure to understand indeed about what it is... Consequently, I failed to properly used the fit.contrast function (gregmisc

Full_Name: Hin-Tak Leung Version: R 2.2.0 OS: x86 linux Submission from: (NULL) (131.111.126.242) Rgnome depends libglade 0.x and libxml 1.x . They are no longer shipped with Redhat EL4. (I know they are still in fedora core 4). Just a wish-list.

Hi, I have got a curious issue with an R package which uses zlib, against the official binary here: http://cran.r-project.org/bin/linux/redhat/el4/i386/R-2.4.1-1.rh4AS.i386.rpm on a Redhat EL4 i686 system. The problem is that at the end of reading a gzip'ed file within my C code, gzgetc() returns -1 (no more to read or error) but gzeof() doesn't return true. Now one can probably jump to

I am using dgc.genetics to perform TDT analysis on SNP data from a cohort of trios. I now have a file with about 6008 variables. The first few variables related to the pedigree data such as the pedigree ID the person ID etc. Thereafter each variable is a specific locus or marker. The variables are named by a pattern such as "Genotype.nnnnn" with nnnnn corresponding to a number which

Hi, (somebody would probably yell at me for not checking 2.6.0rc, for which I can only apologize...) Our R package (snpMatrix in http://www-gene.cimr.cam.ac.uk/clayton/software/) is broken rather badly in 2.6.0 ; I have fixed most of it now so a new release is imminent; but I'd like to mention a few things, mostly to summarize my experience and hopefully the 'writing R extensions'

Hi, I'm having difficulty specifying contrasts for a within subjects factor with 3 levels. I can do it correctly for my factors with 2- levels, but i'm not getting the correct results for a 3-level factor. My design has 1 between subjects factor (gp) and 3 within subjects factors. Within factor "w" has 2 levels, within factor "x" has two factors, and within

My issue seems simple but I can't find any method to work... I've got a matrix with 4 columns and I want to extract the values from each row to make a new 2x2 matrix for each row. From these resultant matrices, I want to get the Eigen Values and output them from the function as usable numbers (not a text string). Here's the function... ### Ef = function(MAT) {for(i in

Dear experts: Is it possible to create a new function based on stats:::model.matrix.default so that an alternative factor coding is used when the function is called instead of the default factor coding? Basically, I'd like to reproduce the results in 'mat' below, without having to explicitly specify my desired factor coding (identity matrices) in the 'contrasts.arg'. dd

Thanks, Mark. Um, how's about from the commandline or how do I get, I guess we're still using X11, windows to load. Sorry, :<}}} Clayton >-----Original Message----- >From: centos-bounces at centos.org [mailto:centos-bounces at centos.org] On >Behalf Of Mark LaPierre >Sent: Sunday, December 14, 2014 5:50 PM >To: centos at centos.org; Mark LaPierre >Subject: Re:

Hi, sorry if this is explained somewhere but I didn't find anything. How can I use "Anova" from the car package to test a modell without between subject's factors? Suppose I have the following data mat.1 mat.2 mat.3 di ex 1 85 85 88 1 1 2 90 92 93 1 1 3 97 97 94 1 1 4 80 82 83 1 1 5 91 92 91 1 1 6 83 83

On Fri, Mar 10, 2023 at 10:26?AM Chris Clayton <chris2553 at googlemail.com> wrote: > > Hi. > > Is it likely that this fix will be sumbmitted to mainline during the ongoing 6.3 development cycle? > yes, it's already pushed to drm-misc-fixed, which then will go into the current devel cycle. I just don't know when it's the next time it will be pushed upwards, but it

On Mon, Feb 20, 2023 at 11:51 AM Chris Clayton <chris2553 at googlemail.com> wrote: > > > > On 20/02/2023 05:35, Ben Skeggs wrote: > > On Sun, 19 Feb 2023 at 04:55, Chris Clayton <chris2553 at googlemail.com> wrote: > >> > >> > >> > >> On 18/02/2023 15:19, Chris Clayton wrote: > >>> > >>> > >>> On

HI, I am trying to perform ANOVA with 2 factors: material (3), temperature(3). The interaction is significant. I tried something like, ( summary(avt, split=list("temp:mat"=list("15"=1, "70"=2, "125"=3))) is not correct). Thanks, av <- aov(time ~ mat*temp, data=dados) avt <- aov(time ~ temp/mat) summary(avt,

I got some R code that I don't understand. Question as comment in code //where is t comming from, what is phi inverse rAC <- function(name, n, d, theta){ #generic function for Archimedean copula simulation illegalpar <- switch(name, clayton = (theta < 0), gumbel = (theta < 1), frank = (theta < 0), BB9 = ((theta[1] < 1) | (theta[2] < 0)), GIG = ((theta[2] < 0) |

Dear all, In SAS (GLM and MIXED) estimable functions (linear functions of the parameters) can be specified in the ESTIMATE and CONTRAST statements. Has anyone written a similar "utility" for use in connection with lm? Thanks in advance S?ren H?jsgaard -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Personally, I am agnostic. I've just read thru Centos documentation that there is a big effort to remove all upstream personalities from Centos. Personally, I don't see why RH is doing this. I would think that it undermines RH. But I'm still new/old to all of this. It used to be the big argument was between Unix from Berkeley(4.? I think) and SysIII/V. Always always battles for turf.

I am using the 'glmmPQL function in the 'MASS' library to fit a mixed effects logistic regression model to simulated data. I am conducting a series of simulations, and with certain simulated datasets, estimation of the random effects logistic regression model unexpectedly terminates. I receive the following error message from R: Error in lme.formula(fixed=zz + arm.long,random=~1 |

Hi, I noticed there is an "-inplace" option on rsync web site's man page. So this seems to mean "copy then move" is the default behavior. Has this been the default all the time? I am running a pretty old version, see below, and I need to have the "copy then move" feature to ensure partial files are not accessed while they are being updated. [qt_prod@cflxqs55p

Hello, I have a few questions, regarding the smoothScatter function. I have a scatter plot with more than 500,000 data points for two samples. So, I am wanting to display the density in colors to convince people that my good correlation coefficient is not due to an "influential point effect" and plus, I also want to make my scatter plot look pretty. Anyway ... I have been able to