search for: shengqiao

How to use sub, gsub, etc. to replace "\" in a string to "/"? For example,convert "C:\foo\bar" to "C:/foo/bar". Thanks, Shengqiao Li

...erent. I have many predictors (p>20) and I really want try knn with a given k. ksmooth and loess use band width to define neighborhood size. This contrasts to knn's variable band width via fixing a k. Are there any such functions I can use in R packages? Your help is highly appreciated. Shengqiao Li

...ed in 10M numbers generated by runif using Wichmann-Hill. The other three methods had cycle length of 2^32. So, anybody can explain this? And any improvement to the implementation can be made to increase the cycle length like the Wichmann-Hill method? ======================================== Shengqiao Li Research Associate The Department of Statistics West Virginia University Morgantown, WV 26506-6330

Unfortunately, RNGkind is buggy. It will not generate warnings except the full name "Buggy Kinderman-Ramage" is supplied for normal.kind. match.arg is supposed to be called before "==" comparison. ======================================== Shengqiao Li Research Associate The Department of Statistics PO Box 6330 West Virginia University Morgantown, WV 26506-6330

...ey should be close. So how to=20 explain above problem? I need generate a large sample without any ties, it seems to me=20 "Wichmann-Hill" is only choice right now. =3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D= =3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D Shengqiao Li The Department of Statistics PO Box 6330 West Virginia University Morgantown, WV 26506-6330 =3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D= =3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D On Thu, 14 Aug 2008, Peter Dalgaard wrote: > Shengqiao Li wrote: >> H...

...hmann's correction in 1984). Thus M <= 6.9536e12. D >= 7.19052 when N=1e7 and D>= 179.763 for N=5e7. But in my tests, duplications were surprisingly not observed. It seems that Wichmann-Hill has a much larger cycle than the one documented! Anybody can solve this puzzle? Regards, Shengqiao Li Department of Statistics West Virgina Unversity ==============Testing=================== > RNGkind(kind="Knuth-TAOCP"); > RNGkind(); [1] "Knuth-TAOCP" "Inversion" > sum(duplicated(runif(1e7))); [1] 46216 > > RNGkind(kind="Knuth-TAOCP-2002&quot...

...r computer. I'm using a supercomputer which require to submit jobs to PBS queue for dispatching. I tried use mpirun in my PBS script. But all my Rslaves are spawned to the same node. This is not desired. Any suggestions are welcome! Thanks in advance. ======================================== Shengqiao Li Research Associate The Department of Statistics PO Box 6330 West Virginia University Morgantown, WV 26506-6330 Phone: (304) 285-5960

Please remove me from the mailing list -----Original Message----- From: r-help-bounces at stat.math.ethz.ch To: Shengqiao Li Cc: r-help Sent: 24/03/2004 10:10 AM Subject: Re: [R] Rmpi and PBS On Tue, 2003-12-30 at 08:39, Shengqiao Li wrote: > Hello: > > Anybody knows how to run Rmpi through PBS (Portable Batch System) on a > cluster computer. I'm using a supercomputer which require to submit jobs &gt...

Shengqiao Li wrote: > Hello all, > > I am generating large samples of random numbers. The RNG help page says: > "All the supplied uniform generators return 32-bit integer values that are > converted to doubles, so they take at most 2^32 distinct values and long > runs will return du...

Hi, R guys: I'm using L-BFGS-B method of optim for minimization problem. My function called besselI function which need non-negative parameter and the besselI will overflow if the parameter is too large. So I set the constraint box which is reasonable for my problem. But the point outside the box was test, and I got error. My program and the error follows. This program depends on CircStats

Shengqiao Li wrote: > Hello all, > > I am generating large samples of random numbers. The RNG help page > says: "All the supplied uniform generators return 32-bit integer > values that are converted to doubles, so they take at most 2^32 > distinct values and long runs will return du...

CRAN (and crantastic) updates this week New packages ------------ * disclapmix (0.1) Maintainer: Mikkel Meyer Andersen Author(s): Mikkel Meyer Andersen and Poul Svante Eriksen License: GPL-2 http://crantastic.org/packages/disclapmix disclapmix makes inference in a mixture of Discrete Laplace distributions using the EM algorithm. * EstSimPDMP (1.1) Maintainer: Unknown Author(s):