similar to: logistic regression by glm

I made an ANOVA with repeated mesures (aov(Mesure~Distance*Genre*Correct+Error(Sujet/(Distance*Genre*Correct)), data)) and I would like to extract the pvalues. The output is: ----------------------------------------------------------- Error: Sujet Df Sum Sq Mean Sq F value Pr(>F) Residuals 21 97.082 4.623 Error: Sujet:Distance Df Sum Sq Mean Sq F value Pr(>F) Distance

dear all, I am a new comer to R and statistic. Now I have a little confuse about the package pls. I have to use 5 components to form a model. There are strong relationship between some of the components, which leads to the changes of the sign of each coeficeince, of course this is unwanted when using the normal regression way. So I choose the way of PLS, which is good at solve this kind of

I tried to compare if cch() and coxph() can generate same result for same case cohort data Use the standard data in cch(): nwtco Since in cch contains the cohort size=4028, while ccoh.data size =1154 after selection, but coxph does not contain info of cohort size=4028. The rough estimate between coxph() and cch() is same, but the lower and upper CI and P-value are a little different. Can we

On Fri, 2017-12-15 at 10:16 +1300, Garming Sam via samba wrote: > Hi, > > We've just had another person notice this issue and we're currently > looking into it. If we need someone to test the patch fix, or if we need > more information, we can let you know. > > I've filed a bug as well: > > https://bugzilla.samba.org/show_bug.cgi?id=13191 Attached is a

----- begin included message In case cohort study, we can fit proportional hazard regression model to case-cohort data. In R, the function is cch() in Survival package Now I am working on case cohort analysis with time dependent covariates using cch() of "Survival" R package. I wonder if cch() provide this utility or not? The cch() manual does not say if time dependent covariate is

Hi, the following strange error appears when I use qda: > qda1 <- qda(as.data.frame(mfilters[cvtrain,]),as.factor(traingroups)) Error: function is not a closure That's also strange: > qda1 <- qda(mfilters[cvtrain,],as.factor(traingroups)) Error in qda.default(mfilters[cvtrain, ], as.factor(traingroups)) : length of dimnames must match that of dims Some backgroud: >

Full_Name: Kevin Wright Version: 1.4 OS: Windows 95 Submission from: (NULL) (170.54.59.160) Note: This was the 1.4 build for Windows that Brian Ripley made available. The first example in the help for reshape doesn't work for me. When I cut and paste, this is what happens: > data(Indometh,package="nls") > summary(Indometh) Subject time conc

Hola a todos, Me gustaria saber como puedo evaluar / obtener el resultado guardado en un objeto cuando este se ejecuta. Quiero decir, el resultado de la siguiente expresion: > paste("res_UP_sep_dec_","16_17_18",sep="") [1] "res_UP_sep_dec_16_17_18" es un objeto (una lista) cuyo nombre es precisamente res_UP_sep_dec_16_17_18 >

>> And I have a quite fast algo that I believe is simpler than [Budimlic02] >> and I can share it with you :) > > Do you have a paper on this? I'd be interested in seeing it. > Yes, I have a tech report on this page: http://compilers/fernando/projects/soc/ and I have submitted a paper to SAS, and now I am waiting for the review. The coalescing algorithm is described in

Who's your advisor? -scooter (aka "Dr. B. Scott Michel, UCLA CS 2004" :-) On Apr 12, 2007, at 4:39 PM, Fernando Magno Quintao Pereira wrote: > >>> And I have a quite fast algo that I believe is simpler than >>> [Budimlic02] >>> and I can share it with you :) >> >> Do you have a paper on this? I'd be interested in seeing it.

> And I have a quite fast algo that I believe is simpler than [Budimlic02] > and I can share it with you :) Do you have a paper on this? I'd be interested in seeing it. -Tanya > > Fernando > _______________________________________________ > LLVM Developers mailing list > LLVMdev at cs.uiuc.edu http://llvm.cs.uiuc.edu >

Dear all, I am using poly() in lm() in the following form. 1> DelsDPWOS.lm3 <- lm(DelsPDWOS[,1] ~ poly(DelsPDWOS[,4],3)) 2> DelsDPWOS.I.lm3 <- lm(DelsPDWOS[,1] ~ poly(I(DelsPDWOS[,4]),3)) 3> DelsDPWOS.2.lm3 <- lm(DelsPDWOS[,1]~DelsPDWOS[,4]+I(DelsPDWOS[,4]^2)+I(DelsPDWOS[,4]^3)) 1 and 2 lead to identical but wrong results. 3 is correct. Surprisingly (to me) the residuals

I am analyzing data from an experiment with two factors: Carbon (+/-) and O3 (+/-), with 4 replicates of each treatment, and 4 harvests over a year. The treatments are assigned in a block design to individual Rings. I have approaches this as a repeated measures design. Fixed factors are Carbon, O3 and Harvest, with Ring assigned as a random variable. I have performed repeated measures analysis

Hello, I?m using zero-inflated Poisson regression via the zeroinfl() function to analyze data on seed-set of plants, but for some reason, I don?t seem to be getting the output for all three levels of my two categorical predictors. More about my data and model: My response variable is the number of viable seeds (AVInt), and my two categorical predictors are elevation (Elev) and Treatment

Hello Guys, I'm using the package RankProduct to determine the differential expression of genes in 2 classes. The importance of the thing is to know: - How can I determine which class is the class of up-regulated genes and which to the down regulated gene - How can i interpret the results : $pfp, $RPs, $RPrank, Orirank,

All, Given the following commands: > ag.m35<-read.table("m35.txt",header=TRUE,sep=",") > ag.m35.lp<-subset(ag.m35, race=="lp") > aov.m35.lp=aov(year~time,data=ag.m35.lp) > anova.m35.lp=anova(aov.m35.lp) > anova.m35.lp Analysis of Variance Table Response: year Df Sum Sq Mean Sq F value Pr(>F) time 1 11.56 11.5649 1.7282

On Thu, 30 Jun 2005, Tanu Sharma wrote: > I am compiling SPEC 2000 benchmarks with llvm .Got stuck with > calculating "execution time" of all the .bc and native files. > > The log for nightly test itself gives execution times but I am passing > the bytecode files to my pass which gives another bytecode file.I have > to calculate execution time of such bytecode and

Hello , I am compiling SPEC 2000 benchmarks with llvm .Got stuck with calculating "execution time" of all the .bc and native files. The log for nightly test itself gives execution times but I am passing the bytecode files to my pass which gives another bytecode file.I have to calculate execution time of such bytecode and native files as well.If i simply do this: time lli

Thanks Andrew, much better than vanilla 4.7.3 but it is still performs much worse than vanilla 4.7.0 especially for DLZ zones. $ dig @127.0.0.1 google.com $ perl -e "foreach(1..500) { system('dig', '@127.0.0.1', 'google.com'); }" | grep time > 4.7.0_external.log $ perl -e "foreach(1..500) { system('dig', '@127.0.0.1',

Hello All, Thanks for the reply.I can generate the reports by compiling Spec through llvm, but that couldn't resolve my problem. I m trying to determine execution time for the bytecode and native files , which are obtained as a result of running my pass over the original bytecode .I am running these experiments on spec benchmark. In SPEC we have command line tools such as runspec where