similar to: A question about R image function

Dear friends, My dataset is like the following: x y mcpvalue 0.4603578 0.6247629 1.001 0.4603715 0.6247788 1.001 0.4603852 0.6247948 1.001 0.4110561 0.5664841 0.995 The x and y variables are unsorted. I use the function image(x,y,mcpvalue) to generate a plot, but the error is that "increasing 'x' and

I've data that look like: Outcome V1_AA V1_EU V1_NA V2_AA V2_EU V2_NA 0 0.046 1.001 0.954 0.045 1.001 0.954 0 0.007 1 0.993 0.007 1 0.993 1 1.774 0.217 0.009 1.774 0.217 0.009 1 0.004 1.996 0 0.004 1.996 0 1

Hello everyone, I use latex() (Hmisc) for report generation and thus have been affected by the problem with rounding decimals described, for example, in this post: http://thread.gmane.org/gmane.comp.lang.r.general/73287/focus=73287 In short, numbers often are printed with 15 or so decimals even if there far less significant digits. The problem has been confirmed by Frank Harrell and Thomas

Dear R-help readers, I am writing to you in order to ask you a few questions about distribution fitting in R. I am trying to find out whether the set of event interarrival times that I am currently analyzing is distributed with a Gamma or General Pareto distribution. The event arrival granularity is in minutes and interarrival times are in seconds, so the values I have are 0, 60, 120, 180, and

"make check" is failing on reg-test3.R with a message that survival was built with an older version of R. (On my Ubuntu 32 bit and Ubuntu 64 bit machines). Why would "make check" be looking anywhere that it would find something built with an older version of R? ~/RoboAdmin/R-3.0.0/tests$ tail reg-tests-3.Rout.fail > print(1.001, digits=16) [1] 1.001 > ## 2.4.1 gave

Hi, I am trying to recode the output from a matrix(here is a small snippet of it): HGlt10RawPerc2008[1:20] [1] -5 0 -1 -1 0 2 3 -5 -2 0 2 0 1 -2 3 0 4 1 4 2 Here is the code I am using to recode it: HGBlt10Points2008 = recode(HGlt10RawPerc

Dear Users, In case of ridge logistic regression, i want to calculate the optimum penalty using aic and bic criteria. Here is the sample code: fit <- lrm(RES ~CAT01+NUM01+NUM02+CAT02+CAT03+CAT04+NUM03+CAT05+CAT06+NUM04+ CAT07+CAT08+NUM05+NUM06, data = train.data, x = TRUE, y = TRUE) pentrace(fit, penalty = list(seq(.001, 5, by=.1))) output: Best penalty: penalty df 1.001

I am trying to write data to csv but I am having issues with the separations. Basically I have some results that I get with R that I copied and pasted into word and then saved as .txt I want to write the results to csv because it's easier to make tables in word (all I would have to do is copy and paste into a table, instead of typing everything out). I am able to write the data to the

Hi, dear all: I am a beginner. I appreciate any help or hint from you. I am trying to do calculation with matrices. I have 3 matrices. One is matrixA, 2nd is matrixB, and last is matrixC. Here is matrixA: 1.8511.40.0831.001 0.8771.30.1161.33 1.9021.21.1020.302 0.8640.1261.110.252 1.8230.2161.0020.307 Next is matrixB: 0.8761.770.1930.328 0.8911.0090.2381.004

Hi, Can anybody tell me what the message below means and how to overcome it. Thanks, Neil Warning message: X matrix deemed to be singular; variable 2 in: coxph(Surv(age_at_death, death) ~ project$pluralgp + project$yrborn + ......... >

Hi All, I'm trying to do the following constrained optimization example. Maximize x1*(1-x1) + x2*(1-x2) + x3*(1-x3) s.t. x1 + x2 + x3 = 1 x1 >= 0 and x1 <= 1 x2 >= 0 and x2 <= 1 x3 >= 0 and x3 <= 1 which are the constraints. I'm expecting the answer x1=x2=x3 = 1/3. I tried the "constrOptim" function in R and I'm running into some issues. I first start off

Dear all, I am trying to use SAMr-library(samr), it gives me p-value = 1.001, any idea why? Many thanks, Amor __________________________________________________ z gegen Massenmails. [[alternative HTML version deleted]]

Since the fact that gfortran performance has improved over the major releases, I decided to benchmark the current releases on a MacPro with the Polyhedron 2005 benchmarks using -ffast-math -funroll-loops -msse3 -O3. The results are... gcc release gcc 4.2.4 gcc 4.3.3 gcc 4.4-pre gcc 4.3.3/ gcc 4.4-pre/

Does anyone know if there's a reason that proc.time() uses cat() rather than message() to print the output when there has been an error in the process of timing? line 31 of time.R, https://github.com/wch/r-source/blob/e5b21d0397c607883ff25cca379687b86933d730/src/library/base/R/time.R#L31 on.exit(cat("Timing stopped at:", ppt(proc.time() - time), "\n")) This means that

On Tue, Sep 06, 2005 at 09:56:06PM -0500, Jay Krivanek wrote: > Yea, don't get me wrong, ogg is definitely an awesome format but aacplus v2 > is currently the leader in high fidelity quality at low bitrates. 24kbps > aacplus even kicks ass in my view. Is that 24kbps AACPLUS a full 44.1kHz stereo signal with frequencies preserved up to at least 15kHz? (Just curious, not

I've got a few theora vids here... and looking at the headers it seems to me they say they play at 2 frames per second... but watching them they appear to to play at a faster rate. Theres also some that say they lpay at 1 fps. Am i missing something here ? I'm getting the info from here.... http://wiki.xiph.org/TheoraSpec Is this header structure still valid ? Heres a dump of the BOS

I am trying to simulate the trajectory of the pension assets of one person. In C-like syntax, it looks like this: daily.wage.growth = 1.001 # deterministic contribution.rate = 0.08 # deterministic 8% Wage = 10 # initial Asset = 0 # initial for (10,000 days) { Asset += contribution.rate * Wage

I am using R 1.1 with Redhat 6.2 and RW 1.001 with Win98 (the upkey doesn't work on my IBM either as has been previously reported by others). The function aov doesn't return either the residuals or the residual degrees of freedom for split-plot designs. If you use the following code from Baron and Li's "Notes on the use of R for psycology experiments and questionnaires"

There's something weird going on for certain non-integer values of argument 'timeout' to base::socketSelect(). For such values, there is no timeout and you effectively end up with an infinite timeout. I can reproduce this on R 3.3.1 on Ubuntu 16.04 and RedHat 6.6, but not on Windows (via Linux Wine). # 1. In R master session > con <- socketConnection('localhost', port

Dear all, I tried to use "step" function to do model selection, but I got an error massage. What I don't understand is that data as data.frame worked well for my other programs, how come I cannot make it run this time. Could you please tell me how I can fix it? ***************************************************************************************************