[R] lmodel2 - regression through origin?

Hi all, is regression through the origin possible at all using lmodel2? I tried the standard technique of inserting a minus 1 into the formula to get regression through the origin in R, but the output still gives me an intercept? For comparison, this is me trying to do regression through the o

[R] MH algorithm syntax help

Hi all, I am trying to write a random walk metropolis hastings algorithm, I using the latest debian distribution of R. Can anyone tell me what I need to insert in the below code? I have tried putting various combinations of curly brackets and punctuation between the "acc=1" and "if" statements,

Re: [R] MH algorithm syntax help

Well spotted, the b is a muck up (what happens when you are basing stuff on someone else's code). The bit you though might be a dimension mismatch seems to work ok, but the bit that I was worried about from the start doesn't; the reason I have sqrt of var y/ var x is because my posterior is for a

[R] Gibbs sampler...did it work?

I am writing a Gibbs sampler. I think it is outputting some of what I want, in that I am getting vector of several thousand values (but not 10,000) in a txt file at the end. My question is, is the error message (see below) telling me that it can't output 10,000 values (draws) because of a limitat

[R] logistic regression - unequal groups in R

I am getting a repeated error when I try to run a logistic regression in R 2.8.1 >(glm(prop1~x1,data=glm1,family=binomial("logit"),weights=nt1)) Error in model.frame.default(formula = prop1 ~ x1, data = glm1, weights = nt1, : invalid type (list) for variable 'x1' x1 is multistate categorical

[R] R - need more memory, or rejection sampling algorithm doesn't work?

Hi all, I am trying to run rejection sampling for the quantity z11 in the function below. Unfortunately I can't simplify the function further so that z11 only appears once. Whenever I run the algorithm, R looks as if it is running it (no error messages or anything), but then nothing happens fo

Re: [R] R - need more memory, or rejection sampling algorithm doesn't work?

riable; z11 is supposed to be a probability of success in a bernoulli trial, so it will range between 0 and 1 (0.5 * 2). Ned Daniel Nordlund-2 wrote: > >> -Original Message- >> From: r-help-boun...@r-project.org >> [mailto:r-help-boun...@r-project.org] On Behalf

[R] bayesm package not downloading via any mirror or repository

I am a pretty new R user, I am running the latest linux version on xandros, updated with some extra debian packages, and I also run the latest windows version, but prefer linux. I am having trouble downloading "bayesm", it won't do it all from any of the sites on the web, I resorted to this one,

Re: [R] bayesm package not downloading via any mirror or repository

e me I've tried many different ones!) --- Please select a CRAN mirror for use in this session --- Loading Tcl/Tk interface ... done Warning in download.packages(unique(pkgs), destdir = tmpd, available = available, : no package 'bayesm' at the repositories E Uwe Ligge

[R] X11 not being detected whilst installing "rgl" package

Any ideas? I certainly have X11, I have tried copying it into other directories. At present it is in /etc/X11. This is my error message. /home/user> R CMD INSTALL -l [/usr/local/R/site-library] rgl_0.81.orig.tar.gz * Installing *source* package 'rgl' ... checking for gcc... gcc -std=gnu99 checki