similar to: Expression for pedices

I need to extract the p value from a ANOVA done with lm model fitting <- lm(var ~ group) Sfitting <- summary(fitting) Sfitting[10][1] gives the F value and the degrees of freedom but I am not able to get the p value. The function df should give a p value given a F but I am not able to make it work. I found only something about aov in the R help and I am not able to make it work Massimo

I am writing some code to obtain publication-like plots (like the ones can be obtained with SigmaPlot). I am not able to find the dimensions of a point. Parameters like cex make you able to make it bigger or smaller but I need to know the exact dimension (in relative coordinates : if my plot is 640x480 a dot plotted with cex=1.5 will be...) Massimo Cressoni

I am not able to compile R under ubuntu Edgy. ./configure error message is configure: error: --with-x=yes (default) and X11 headers/libs are not available I do not know which library is missing (I tried to install some X11 libs but were already on the system). Massimo Cressoni

I am trying to write a function that fits a sigmoid given a X and Y vector guessing the start parameters. I use nls. What I did (enclosed) seems to work well with many data points but if I want to fit small vectors like : pressure <- c(5,15,9,35,45) gas <- c(1000,2000,3000,4000,5000) it do not work. The help page says that it do no not work on zero residual data. Massimo Cressoni

I need to perform the Exact Wilcoxon Mann-Whitney on a subset of my database. Assuming that IPPO is my data frame and IPPOBIS is the subset my variable still have 3 different levels and the function wilcox_test (package "coin") does not accept it. I do not know how to overcome this problem. ippo <- c(rep("A",10),rep("B",10),rep("C",10)) ippo2 <-

I have 2 DCs running 4.1.12 that have stopped replicating and are producing these error messages: from dc1: Failed to commit objects: WERR_GENERAL_FAILURE/NT_STATUS_INVALID_NETWORK_RESPONSE [2014/09/23 10:43:35.530000, 0] ../source4/rpc_server/drsuapi/getncchanges.c:1646(dcesrv_drsuapi_DsGetNCChanges) ../source4/rpc_server/drsuapi/getncchanges.c:1646: DsGetNCChanges 2nd replication on

Dear R users, When I run the code below, I get the error "Error in integrate(integrand, 0, Inf) : non-finite function value". The code works if the function returns only "sum(integ)". However, I want to add "cmh" to it. When I add "cmh" I get that error. I can't figure out why this is happening because my integrate function has nothing to do with

In SAS, for a two-way (or 3-way, stratified) table, the CMH option in SAS PROC FREQ gives 3 tests that take ordinality of the factors into account, for both variables, just the column variable or neither. Is there an equivalent in R? The mantelhaen.test in stats gives something quite different (a test of conditional independence for *nominal* factors in a 3-way table). e.g. I'd like to

Dear R users, I have a microRNA (similar to microarray) data set with gene expressions and group variables. I'd like to perform (class) prediction analysis to see which genes or gene sets predict the group variable well. I'm at the beginner level. Is there any free software in R or in point/click format that is commonly used and highly recommended? I googled this topic and found

Dear R users, I'm wondering how I can generate an arbitrary number of loops in R. For instance, I can generate two "for" loops to get ICC among any two-way combination among 10 variables. Here is the code n<-10 for (i in 1:(n-1)) { for (j in (i+1):n) { icc(cbind(DATA[,i],DATA[,j])) } } If I need three-way combination, then a code with three "for" loops will be:

Hello, I've tried to analyze some data with a CMH test. My 3 dimensional contingency tables are 2x2xN where N is usually between 10 and 100. The problem is that there may be 2 strata with opposite counts (the 2x2 contigency table for these are reversed), producing opposite odds ratios that cancle out in the overall statistics. These opposite counts are very important for my analysis, since

Dear R users, I'm wondering how I can generate an arbitrary number of loops in R. For instance, I can generate two "for" loops to get ICC among any two-way combination among 10 variables. Here is the code n<-10 for (i in 1:(n-1)) { for (j in (i+1):n) { icc(cbind(DATA[,i],DATA[,j])) } } If I need three-way combination, then a code with three "for" loops will be:

I would really appreciate it if someone can give suggestions on how to do spectra fitting in R using ordinary least square fitting and non-negativity constraints. The lm() function works well for ordinary least square fitting, but how to specify non-negativity constraints? It wouldn't make sense if the fitting coefficients coming out as negative in absorption spectra deconvolution. Thanks.

Hi all, Does anyone know if an R function for the Mantel-Haenszel chi-square for a 2x2 table exists? I’ve also seen it called the randomization Q statistic. Note that I’m not looking for the Cochran-Mantel-Haenszel…I did see that out there as cmh.test. Thanks in advance, Jeff Internal Virus Database is out-of-date. Checked by AVG Free Edition. 11:27 AM [[alternative HTML version

Thanks. Eric Leifer, Ph.D. Office of Biostatistics Research National Heart, Lung, and Blood Institute 6701 Rockledge Drive, MSC 8217 Bethesda, MD 20892-8217 Email: LeiferE@NHLBI.NIH.GOV Phone: 301-435-0436 FAX: 301-480-1862

Well, I've "fixed" two problems so far with getting things built on OS X. If whoever did the OS X work is reading this, feel free to tell me where/how I've screwed up =) The first problem is that it couldn't find MacTypes.h in the path. MacTypes.h is simply a one line file that includes <CoreServices/CoreServices.h>, so I just replaced the MacTypes.h with that. In

Just a general question concerning the woolf test (package vcd), when we have stratified data (2x2 tables) and when the p.value of the woolf-test is below 0.05 then we assume that there is a heterogeneity and a common odds ratio cannot be computed? Does this mean that we have to try to add more stratification variables (stratify more) to make the woolf-test p.value insignificant? Also in the

Hi there, is there any way to compare 2 odds ratios? I have two tests that are supposed to detect a disease presence. So for each test, I can compute an odds ratio. My problem is how can I compare the 2 tests by testing whether the 2 odds ratios are the same? Appreciate

Dear all, I am trying to look for a built in function that performs the cochran Q test. that is, cochranq.test(X) where X is a contingency table (maybe a matrix or data.frame). The output will naturally be the test statisitcs, p-value, etc. A quick search on Google gives me the cochran.test in the 'outlier' package, but I had a look at the description of the test and it doesn't look

Hello group - Never heard a peep about this....wasn't sure it was actually posted. Anyone seen this before? Thx RSV869 Begin forwarded message: Date: Sat, 26 Oct 2013 11:45:35 -0400 From: Gmail <rsv869 at gmail.com> To: nouveau at lists.freedesktop.org Subject: Dual monitor question Hello Nouveau group - I am an OpenSuSE 12.3 user that's been struggling with a specific