similar to: Partial R-square in multiple linear regression

I am trying to figure out how to calculate p-values for the difference in prevalence of a risk factor between men and women. For example, I find that 277 out of 710 male patients and 125 out of 305 female patients have obesity, what is the p-value for their difference? If there is a package that can calculate this in bulk, I would appreciate to learn about it! Thank you [[alternative HTML

Dear users of metafor, I am working on a meta-analysis using the metafor package. I have a excel csv database that I am working with. I am interested in pooling the effect measures for a particular subgroup (European women) in this csv database. I am conducting both sub-group and meta-regression. In subgroup-analyses, I have stratified the database to create a separate csv file just for European

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

Hello: Forgive me, this is surely a simple question but I can't figure it out, having consulted the help archives and "Data Manipulation With R" (Spector). I have a list of 11 data frames with one common variable in each (prov). I'd like to use lapply to go through and recode one particular level of that common variable. I can get the recode to work, but it only returns the

Hey I'm not really sure what I should put on here, but I am having trouble with my R code. I am trying to get the p-values, R^2s etc for a number of different groups of variables that are all in one dataset. This is the code: #Stand counter st<-1 #Collections stands<-numeric(67) slopes<-numeric(67) intercepts<-numeric(67) mses<-numeric(67) rsquares<-numeric(67)

Dear all, I have a question regarding the output of the coxph function. What does the 'max possible' exactly mean in the output below? Many thanks. coef exp(coef) se(coef) robust se z Pr(>|z|) smocc_zyban -0.4384 0.6451 0.8667 0.9473 -0.463 0.644 self 1.1857 3.2728 0.1405 0.1443 8.216 2.22e-16 *** smocc_zyban:self

Hi, I have something like this: > x <- 1:10 > y2 <- 30+5*x+rnorm(x,sd=3) > y <- c(y1,y2) > x <- c(x,x) > plot(x,y) > x <- 1:10 > y1 <- 1+5*x+rnorm(x,sd=2) > y2 <- 30+5*x+rnorm(x,sd=5) > y <- c(y1,y2) > x <- c(x,x) > f <- factor(rep(c("a","b"),c(10,10))) > m <- lm(y~x+f) > anova(m) Analysis of Variance

Hello there, I?m not sure if this is the right mailing list, but I?d like to see if we can add Eucalyptus to the IaaS page here: http://libvirt.org/apps.html#iaas Please let me know what is required in order to add us? Thanks very much, *_________________________* *Eric Choi - *Product Marketing *Eucalyptus Systems* www.eucalyptus.com +1 508 348-9778 Follow us on Twitter

Dear R-users, This should be simple but still eludes me: Given the following tmp<-as.data.frame(matrix(c(44, 10, "abc", 1, 44, 10, "def", 1, 44, 12, "abc", 2), 3, 4, byrow=T)) I want to expand the data to the following form: V1 V2 V3 V4 1 44 10 abc 1 2 44 10 def 1 3 44 12 abc 1 4 44 12 abc 1 The last row of the original df was duplicated the row by the

Hello, I found that round() does not behave as I expected. Have you had similar experience as following? > x<-seq(0.5,10.5,by=1) > x [1] 0.5 1.5 2.5 3.5 4.5 5.5 6.5 7.5 8.5 9.5 10.5 > round(x) [1] 0 2 2 4 4 6 6 8 8 10 10 > cbind(x,round(x)) x [1,] 0.5 0 [2,] 1.5 2 [3,] 2.5 2 [4,] 3.5 4 [5,] 4.5 4 [6,] 5.5 6 [7,] 6.5 6

Amigos de la erre. He creado mi primer bucle con for para entrenar unos modelos con GAM. La respuesta es quasipoisson porque estoy trabajando con densidades de peces. Sin embargo, tengo un problema, no se muy bien como añadir el offset a la formula siguiente cuando creo el bucle. GAM.A1 <-gam ((DYO)~s(DMA,k=4)+ s(WOD,k=4)+s(CIN,k=4)+s(DRA,k=4)+s(DBR,k=4)

Dear all, My machine is SUN Java Workstation 2100 with 2 AMD Opteron CPUs and 16GB RAM. R is compiled as 64bit by using SUN compilers. I trying to fit quantile smoothing on my data and I got an message as below. > fit1<-rqss(z1~qss(cbind(x,y),lambda=la1),tau=t1) Error in as.matrix.csr(diag(n)) : cannot allocate memory block of size 2496135168 The lengths of vector x and y are

Hello all, I have a data frame with year, month, species, fishing gear and catch (Y, M, S, F, C) and I want the sum of C by Y for species "A" and fishing gear "trawl". I tried things like aggregate(C[S=="A" & F=="trawl"], list (Year = Y[S=="A" & F=="trawl"]), fun=sum), but it didn't worked. To overcome this problem I did

Dear list: I'm using bwplot to compare concentrations by location and treatment as in: # using built in data bwplot( conc ~ Type : Treatment, data = CO2 ) I would like the order of the plots to be: 3,4,1,2. I can't seem to figure this out with index.cond or permc.cond. Any help is appreciated! Brad Eck [[alternative HTML version deleted]]

Hi, I want to have a horizontal line on a mosaic plot with "vcd" package. This would give me an idea where is 0.5 proportion in a cell. Using "mosaicplot" function of "graphics" package, I can draw a line using "abline." But, with "mosaic" function of "vcd" package, I have tried to use "abline" function, which complains

Hi, I am trying to run the command "forward.sel.par," however I receive the error message: "Error: could not find function 'simpleRDA2'." I have the vegan library loaded. The documentation on "varpart" has not helped me to understand why I cannot call this function. Maybe I am missing something obvious because I am still an 'R' novice. Below is a

Hi there, As you know, $HOME is generally located at "/home/$username" by default. I would like to re-locate all users' $HOME directories to something like "/export/home/$username" without having a hassle/trouble. Initially, I've thought of just copying them to the new directory (under /export/home/xxx), but guessed it might trouble for the normal use (I'm pretty

> On Sep 16, 2016, at 2:23 PM, Sanjoy Das <sanjoy at playingwithpointers.com> wrote: > > Hi Mehdi, > > Mehdi Amini wrote: > > > >> On Sep 16, 2016, at 1:56 PM, Kevin Choi via llvm-dev > >> <llvm-dev at lists.llvm.org <mailto:llvm-dev at lists.llvm.org>> wrote: > >> > >> > int mat[9][9][9]; > >> > for (p =

Hi Quentin, MyInst is a custom instruction that has implicit-defs of fixed registers. The implicit-defs are seen at the end of Instruction Selection. I'd like to add a report, but I am working on an out-of-tree backend based on 7.0. I can try to help reduce the testcase down. Filed https://bugs.llvm.org/show_bug.cgi?id=41790 Regards, Kevin On 2019-05-07 3:45 p.m., Quentin Colombet wrote:

Dear R users, I tried to find a solution in the search list, but I cannot find it. I would like to read a .txt file with, let say, three variables, with two of which have repeated values in a number a columns. An example: The variables: Treat, x1, x2. The values: A 2.5 3.4 2.7 5.6 5.7 5.4 10.1 9.4 ... B 5.3 5.4 6.5 7.5 1.3 4.5 10.5 4.1 ... ... In the first column, the letters represent the