similar to: How to get the tuning parameter lamda in storey's qvalue package

I am using storey's qvalue package but I keep on getting errors. Why is this? > qvalue(p, lambda=0.5)$pi0 [1] "ERROR: p-values not in valid range." Error in qvalue(p, lambda = 0.5)$pi0 : $ operator is invalid for atomic vectors -- Thanks, Jim. [[alternative HTML version deleted]]

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

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

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,

Dear all I am so glad the R can provide the efficient calculate about eigenvector and eigenvalue. However, i have some puzzle about the procedure of eigen. Fristly, what kind of procedue does the R utilize such that the eigen are obtained? For example, A=matrix(c(1,2,4,3),2,2) we can define the eigenvalue lamda, such as det | 1-lamda 4 | =0 | 2 3-lamda | then

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

I am trying to it a particular nonlinear model common in Soil Science to moisture release data from soil. I have written the function as shown below according to the logist example in Ch8 of Pinheiro & Bates. I am getting the following error (R version 2.1.1) *Error in qr(attr(rhs, "gradient")) : NA/NaN/Inf in foreign function call (arg 1)* Below is the function and data. /#

Hi I have been having issue with a ML estimator for Jump diffusion process but know I am get little error I didn't notice before like I am try to create a vector > #GBMPJ MLE Combined Ph 1 LR > # > n<-length(combinedlrph1) > j<-c(1,2,3,4,5,6,7,8,9,10) Error in c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10) : unused argument(s) (3, 4, 5, 6, 7, 8, 9, 10) >

Hi all, I am a bit mystified by this error message that I get when I try to apply lmer() to a simple dataset with one between factor (age) and one within factor (item): "$ operator is invalid for atomic vectors" I'll just provide the code, because I don't see where the problem is: library(lme4) options(contrasts=c("contr.helmert","contr.poly")) data =

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, ...)

I have a set of data for US t-bill returns and US stock returns frm 1980-2012. I am trying to bootstrap the data and obtain the minimum variance portfolio and repeat this portfolio 1000 times. However I am unable to get the correct code function for the minimum variance portfolio. When I tried to enter Opt(OriData+1, 1, 5, 0), I get "error:subscript out of bounds" Please help!

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

Hello guys.I am working on implementing the PL2 weighting scheme of the DFR framework by Gianni Amati. It uses the Poisson approximation of the Binomial as the probabilistic model (P), the Laplace law of succession to calculate the after effect of sampling or the risk gain (L) and within document frequency normalization H2(2) (as proposed by Amati in his PHD thesis). The formula for w(t,d) in

Hi, I have been using maxLik to do some MLE of Geometric Brownian Motion Process and everything has been going fine, but know I have tried to do it with jumps. I have create a vector of jumps and then added this into my log-likelihood equation, know I am getting a message: NA in the initial gradient My codes is hear # n<-length(combinedlr) j<-c(1,2,3,4,5,6,7,8,9,10)

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

Dear Sir/madam, I'm getting a problem with a R-code which calculate Fisher Information Matrix for Hybrid Censored Weibull Distribution. My problem is that: when I take weibull(scale=1,shape=2) { i.e shape>1} I got my desired result but when I take weibull(scale=1,shape=0.5) { i.e shape<1} it gives error : Error in integrate(int2, lower = 0, upper = t) : the integral is probably

Hello. Here is my R code. I used the functional data . Now I need to use the functional data by applying the kernels instead of the xi, yi functions. Bonjour. Voici mon code en R . J'ai utiliser les donn?es fonctionnelles . Maintenant j'ai besoin d'utiliser les donn?es fonctionnelles en appliquant les noyaux ? la place des fontions xi, yi library(MASS)

Hi list, I have a group of data. It looks like they follow a exponential distribution. In R, how can I esimate lamda, that is the rate in pexp, of the distribution and can I use Kolmogorov-Smirnov for hypothesis testing in such a situation? I have read the "8.2 Examing the distribution of a set of data" of "An Introduction to R" but I did not find any clues on this issue.

Good evening dear administrators, It is with pleasure that I am writing to you to ask for help to finalize my R programming algorithm. Indeed, I attach this note to my code which deals with a case of independence test statistic . My request is to introduce the kernels using the functional data for this same code that I am sending you. So I list the lines for which we need to introduce the