similar to: Error: cannot set length of non-vector

The functions metrop and temper in the mcmc package have a debug = FALSE argument that when TRUE adds a lot of debugging information to the returned list. This is absolutely necessary to test the functions, because one generally knows nothing about the simulated distribution except what what one learns from MCMC samples. Hence you must expose all details of the simulation to have any hope of

Full_Name: Charles J. Geyer Version: 1.8.0 OS: i686-pc-linux-gnu (Suse 8.2) Submission from: (NULL) (134.84.86.22) Two bugs (perhaps related, perhaps independent) revealed by the same Poisson regression with offset mydata <- read.table(url("http://www.stat.umn.edu/geyer/5931/mle/seeds.txt")) out.fubar <- glm(seedlings ~ burn01 + vegtype * burn02 + offset(log(totalseeds)),

Just built 3.6.3 from source and tcl doesn't work. Worked fine with the same laptop in 3.6.2. Here's the exact error. blurfle$ R --vanilla R version 3.6.3 (2020-02-29) -- "Holding the Windsock" Copyright (C) 2020 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to

OK. I give up. I'll ask a stupid question. How do I get the $!#@*$ signif stars line printed by summaries to not look extremely bizarre in the latex produced by Sweave? For example, see p. 7 of http://www.stat.umn.edu/geyer/aster/library/aster/doc/tutor.pdf I can see what the problem is. R emits non-ascii characters (as it is supposed to do), Sweave puts them in the tex file, and

I knew I could work around. But this shouldn't happen. And yes. Same problem with your example. blurfle$ R --vanilla R version 3.6.3 (2020-02-29) -- "Holding the Windsock" Copyright (C) 2020 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain

>>>>> Charles Geyer writes: > In that case, perhaps the question could be changed to could CC0 be > added to the list of R licences. Right now the only CC licence that > is in the R licenses is CC-BY-SA-4.0. Hmm, I see Name: CC0 FSF: free_and_GPLv3_compatible (https://www.gnu.org/licenses/license-list.html#CC0) OSI: NA (https://opensource.org/faq#cc-zero) URL:

is it possible to set an option inside a function ((I want to set na.action = na.fail) and have the previous state restored if there is an error so that the function doesn't change the option behind the user's back? Sorry if this has been answered before, but this subject is hard to Google. -- Charles Geyer Professor, School of Statistics University of Minnesota charlie at stat.umn.edu

I would like the unlicense (http://unlicense.org/) added to R licenses. Does anyone else think that worthwhile? -- Charles Geyer Professor, School of Statistics Resident Fellow, Minnesota Center for Philosophy of Science University of Minnesota charlie at stat.umn.edu

I have hit a problem with the design of the mcmc package I can't figure out, possibly because I don't really understand the R function call mechanism. The function metrop in the mcmc package has a ... argument that it passes to one or two user-supplied functions, which are other arguments to metrop. When the two functions don't have the same arguments, this doesn't work.

I have a design problem. I'm trying to design a package that does something like glm, but 1. the response is multivariate, so I can't be just like glm and get the response variables out of the formula. I decided (perhaps incorrectly) to just supply the response variable names as an argument "response". 2. I have the usual predictor variables. 3. I

Hello, klibc-1.4 and klibc-1.4.29 don't export the creat function: $ klcc -static -s -Wall rtfs.c -o rtfs rtfs.c: In function 'move_ent': rtfs.c:318: warning: implicit declaration of function 'creat' rtfs.o: In function `move_ent': rtfs.c:318: undefined reference to `creat' $ grep creat\\b /usr/lib/klibc/include/ -r /usr/lib/klibc/include/zlib.h: descriptors are

Can anyone tell me if the following C code (which proved very useful when using the snow package -- use it to nice slaves) compiles and dyn.loads under Windoze or Mac? It is (apparently) POSIX, so I suppose it is fine in OS X, but does Windoze's advertized POSIX compliance mean anything here? If it doesn't work on some platform, how do I deal with that? I want to submit as contributed

The note to 1004 says "fixed for 1.3.1" Uh. No. It ain't. The problem was more serious than guessed as even the simplest testing would show. For example, Example 5.4 in Hollander and Wolfe (Nonparametric Statistical, Methods, 2nd ed., Wiley, 1999, pp. 180-181) R Version 1.3.1 (SuSE Linux 7.1) > X <-

Sorry if I was supposed to file a bug report, but I don't know whether we're supposed to do that on R devel. I just built R devel from source (2010-06-15 r52280) and tried to check a package with it and Sweave failed on the vignette. It puts the line \usepackage{/HOME/faculty/charlie/local/devel/lib64/R/share/texmf/Sweave} in the *.tex file but, that's not where Sweave.sty is oak$

R CMD check --use-valgrind <packagename> used to run valgrind on the tests in the tests directory of the package. But it seems to have stopped. R-2.9.1 doesn't -- at least on my box -- and neither does R-2.10.0 (devel). I am not sure when this stopped. I think 2.8.x did this. The only old R I have around is 2.6.0 and it certainly does. R CMD check --help for 2.9.1 says (among other

Dear All, Apologies for not providing a reproducible example, but if I could, then I would be able to answer myself my question. Essentially, I am trying to fit a very complicated custom probability distribution to some data. Fitdistrplus does in principle everything which I need, but if require me to specify not only the density function d, but also the cumulative p and and inverse cumulative

There is a fairly new bug in wilcox.test in R-2.2.1 (stable). It wasn't there when I last taught nonparametrics in fall 2003. Line 86 of wilcox.test.R achieved.alpha<-2*psignrank(trunc(qu),n) It should be achieved.alpha<-2*psignrank(trunc(qu)-1,n) If you don't see why, decode the cookbook instructions p. 56 in Hollander and Wolfe (2nd ed.) or see

I think I am using objects according to the man page. This seems to be a valid regular expression. But whether I know what I'm doing or no, it still shouldn't be doing what valgrind seems to be saying it's doing. (IMHO) ---------- start of script ---------- Script started on Mon 01 Aug 2005 02:09:00 PM PDT linux$ printenv VALGRIND_OPTS --tool=3Dmemcheck linux$ cat bar.R foo <- 1

So R is on subversion right? So why doesn't somebody authoritative get r-project.org on the subversion testimonials page? http://subversion.tigris.org/testimonials.html Just a suggestion. -- Charles Geyer Professor, School of Statistics University of Minnesota charlie@stat.umn.edu

I am puzzled at the use of regression. I have a categorical variable ClassePop33000 which factors a Population variable into 3 levels. I want to investigate whether that categorical variable has some relation with my dependent variable, so I go : lm(Cout.ton ~ ClassePop33000, data=ech2) Call: lm(formula = Cout.ton ~ ClassePop33000, data = ech2) Residuals: Min 1Q Median 3Q