similar to: isotone

Hi All, How can I get the expected value from a discrete cdf? Is there any R function that can do this? Thanks, Annie [[alternative HTML version deleted]]

Hi, All, How can I get the indices of the minimum elements in a matrix without using a loop? For example, if the matrix is 4 5 2 2 8 9 5 2 3 Then I want to output (1,3), (2,1), (3,2). Thanks, Annie [[alternative HTML version deleted]]

Hi All, I have a 2-dim density defined on 0<x<1, 0<y<1, x<y. I know the exact formula of the density. How can I visualize it? What plot functions can I use? Thanks, Annie [[alternative HTML version deleted]]

Hi, R users, I am trying to use glm to do logistic regression. I know generally when I have two covariates, say x1 and x2, then I do fit <- glm(y~x1+x2,famliy='binomial') But now my covariates form a n*p matrix, say x, so actually each column is a covariate. So I think I should do fit <- glm(y~x,family='binomial') Then I need to predict new data. How should I write the

To Whom It May Concern: I have an R package and want to put this package be part of R-project and available to anyone who is interested in. The R package is created for my paper, titled "Acceptance Sampling Plans from Truncated Life Tests Based on the Birnbaum-Saunders Distribution for Percentiles". The paper has been accepted by Communications in Statistics: Simulation and

The homals package http://www.jstatsoft.org/v31/i04 will get a major programming overhaul. This will take some time, but what's a few years on a 40-year project. Suggestions from the audience are welcome. homals() has a core loop over the m variables in which 1. tapply is used to compute category quantification (centroids) 2. category quantifications are then adjusted to satisfy the

The homals package http://www.jstatsoft.org/v31/i04 will get a major programming overhaul. This will take some time, but what's a few years on a 40-year project. Suggestions from the audience are welcome. homals() has a core loop over the m variables in which 1. tapply is used to compute category quantification (centroids) 2. category quantifications are then adjusted to satisfy the

Hi, R users, What may be the best function in R to do variable selection in logistic regression? I have the same number of variables as the number of samples, and I want to select the best variablesfor prediction. Is there any function doing forward selection followed by backward elimination in stepwise logistic regression? Thanks, Annie [[alternative HTML version deleted]]

After experiencing these symptoms through many of the release candidates, we've been at v1.0.0 now for a few weeks, and they're still there. Perhaps a config issue with my setup (mbox locks?). It's really too disrupting to continue using this way. I've included symptoms, client info, mta info, log entries and config file. What other information can I provide? Please

Hi In rtc_set_time, mktime is called to calculate seconds since 1970/01/01, input parameters of mktime are required to be in normal date format. Such as: year=1980, mon=12, day=31, hour=23, min=59, sec=59. However, the current input parameter of mktime is tm->tm_year, and it is the number of years since 1900. (For example, if current time is 2012/12/31, and tm->tm_year is 112). This is

Hi, I'm having a problem getting my Windows machines to access shares in Samba. When they browse to the Samba box it sometimes gives them an error saying that they don't have permission or that the server is unavailable. However this doesn't always happen and other times it lists the shares. When I try to access the shares it just prompts for the username/password over and over.

Hi, R users, Is there any package for penalized logistic regression with more than two response classes? I read the manual for stepPlr, but it seems it's only for binary case. Thank you, Annie [[alternative HTML version deleted]]

Hi, ALL, Is there any function in R that does the exact test for the matched pair proportions (one sided), which I assume is binomial(b+c, .5). Thanks, Annie [[alternative HTML version deleted]]

Hi R-help, ? I have fit a cox ph model to my data, but have beenreceiving an error when trying to fit a model to the nomogram. Here is the codeand corresponding error: ? ? >nomogramCF = nomogram(cph_age6_40avp4, +????????????????????lp.at= seq(-10,10,by =0.5),lp = TRUE, +?????????????????????? +??????????????????????funlabel="5year survival",

Dear Soren and R users: I am trying to use the summaryBy function with weights. Is this possible? An example that illustrates what I am trying to do follows: library(doBy) ## make up some data response = rnorm(100) group = c(rep(1,20), rep(2,20), rep(3,20), rep(4,20), rep(5,20)) weights = runif(100, 0, 1) mydata = data.frame(response,group,weights) ## run summaryBy without weights:

Well, I could argue that it's not *completely* OT since my question is motivated by an enquiry that I received in respect of a CRAN package "Iso" that I wrote and maintain. The question is this: Given observations y_1, ..., y_n, what is the solution to the problem: minimise \sum_{i=1}^n (y_i - y_i^*)^2 with respect to y_1^*, ..., y_n^* subject to the "isotonic"

Dear Sir, Very glad to send E-mali to you , Im Annie. I have encountered some problems on the way of learning opus, so I would like to consult you. I found "global_stack" was not released after malloc, and malloc memory was too large.So what I want to ask you is whether "global_stack" malloc can have less memory, since the value "GLOBAL_STACK_SIZE" can be

dear R-help: one of my students is struggling to test an ordered alternative hypothesis on a set of groups (e.g., mu_a <= mu_b <= mu_c). There has been some literature on this topic -- a lot of this goes back to Bartholomew (1961); Gaines and Rice (see refs below) are the ones who've popularized it in the ecology community. The topic is closely related to isotonic regression, but

Hi, I am a Ph.D. student from Québec, Canada. I’m a beginner with R and Bioconductor. Until now the only experience I have is in analyzing microarray data using affy and limma packages. Now I am trying to analyze Rat Gene 10 st arrays and I would like to run RMA analysis and Smyth moderated t test on those arrays. Since no cdf official package is available for those arrays, after reading many

Hello, I'm doing some graphics for a paper and a need customize such with filled region above the density curve. My attempts I get something very near what I need, but I don't solve the problem of use semitransparent filled. Below a minimal reproducible code. Someone has any idea? require(lattice) # toy data... dt <- expand.grid(A=1:2, B=1:3, y=1:50) dt$y <- rnorm(nrow(dt), dt$B,