similar to: qvalue

Dear colleagues, This is not strictly a R question, but I hope it is ok to ask on the list. I fed a vector of p-values from about 20 million anova tests to the package q-value and obtained this output: > qsummary(asso_p.qvalue) Call: qvalue(p = asso_p.vec) pi0: 1 Cumulative number of significant calls: <1e-04 <0.001 <0.01 <0.025 <0.05 <0.1 <1

Dear All, In Storey's estimator of the proportion of true nulls, the estimator depends on the tuning parameter lamda. Suppose now that an estimator of this proportion has been obtained by the qvalue package, what is the lamda that corresponds to the estimate? How to get this lamda? Thanks, -Chee [[alternative HTML version deleted]]

I couldnt find a previous posting on this in the archives, maybe it has already been mentioned. If you use a calculation to generate n observations in random number generators and you don't round to the nearest integer you may be generating n-1 numbers not n numbers as you thought depending on the storage precision of the calculation. e.g. > m <- 1000 > pi0 <- 0.9 >

Good Morning Everyone, I was going through my list of FTBFS packages today and I fixed all my R packages but 2: qvalue and xtable. qvalue requires ggplot2 xtable requires: lsmeans, spdep, splm, sphet, plm I am not doing any R anymore these days and in fact spot has been the one maintaining most of my R packages these days (thanks spot!!), so I am not really interested in maintaining more R

Mark, there is a fdr website link via Yoav Benjamini's homepage which is: http://www.math.tau.ac.il/%7Eroee/index.htm On it you can download an S-Plus function (under the downloads link) which calculates the false discovery rate threshold alpha level using stepup, stepdown, dependence methods etc. Some changes are required to the plotting code when porting it to R. I removed the

Had some packages fail install so I updated to today's R-beta release. On updating packages the following packages still fail. sessionInfo follows. Mark Warning messages: 1: In install.packages(update[instlib == l, "Package"], l, contriburl = contriburl, : installation of package 'MANOR' had non-zero exit status 2: In install.packages(update[instlib == l,

Hello all, I'm trying to use the package qvalue. Upon trying to install it I'm running into trouble. I re-downloaded and installed the most recent CRAN distribution of R which announced itself to have tcltk included. Now when I try to load qvalue I'm getting errors that seem to indicate that I need tcltk2. I tried going to get that package and loading it but both qvalue and tcltk seem

You asked the same question on the Bioconductor mailing list back in August. At that time, you suggested yourself a solution for how the adjusted p-values should be interpreted. I answered your query and told you that your interpretation was correct. So I'm not sure what more can be said, except that you should read the article Wright (1992), which is cited in the help entry for p.adjust(),

Dear all, I need to do some calculation where the code used are below. I get error message when I choose k to be large, say greater than 25. The error message is "Error in integrate(temp, lower = 0, upper = 1, k, x, rho, m) : the integral is probably divergent". Can anyone give some help on resolving this. Thanks. Hannah m <- 100 alpha <- 0.05 rho <- 0.1 F0

Hello, I would like to use the parfm package: https://cran.r-project.org/web/packages/parfm/parfm.pdfhttps://cran.r-project.org/web/packages/parfm/parfm.pdf in my work. This package fits parametric frailty models to survival data. To ensure I was using it properly, I started by running some small simulations to generate some survival data (without any random effects), and analyse the data using

Dear list, I have been trying to run the function "qvalue" under the package qvalue on a vector with about 20 million values. > asso_p.qvalue<-qvalue(asso_p.vector) Error: cannot allocate vector of size 156513 Kb > sessionInfo() Version 2.3.1 (2006-06-01) i686-pc-linux-gnu attached base packages: [1] "methods" "stats" "graphics"

Hi, I am trying to install qvalue, however its giving installation error -> Error : package 'tcltk' does not have a name space ERROR: lazy loading failed for package ?qvalue? * removing ?/home/sbw/R/x86_64-unknown-linux-gnu-library/2.12/qvalue? The downloaded packages are in ?/tmp/RtmpKnS1X4/downloaded_packages? Warning message: In install.packages(pkgs = pkgs, repos = repos, ...)

Hi, I have a huge dataset (53 million records). I have to calculate the p and q values of my data. How can I do it in R or perl? I have downloaded R (I'm completely new to R). and the package qvalue but I don't understand how can I call/use qvalue package with R. When I type library(qvalue), it gives me an error that this package doesn't exist. What should I do? Thanks! -- View this

Hi, All I got trouble on installing the qvalue package. Error message: package 'tcltk' does not have a name space [cchen1 at ibibmem Yale_CB]$ R CMD INSTALL qvalue.tar.gz * installing to library '/cchome/cchen1/R/x86_64-unknown-linux-gnu-library/2.10' * installing *source* package 'qvalue' ... ** R ** data ** inst ** preparing package for lazy loading Error : package

Hi all, I want to draw a histgram for each row of a matrix and compare them. However the plot I got does not have the same y range and x range, which makes it difficult to make the comparison. Is there a easy way to fix the x range and y range in a xy plot for several plots, instead of specifying them for each plot. The following is my code for generalizing the matrix and draw the histogram.

Hi, I am writing some basic smoothers in R for cleaning some spectral data. I wanted to see if I could get close to matlab for speed, so I was trying to compare SparseM with Matrix to see which could do the choleski decomposition the fastest. Here is the function using SparseM difsm <- function(y, lambda, d){ # Smoothing with a finite difference penalty # y: signal to be smoothed #

Hi all, Sorry for the re-post, I sent my previous e-mail before it was complete. I am trying to model seroprevalence using the differential equation: dP/dt = beta*seronegative*.001*(seropositive)-0.35*(0.999)*(seropositive)-r*seropositive. I would like to estimate my two parameters, beta and r, using maximum likelihood methods. I have included my code below:

Hi R users, I have a couple of questions about some problems that I am facing with regard to numerical integration and optimization of likelihood functions. Let me provide a little background information: I am trying to do maximum likelihood estimation of an econometric model that I have developed recently. I estimate the parameters of the model using the monthly US unemployment rate series

I?m working on an R-implementation of the simulation-based finite-sample null-distribution of (R)LR-Test in Mixed Models (i.e. testing for Var(RandomEffect)=0) derived by C. M. Crainiceanu and D. Ruppert. I'm in the beginning stages of this project and while comparing quick and dirty grid-search-methods and more exact optim()/optimize()-based methods to find the maximum of a part of the

I want to fit some p-values to a beta distribution. But the problem is some of the values have 0s and 1's. I am getting an error if I use the MASS function to do this. Is there anyway to get around this? -- Thanks, Jim. [[alternative HTML version deleted]]