similar to: Regression versus functional/structural relationship?

Hello R-devel, I am currently attempting to implement an API similar to data.table wherein single bracket subsetting can accept an unquoted expression to be evaluated in the context of my object. A simple example from the data.table package looks like this: DT <- data.table(col1 = c('a', 'b', 'c'), col2 = c('x', 'y', 'z')) DT[col1 ==

I know that there are quite a few packages out that there for cluster analysis. The problem that I am facing is finding a package that will not incorporate all my samples into clusters but just the samples that fit a threshold (that I have not set yet and may need help finding the right level) for genotyping. It should be able to "no call" samples outside the clusters. It also needs to

Hi. I have data which i would want to assess the degree of agreement between two assays, e.g., to evaluate reproducibility or for inter-rater reliability. I have used the Pearson product-moment correlation coefficient. It looks good ranginging between 0.90 to 0.998. Though this looks good. I am told the Concordance correlation coefficient will give a better picture of how reproducible the assay

Dear all, Is there anyway too generate MA plot for 2 qPCR assays (an array of 2x 400). -- View this message in context: http://r.789695.n4.nabble.com/Bioconductor-MA-plot-for-qPCR-array-tp4182805p4182805.html Sent from the R help mailing list archive at Nabble.com.

Greetings ~ I need some assistance determining an appropriate approach to analyzing multivariate binary response data, and how to do it in R. The setting: Data from an assay, wherein 6-hours-post-fertilization zebrafish embryos (n=20, say) are exposed in a vial to a chemical (replicated 3 times, say), and 5 days later observed for the presence/absence (1/0) of defects in several organ systems

Dear R-helpers, I am developing a Mixed-Effects model for a study of immunoassays using 'lme4' library. The study design is as follows: 10 samples were run using 7 different immunoassays, 3 times each, in duplicates. Data attached. I have developed the following model: c.lme <- lmer(Result~SPL + (SPL|Assay/Run) -1, data=data) This model has excellent predictions - the Rsquared of

Hi dear R-listers, I'm trying to fit a 3-level model using lme in R. My sample size is about 2965 and 3 factors: year (5 levels), ssize (4 levels), condition (2 levels). When I issue the following command: > lme(var~year*ssize*condition,random=~ssize+condition|subject,data=smp,method ="ML") I got the following error: Error in logLik.lmeStructInt(lmeSt, lmePars) :

Hi list: I work with a lot of laboratory analytical data and I often have inconsistent names of files and variables within those files so I wrote this sourcearg function to facilitate handling file and variable names as both character and R names. The source of the function is given below sourcearg=function(arg){ su=sprintf("%s\n",arg) sink("TMPTMPRarg")

The Gottardo lab at the Fred Hutchinson Cancer Research Center in Seattle has an opening for a programmer. To apply: https://erecruit.fhcrc.org/psp/RECRUIT/EMPLOYEE/HRMS/c/HRS_HRAM.HRS_CE.GBL?Page=HRS_CE_JOB_DTL&Action=A&JobOpeningId=24122&SiteId=2&PostingSeq=1 The position description is below: About Us Fred Hutchinson Cancer Research Center, home of three Nobel laureates,

Hi: I have data from an assay in the form of two vectors, one is response and the other is a predictor. When I attempt to fit a 5 parameter logistic model with "nls", I get converged parameter estimates. I also get the same answers with "gnls" without specifying the "weights" argument. However, when I attempt to use the "gnls" function and try to

There may be benefits to having a machine learning method that explicitly targets covariate balance. We have experimented with optimizing the weights directly to obtain the best covariate balance, but got some strange solutions for simple cases that made us wary of such methods. Machine learning methods that yield calibrated probability estimates should do well (e.g. those that optimize the

Hi, Given a positive integer N, and a real number \lambda such that 0 < \lambda < 1, I would like to generate an N by N stochastic matrix (a matrix with all the rows summing to 1), such that it has the second largest eigenvalue equal to \lambda (Note: the dominant eigenvalue of a stochastic matrix is 1). I don't care what the other eigenvalues are. The second eigenvalue is

R Package Developer / Consulting Statistician Burlington, Vermont Precision Bioassay is a small and growing statistical consulting firm specializing in the design and analysis of biological assays of pharmaceutical proteins and vaccines. We are building a turnkey software solution for cell culture bioassay that integrates optimal design of bioassay experiments, implementation of designs using

Hi, Where can I find the source code for nlm()? I dowloaded the R2.8.1.tar.gz file and looked at all the .c and .f files, but couldn't find either nlm.c or nlm.f There is an nlm.r file, but that is not useful. Thanks for any help, Ravi. ---------------------------------------------------------------------------- ------- Ravi Varadhan, Ph.D. Assistant Professor, The Center on Aging

Hi, I am trying to do the following computation: p <- rep(0, n) coef <- runif(K+1) U <- matrix(runif(n*(2*K+1)), n, 2*K+1) for (i in 0:K){ for (j in 0:K){ p <- p + coef[i+1]* coef[j+1] * U[,i+j+1] } } I would appreciate any suggestions on how to perform this computation efficiently without the "for" loops? Thank

Hi, I have a simple problem that I would appreciate getting some tips. I am using the "truehist" function within an "apply" call to plot multiple histograms. I can't figure out how to get truehist to use the column names of the matrix as the labels for the x-axis of the histograms. Here is a simple example: X <- matrix(runif(4000),ncol=4) colnames(X)

Hi, We (Paul Gilbert and I) have just released a new R package on CRAN called "BB" (stands for Barzilai-Borwein) that provides functionality for solving large-scale (and small-scale) nonlinear system of equations. Until now, R didn't have any functionality for solving nonlinear systems. We hope that this package fills that need. We also have an implementation of the

Hello everyone: I'm a new member of this group. I have a question about "oneway.test". When I use "anova(lm(....))" to analysis the ANOVA, I can get the information about Sum Sq and Mean Sq. (The R code and the results are as follows.)

Hi, We have some sample data from the US census, and we know the marginal totals for the population. We need to make the population estimates add up to the correct sums. I have two questions: Is there some package in R which does this adjustment, by any means? Is there some more modern reference for this problem than Deming's 1943 monograph, ``Statistical Adjustment of Data''?

Hi, Is there an R implementation of least angle regression for binary response modeling? I know that this question has been asked before, and I am also aware of the "lasso2" package, but that only implements an L1 penalty, i.e. the Lasso approach. Madigan and Ridgeway in their discussion of Efron et al (2004) describe a LARS-type algorithm for generalized linear models. Has