similar to: Improvement of [dpq]wilcox functions

Full_Name: Mark J. Lamias Version: 1.7.0 OS: Windows 2000 Pro Submission from: (NULL) (65.222.84.72) I am running the qwilcox procedure and it is producing incorrect results. For example, dwilcox(.025, 3, 5) should equal 6, but it is equal to 1. Similarly, dwilcox(.025, 3, 6) should equal 7, but it equals 2. The critical values are not set being returned with the correct values. I've

I?ve been looking at this for a couple hours--it ended up being trickier than I expected to implement well. I?ve attached a new patch here. This version scales significantly better than the existing method for both `pwilcox` and `dwilcox`. Examples are included below. I can?t think of any other ways to improve runtime at this point. That?s not to say there aren?t any, though?I?m hopeful

Has there been a recent change in the behavior of binary operators? In SparseM it was, until quite recently ok to do scalar multiplication but now, > A*4 Error in A * 4 : non-numeric argument to binary operator > 4*A Error in 4 * A : non-numeric argument to binary operator > A%*%A An object of class "matrix.csr" Slot "ra": [1] 1 1 1 Slot "ja": [1] 1 2 3

I took some time and liberty and tried to improve existing implementation of SignRank functions in R. (dsignrank, ...) As I have seen they've been based on csignrank. So I modified csignrank and, I believe, improved calculation time and memory efficiency. The idea is basically the same. I use the same recursion as original author used with one slight modification. I am generating Wilcoxon

Hello, can someone please help me with coding a function for a bootstrap WMW test (package boot, R under Windows, version 1.6.2)?

Hi again ! Thanks for the info on updating the config.site file which I have done. I have also added -lm in the Makeconf manually because this is needed explicitly for DEC cc. However, there are still a few problems when linking some of the files as you can see from the enclosed log. Ciao, Andreas ------------------------------------------------------- R-0.63.3>make make[1]: Entering

Hi, Just downloaded the R-0.62.4 of R and tried to install it. With the standard procedure : ./configure make At the end of the compilation I got the following message : ld: ../lib/libunix.a(system.o): main: multiply defined fort: Severe: Failed while trying to link. *** Exit 1 Stop. *** Exit 1 Stop. *** Exit 1 Stop. *** Exit 1 Stop. I attach the printout after the ./configure and make

>>>>> akshay kulkarni >>>>> on Wed, 31 May 2023 20:55:33 +0000 writes: > dear members, > I am using arfima() from forecast package to model a time > series. The following is the code: >> LYGH[[202]] > [1] 45.40 3.25 6.50 2.15 >> arfima(LYGH[[202]]) > Error in .fdcov(x, fdf$d, h, nar = nar, nma = nma,

Dear Martin, Sad that the bug is beyond your ken... Fortunately, the error happens only rarely...The length of LYGH was 719 and there were only two such errors..I will just replace them with NA and make do. By the by, what if I send LYGH as an attachment to your actual mail ( not the r-help mail)? Will it help? Can you then pinpoint the cause? Or should I raise a bug

hi I'm puzzled as to the relation between the REML score computed by gam and the formula (4) on p.4 here: http://opus.bath.ac.uk/22707/1/Wood_JRSSB_2011_73_1_3.pdf I'm ok with this for poisson, or for quasipoisson when phi=1. However, when phi differs from 1, I'm stuck. #simulate some data library(mgcv) set.seed(1) x1<-runif(500) x2<-rnorm(500)

> > > Performance statistics are interesting. If we assume the two populations > have a total of `m` members, then this implementation runs slightly slower > for m < 20, and much slower for 50 < m < 100. However, this implementation > works significantly *faster* for m > 200. The breakpoint is precisely when > each population has a size of 50; `qwilcox(0.5,50,50)`

Hello, I've allways had some problems compiling R (R-62.x and R-63.x) on a AIX 4.2.1 System to the end: make produces following output: ... ... mkdir ../../../library/eda/R ld: 0711-327 WARNING: Entry point not found: __start ld: 0711-244 ERROR: No csects or exported symbols have been saved. gnumake[4]: *** [eda.so] Error 8 gnumake[3]: *** [all] Error 1 gnumake[2]: *** [R] Error 1 gnumake[1]:

Hello, I've allways had some problems compiling R (R-62.x and R-63.x) on a AIX 4.2.1 System to the end: make produces following output: ... ... mkdir ../../../library/eda/R ld: 0711-327 WARNING: Entry point not found: __start ld: 0711-244 ERROR: No csects or exported symbols have been saved. gnumake[4]: *** [eda.so] Error 8 gnumake[3]: *** [all] Error 1 gnumake[2]: *** [R] Error 1 gnumake[1]:

dear members, I am using arfima() from forecast package to model a time series. The following is the code: > LYGH[[202]] [1] 45.40 3.25 6.50 2.15 > arfima(LYGH[[202]]) Error in .fdcov(x, fdf$d, h, nar = nar, nma = nma, hess = hess, fdf.work = fdf$w) : NA/NaN/Inf in foreign function call (arg 5) I tried viewing .fdcov() with the following code:

This is a multipart MIME message. --==_Exmh_981436288450 Content-Type: text/plain; charset=iso-8859-1 Content-Transfer-Encoding: quoted-printable Hi ! I have problems compiling R successfully on a DEC-UINX 4.0E. I have applied the recommended config.site, which I enclose. As can be seen from the compilation log there are linking errors... I did a 'make check' which fails for the R

I've discovered an infelicity (I guess) in qbeta(): it's not a bug, since there's a clear warning about lack of convergence of the numerical algorithm ("full precision may not have been achieved"). I can work around this, but I'm curious why it happens and whether there's a better workaround -- it doesn't seem to be in a particularly extreme corner of parameter

Can anyone help? I am trying to install R-0.62.2 on a Sun with OS 5.6 (Solaris 2.6), and gcc/f77. Unfortunatley the simple commands ./configure --prefix=/opt/R make did not work, complaining when building help/docs. It claims to have problems perl (which has worked fine so far) but in practice it is not only a documentation problem: if I force "make install" and run R, then it

Hi, I'm using the function fracdiff and can not figure out how to get the estimated values for sigma2 or confidence intervals for the parameter estimates. Does anyone know how to obtain these values? Thanks, Melissa

Full_Name: M. Welinder Version: 1.8.1 OS: Solaris Submission from: (NULL) (65.213.85.227) Currently R has... #define R_D_Lval(p) (lower_tail ? (p) : (1 - (p))) /* p */ #define R_D_val(x) (log_p ? log(x) : (x)) /* x in pF(x,..) */ #define R_DT_val(x) R_D_val(R_D_Lval(x)) /* x in pF */ ...which is sub-optimal in the lower_tail==FALSE && log_p==TRUE case. Something like this ought

I am a newbie to C (although I did some programming in Perl and Pascal) so forgive me for the language that follows. I am writing because I have a question concerning the *implementation of *(the internal function)* cwilcox* in https://svn.r-project.org/R/!svn/bc/81274/branches/R-DL/src/nmath/wilcox.c. This function is used to determine the test statistics of the Wilcoxon-Mann-Whitney U-test.