similar to: A function for random test based on longest run (UNCLASSI FIED)

Classification: UNCLASSIFIED Caveats: NONE Hello, R users, Has anybody written a function for random test based on the length of longest run of same events. I really appreciate your help. Kyong Park Classification: UNCLASSIFIED Caveats: NONE [[alternative HTML version deleted]]

Classification: UNCLASSIFIED Caveats: NONE Dear R users, I'm interested in finding a random effect of the Block in the data shown below, but 'lme' does not work without the random effect. I'm not sure how to group the data without continuous value which is shown in the error message at the bottom line. If I use 'aov' with Error(Block), is there a test method comparing

Classification: UNCLASSIFIED Caveats: NONE Dear R users, How can I retrieve an evaluated gradient value from "deriv" function provided below? I want to retrieve b0 value. Appreciate your help. Kyong junk1<-deriv(~(1/b1)*k-b0/b1,"b0",c("k","b0","b1"),formal=T) junk1(k=0,b0=-14.0236,b1=2.44031) [1] 5.746647 attr(, "gradient"):

Hello, R users, I am a very beginner of R and tried read.csv to import an excel file after saving an excel file as csv. But it added alternating rows of fictitious NA values after row number 16. When I applied read.delim, there were trailing several commas at the end of each row after row number 16 instead of NA values. Appreciate your help. Kyong [[alternative HTML version deleted]]

Classification: UNCLASSIFIED Caveats: NONE R users, This is a summary for the responses I got from Sandra Dorai-Raj and Simon Blomberg followed by my question. Sandra suggested using lm instead of lme for a model without a random effect, and Simon suggested anova.lme instead of anova.lm for comparing two models with and without a random effect. For the anova test, both anova.lme and anova

Hi R users, I’d like to get an expected value for a minimum value from order statistics of sample size, say, 5 for standard normal distribution, and the formula would be 5* Integral (from -inf to Inf) x*f(x)* (1-F(x))^4,where f(x) and F(x) are a standard normal density and a cumulative distribution function respectively. After searching R posts, I applied the following double integral

Dear R users, I'm using R 2.1.0 with Windows 2000. In the middle of a R session, an error message pops up saying Rgui.exe generated errors and shuts down the R session. How can I correct this problem? Appreciate your help. Kyong [[alternative HTML version deleted]]

Hello, I am having trouble installing the package from CRAN. As the package compiles I get: * installing *source* package ‘RDieHarder’ ... ** package ‘RDieHarder’ successfully unpacked and MD5 sums checked checking for gcc... gcc checking for C compiler default output file name... a.out checking whether the C compiler works... yes checking whether we are cross compiling... no checking for

Hello, in my rails application I have the following requierement: I want to store the content of all editable fields (form tag) of a page in a session variable when the "link_to" action is selected. (Submit would be easy, for details & background please see below). Thanks a lot Christian I may miss here something, but this is my problem (I tried hard to make it clear - just let

RDieHarder fails its regression tests on x86_64 (aka "amd64") at CRAN (using Debian), and I see the same on Ubuntu 8.10 in 64 bit. No issues on 32bit. One odd thing is that the program behaves well if run via R --no-save < tests/RDieHarder.R but NOT when started using R --slave < tests/RDieHarder.R as R CMD check does. So I tried different things related to startup

I got hold of the 0.0.7 RDieHarder binary package as zip for Windows from CRANextras but don't know how to invoke its main method dieharder() in R with "file_input" as RNG argument - how to specify path to that input file ? the 0.0.7 as opposed to 0.1.0 doesn't seem to have ''inputfile'' arg for path specification. So is "file_input" prior to 0.1.0

Hi, subject more or less says it all. I freely admit to not having bothered to find some of the online papers about method of testing the quality of random number generators -- but in an idle moment I wondered what to expect from something like the following: randa<-runif(1000) randb<-runif(1000) t.test(randa,randb)$p.value var.test(randa,randb)$p.value [repeat ad nauseum] Is the

Hello, R users, I applied segmented regression method contributed by Muggeo and got different slope estimates depending on the initial break points. The results are listed below and I'd like to know what is a reasonable approach handling this kinds of problem. I think applying various initial break points is certainly not a efficient approach. Is there any other methods to deal with segmented

Hi R bayesians, I need an advise how to resolve the two different estimates applying a traditional glm (TG) and a bayes glm (BG), and different results depending on the data formats of response data and the prior specs using bayesglm in R. I'm not familiar with bayes estimate and my colleague asked me to look into this because the EPA from France reported a quite different estimates for

Can someone point me to an implementation of the "runs up and runs down" test, or does such beast exist in R? From web searches the "runs up and runs down" test is commonly used for testing pseudo-random number generators and in simulations. John C. Davis describes its use in geology in his book on geological statistics (2002). Having searched the full site, I can locate only

Dear useRs, A few days ago, the initial version 0.1.0 of a new package 'random' was uploaded to CRAN. The random packages provides convenient access to the non-deterministic random numbers provided by the random.org site created by Mads Haahr (http://www.random.org). While certain hardware and software solutions that provide access to non-deterministic random-numbers exist, few if any

Dear useRs, A few days ago, the initial version 0.1.0 of a new package 'random' was uploaded to CRAN. The random packages provides convenient access to the non-deterministic random numbers provided by the random.org site created by Mads Haahr (http://www.random.org). While certain hardware and software solutions that provide access to non-deterministic random-numbers exist, few if any

I really appreciate all the helpful answers from Richard Müller, Marc Schwartz, Adaikalavan Ramasamy, Vito Ricci. I tried Marc's simple suggestion which was to copy only the data area and paste to a new excel sheet, and that solved my problem. I'll try other suggestions while I'm learning more about R. Unedited answers are provided below followed by my question. Thanks again. Kyong

I appreciate to Roger Koenker, Achim Zeileis and Vito Muggeo for their informative answers. Listed below is unedited replies I got followed by the question I posted. Kyong 1. Roger Koenker: You might try rqss() in the quantreg package. It gives piecewise linear fits for a nonparametric form of median regression using total variation of the derivative of the fitted function as a penalty

Hi. I have a VM with CentOs6.2 on KVM. I mount an image of the VM when it's not working with guestmount: root at diehard:~# guestmount -a named-centos.img -i tmp vi /etc/shadow insert comment line into the file (esc yy p 0 i # esc :wq!), save changes, unmount the image, turn on the VM, and I cannot login. I cannot login even if I recover the shadow file. root at diehard:~# guestmount