similar to: ansari.test.default: bug in call to uniroot?

Full_Name: Wei Xu Version: 1.5.1 OS: WindowsME Submission from: (NULL) (63.215.238.92) The P-value for a two.sided test is not consistent with the confidence interval. For example, P-value=0.1372, but the 95% CI doesn't include the H0 value(1). > x [1] 0.80 0.83 1.89 1.04 1.45 1.38 1.91 1.64 0.73 1.46 > y [1] 1.15 0.88 0.90 0.74 1.21 >

Full_Name: Gael Millot Version: 2.2.0. OS: XP Submission from: (NULL) (195.220.102.20) Hello. I sent an Email in r-help without answer for the moment. I am wondering if it could have a mistake in the code of the ansari.test function. For me, it seems that the function do not recover the p value at the correct side of the normal law N(0, 1) when it use the normal approximation (presence of

The help page for ansari.test() says (emphasis added): By default (if exact is not specified), an exact p-value is computed if both samples contain less than 50 finite values **and there are no ties**. Otherwise, a normal approximation is used. However, this does not appear to be the case in R 2.4.1 or R 2.5.0. In fact, even example(ansari.test) produces a warning Warning message:

hello, I have the following problem setting parameter 'frac' in the sigest function of the kernlab package. ## executing the ?sigest example: library(kernlab) data(spam) srange <- sigest(type~.,data = spam) ## works fine... ## setting 'frac' explicitly ## (in this case even to the default of .25) options(error=recover) srange <- sigest(type~.,data = spam, frac = .25) ##

Hi, I created the model below, which returns me the following warning message: In sem.default(ram = ram, S = S, N = N, param.names = pars, var.names = vars, : Could not compute QR decomposition of Hessian. Optimization probably did not converge. ######### Model ######## mDPDF = data.frame(mj1,mj2,mj3,mj4,mj5,eL1,eL2,eL3,eL4,eL5,aC1,aC2,aC3,aC4,disR1,disR2,disR3,disR4,disR5,

Hi all, This is a statistics question, I hope someone out there will be able to help me. I have one population (oligonucleotide probes spotted on a Nimblegen array). I measured "parameter one" (intensity after hybridisation) and I have selected a subpopulation of the initial (one tenth of the initial) according to a threshold value. I now measure "parameter two" of

Hello. I am probably wrong... I am wondering if it could have a mistake in the code of the ansari.test function. For me, it seems that the function do not recover the p value at the correct side of the normal law N(0, 1) when it use the normal approximation (presence of ties) in a one tailed test. Exemple : quanti<-c(197, 205, 228, 234, 237, 195, 233, 226, 244, 227, 259, 185, 198, 253, 207,

Dear list, I have a long list of two vectors with some matching elements. I would like to line them up in two columns and have an NA in those positions of the second vector where a match is absent. With a simple example, I will explain my problem. (a<-1:6) (b<-c(5,2)) (m1<-match(a,b)) (ab<-cbind(a,m1)) m2<-numeric(length(m1)) for (i in 1:length(m1))

I am trying to get R 1.0 running on the Mac. The main target is MacX. Anyone else working on that? For the recent Macintosh system, I am trying to compile R using MachTen. The R core compiles and runs without any changes, using the Unix make files. The libraries give some problems, presumably due to handling of shared libraries. It seems I am missing something to link. Before I spend my time

I would like to be able to use two functions; qansari and pansari which are found in the function ansari.test. How can I evaluate these functions independently? Thanks. --Dale For example, when I load the function ... qansari <- function(p, m, n) { .C(R_qansari, as.integer(length(p)), q = as.double(p), as.integer(m), as.integer(n))$q } and

I have two datasets: A with columns Open and Name (and many others, irrelevant to the merge) B with columns Time and Name (and many others, irrelevant to the merge) I want the dataset AB with all these columns Open from A - a difftime (time of day) Time from B - a difftime (time of day) Name (same in A & B) - a factor, does NOT index rows, i.e., there are _many_ rows in both A & B with

Hi all, I've analyzed the packet stream resulting from a SIP session through Asterisk, and I found that Asterisk sends a lot of useless messages! In example, it repeats the three way handshake (invite-ok-ack) every time, also when it has just to forward the 'bye' message. Has anyone found the same srange thing? thanks..

Hi all, Got a srange problem and wondered if anyone could shed some light on it. Have set up Samba but I cant log in from a windows workstation (get a message saying that either the password is incorrect or access was denied). I have set all the users up and Linux users (red hat 7.2) and as smb user (had to do this manually though as smbpasswd -a didnt copy them across) and both users accounts

I am writing a package, which contains several functions and variables. The variables are for internal use by functions. Functions are ment to be callable by a user. However, the function cannot access these variables. The package contains only R code, and was created using package.skeleton package.skeleton(name="ac9", list=c(ls(pattern="(AB2)|(ac9)"),

Hi, I'm trying to compile R-1.7.1 from source (on a RedHat 8.0) instead of using the binary version, as it has often been advised. However I don't manage to find a solution to the following error which occurs during the make procedure : /!\-------------------------------------------/!\ building package 'ctest' mkdir -p -- ../../../library/ctest/R mkdir -p --

Hi all, Trying to install R-1.7.1 on a RedHat 8.0 platform, I have a few problems. R gets installed without default packages (but base and ctest) : make script fails at the end of the procedure (configure is made successfully). This is the (translated) log I have : /gcc -I../../../../include /usr/include/mysql -D__NO_MATH_INLINES -mieee-fp -fPIC -g -O2 -c ansari.c -o ansari.o gcc: cannot

Hi, I'm working on the codes below however every time I run them when they get to OpenBUGS I keep getting the error message: array index is greater than array upper bound for hab. Any help would be greatly appreciated, Suzie Codes: ungulate <- read.csv(file.choose ()) #ungulate ungulate <- as.matrix(ungulate);colnames(ungulate)<-NULL;rownames(ungulate)<-NULL

Hi, I tried compiling the latest pre-release for R 3.3.3 for the NixOS Linux distribution [1], but the build fails during the "make check" phase because of the following 2 issues: 1) The "tools" test in "tests/Examples" requires network access, which it doesn't have in our build environment. Therefore, it fails as follows according to

Sorry to all that are angry about the form of my previous mail. I didn't realise what would happen :((. Here it is in (hopefully) plain text (if my mailer doesn't spoil it again): ############## Dear developers, I have a problem with some discrepancy between R 1.2.1 for Windows and R 1.2.2 (and less) for Linux. While trying to correct the wilcox.test (see my previous bug report) I

I have a data.frame object like: > data.frame(x=c('A','A','B','B'), y=c('Ab','Ac','Ba','Bd')) x y 1 A Ab 2 A Ac 3 B Ba 4 B Bd how could I create a tree structure object like this: |---Ab A---| _| |---Ac | | |---Ba B---| |---Bb Thanks, Zech [[alternative HTML version deleted]]