similar to: Difference between the S-plus influence and R empinf functions

Hi I'm using boot package for some analysis on linear regression coeficients. My problem is that I can not compute bca intervals, I get an error message > bca.ci(blm8901,index=1) Error in if (!all(rk > 1 & rk < R)) warning("Extreme Order Statistics used as Endpoints") : missing value where logical needed The problem is the empinf.reg function that is

[I thought I'd submitted this bug report some time ago, but it's never showed up on the bug tracking system, so I'm submitting again.] qr.solve() gives incorrect results when dealing with complex matrices or with qr objects that have been computed with LAPACK=TRUE, whenever the b argument has more than one column. This bug flows from qr.coef(), which has a similar problem. I believe

> Hello, > > I'm running R 1.8.1 on a Mac G4 with Dual 800 MHz CPU, 1.5 GB RAM, and > sufficient disk space. > > I have loaded the boot package version 1.2-15. > > I defined a function to pass to the statistic argument which estimates > a vector (length = 6) of parameter estimates; none are variances. > > I defined a function to pass to the ran.gen argument

Full_Name: Jeremiah Cohen Version: 2.9.0 OS: Windows XP Submission from: (NULL) (129.59.230.235) I believe there is a bug in the seq() function for certain values of the "from" argument. Here are examples: > seq(-.2, .1, .1) [1] -0.2 -0.1 0.0 0.1 > seq(-.3, .1, .1) [1] -3.000000e-01 -2.000000e-01 -1.000000e-01 5.551115e-17 1.000000e-01 > seq(-.4, .1, .1) [1] -0.4 -0.3

Hello, I am trying to call a FORTRAN subroutine from R. The Fortran code is @: http://lib.stat.cmu.edu/apstat/206 It performs a bivariate isotonic regression on a rectangular grid (m X n) matrix. I used the g77 compiler and successfully created a dll file and it also loads successfully from R. But somehow the programs fails to run properly. (I do get the correct result when I compile the

Hi everybody, I'm trying to create a numerical data frame on which to perform PRCC. So far I have created a data frame that consists of function/vector output that displays in numerical form, but when I try and run PRCC (from epiR package) I get the following error message: "Error in solve.default(C) : Lapack routine dgesv: system is exactly singular" It appears this is because

On Sep 29, 2010, at 8:35 AMPDT, Ralf Karrenberg wrote: > Hello everybody, > > I have stumbled upon a test case (the attached module is a slightly > reduced version) that shows extremely reduced performance on linux > compared to windows when executed using LLVM's JIT. > > We narrowed the problem down to the actual code being generated, the > source IR on both systems

Hi everyone, I emailed yesterday thinking I had a problem with reading matrix information as numerical output into epiR. After working on the data frame today, I realise it's a problem with the data itself rather than the syntax. An original data frame of min max T1 1.500000e+01 3.999954e+01 SE1 0.000000e+00 1.000000e+00 PRE

Hello everyone, I am using the following program to get the p-value of some numbers (column 'LR' of the data.dat file). I want to write the 1st and 2nd column of the output file (data.out) as an integer while the program change them. Could anybody please tell me how I can write the code which writes the values of the first two columns as integer? For your convenience, I have attached the

Hi everyone, I have a simple problem but don´t know how to solve it. I read in a table from a text file that looks like that: -3,932200E-01 -2,000000E-01 4,999995E-02 -3,932099E-01 -1,000000E-01 3,999996E-02 -3,932000E-01 0,000000E+00 3,999996E-02 -3,931899E-01 -1,000000E-01 4,499996E-02 -3,931800E-01 -1,000000E-01 4,499996E-02 -3,931699E-01

Hello! I have two matrices of equal dimension ll and lu that I want to do a ks.test on corresponding rows in these matrices (dim(ll) is [1] 101 100). If I do e.g. > ks.test(ll[50,],lu[50,]) just for testing, it displays a lot of numbers, and some more info: [116] -3.000000e-02 -4.000000e-02 -4.000000e-02 -3.000000e-02 -4.000000e-02 [121] -3.000000e-02 -4.000000e-02 -3.000000e-02 -2.000000e-02

Dear R users, I am using 'R' version 2.2.1 and 'S-PLUS' version 6.0; and I loaded the MASS library in 'S-PLUS'. I am running a logistic regression using glm: --------------------------------------------------------------------------- > mydata.glm<-glm(COMU~MeanPycUpT+MeanPycUpS, family=binomial, data=mydata)

Hi, I have a following vector: > smallch [1] 0.0652840 0.1181300 0.0319370 0.0155700 0.0464110 0.0107850 [7] 0.0158970 0.0375900 0.0603090 0.0310300 0.0105920 0.0540580 [13] -0.0177740 0.0039393 Pretty (R 1.5.1) has problems with zero: > pretty(smallch) [1] -2.000000e-02 -3.469447e-18 2.000000e-02 4.000000e-02 6.000000e-02 [6] 8.000000e-02 1.000000e-01 1.200000e-01

Hi, I have a following vector: > smallch [1] 0.0652840 0.1181300 0.0319370 0.0155700 0.0464110 0.0107850 [7] 0.0158970 0.0375900 0.0603090 0.0310300 0.0105920 0.0540580 [13] -0.0177740 0.0039393 Pretty (R 1.5.1) has problems with zero: > pretty(smallch) [1] -2.000000e-02 -3.469447e-18 2.000000e-02 4.000000e-02 6.000000e-02 [6] 8.000000e-02 1.000000e-01 1.200000e-01

Leandra Desousa <sousa <at> ims.uaf.edu> writes: >> > I am using 'R' version 2.2.1 and 'S-PLUS' version 6.0; and I loaded the >> > MASS library in 'S-PLUS'. >> > >> > I am running a logistic regression using glm: >> > >> > >summary(mydata.glm) >> > Call: >> > glm(formula = COMU ~

I have the following output from pretty(): > pretty(c(-.1, 1)) # note 2nd element [1] -2.000000e-01 -2.775558e-17 2.000000e-01 4.000000e-01 6.000000e-01 [6] 8.000000e-01 1.000000e+00 > as.character(pretty(c(-.1, 1))) [1] "-0.2" "-2.77555756156289e-17" "0.2" [4] "0.4" "0.6"

I need a little piece of advice concerning passing data frames into the functions. As I do a lot of similar fits at a time, I'd like to write a small function doing the fits for all relevant variables automatically. However, I usually get error messages of the following kind: (I present here a part of a test code). ##################################################### # Data set:

I've attached the module->dump() that our code is producing. Unfortunately this is the smallest test case I have available. This is before any optimization passes are applied. There are two separate modules in existence at the time, and there are no guarantees about the order the surrounding code calls those functions, so there may be some interaction between them? There shouldn't

Can anyone explain this? I have a matrix with double components. It's taking up a lot of memory, so I want to multiply then turn it to integers. I'm pretty certain that there are only 2 decimal places, but I wanted to check by using modulo. E.g. mat = matrix(11:50/100, ncol=4,nrow=10) #Matrix with values out to the hundredths any((mat * 100)%%1!=0) But oddly enough it doesn't work.

Dear R users, I'm trying to construct a distance matrix based on a nb object created by spdep where sites must have a larger influence in one direction then in the other. Here is an example to better illustrate what I need: Let say I have the following Gabriel connection network library(spdep) library(ade4) data(orbatid) nbgab<-graph2nb(gabrielneigh(as.matrix(oribatid$xy)))