similar to: problem with mfg argument of par

I'm using scan in a script to record a series of responses of the user as a function of some graphs that I put up on the screen. A toy version would be, y <- rep(NA, 3) for (ix in seq( length(y) ) ) { y[ix] <- scan( n = 1 ) } However, if I include any code after this loop, for example, y <- rep(NA, 3) for (ix in seq( length(y) ) ) { y[ix] <- scan( n = 1 ) } y I get an error

New packages psyphy and MLDS are available on CRAN: psyphy ncludes an assortment of functions useful in analyzing data from pyschophysical experiments. It includes functions for calculating d' from several different experimental designs, links for mafc to be used with the binomial family in glm (and possibly other contexts) and selfStart functions for estimating gamma values

New packages psyphy and MLDS are available on CRAN: psyphy ncludes an assortment of functions useful in analyzing data from pyschophysical experiments. It includes functions for calculating d' from several different experimental designs, links for mafc to be used with the binomial family in glm (and possibly other contexts) and selfStart functions for estimating gamma values

Hi, When I compile the construction \begin{Scode}{eval=FALSE} ?HSP \end{Scode} with Sweave and latex, it outputs in the pdf as, > `?` (HSP) which is not incorrect but a bit more formal than I wanted for demonstrating the use of the help shortcut. I would like the output to look like, > ?HSP but I can't seem to make this work. I have also tried the results=verbatim argument.

Hi, I'm fitting a standard nonlinear model to the luminances measured from the red, green and blue guns of a TV display, using nls. The call is: dd.nls <- nls(Lum ~ Blev + beta[Gun] * GL^gamm, data = dd, start = st) where st was initally estimated using optim() st $Blev [1] -0.06551802 $beta [1] 1.509686e-05 4.555250e-05 7.322720e-06 $gamm [1] 2.511870 This works fine but I

Must packages in the Depends field of the DESCRIPTION file have NAMESPACES? I haven't seen this explicitly indicated anywhere. I am writing a small package and find that when I add the abind package to the list of the Depends field, I get an error in R CMD check of the build. * checking package name space information ... OK * checking package dependencies ... ERROR Packages required but not

Hi, I'm trying to redefine the contrasts for a linear model. With a 2 level factor, x, with levels A and B, a two level factor outputs A and B - A from an lm fit, say lm(y ~ x). I would like to set the contrasts so that the coefficients output are -0.5 (A + B) and B - A, but I can't get the sign correct for the first coefficient (Intercept). Here is a toy example, set.seed(12161952) y

Hi, I was taken off guard by the following behavior in a lattice plot. I frequently want to add a predicted curve defined at more points than in the formula expression of xyplot. There have been numerous examples of how to do this on r-help, but I still often struggle to make this work. I just realized that specifying one of the axes on a log scale does not guarantee that the added data for a

i was searching for the source of bca.ci, a function of the package boot. I tried require(boot, keep.source=TRUE) but again the source was not viewable. How should i do? Best regards Meinhard Ploner ---------------- PS > version _ platform i386-apple-darwin8.8.1 arch i386 os darwin8.8.1 system i386, darwin8.8.1 status major 2

Hi I just installed the gmodels package and the installation was successful but when I was trying to load the library I got an error (see below). Interesting, yesterday I wrote to the maintainer of RSQLite apckage because I got the same error. Does somebody knows what is going on ?? thanks, Mayte library(gmodels) Error in dyn.load(x, as.logical(local), as.logical(now)) : unable to load

Hi I just installed the gmodels package and the installation was successful but when I was trying to load the library I got an error (see below). Interesting, yesterday I wrote to the maintainer of RSQLite apckage because I got the same error. Does somebody knows what is going on ?? thanks, Mayte library(gmodels) Error in dyn.load(x, as.logical(local), as.logical(now)) : unable to load

Dear all, I would like to generate 25 numbers from 1 to 100 but I would like to have some numbers that could be more probable to come out. I was thinking of the function runif: runif(25, 1, 100) , but I don?t know how to give more weight to some numbers. Example: each number from 2 to 10 has the probability of 40% to come out but the probability of each number from 11 to 100 to come out is

Hi, there: After I upgraded my R to 2.4.1, it is my first time of trying to use MASS and found the following error message: > install.packages("MASS") --- Please select a CRAN mirror for use in this session --- trying URL 'http://cran.cnr.Berkeley.edu/bin/macosx/universal/contrib/2.4/VR_7.2-33.tgz' Content type 'application/x-gzip' length 995260 bytes opened URL

I have a list of groups of xy positions I want to set to 0 in an array full of 1s. When the assignments are done directly from the list, they are incorrect, while if I use a temporary array derived from the list the assignments are correct. The following example will hopefully make my problem clearer. The matrices z and zz are initialised with 1. The z and zz values at coordinates in

Shuji, I suspect that bruto blows up because your data are linearly separable. To see this (if you didn't already know), try library(lattice) splom(~x, groups = y) and look at the first row. If you are trying to do classification, there are a few methods that would choke on this (logistic regression) and a few that won't (trees, svms etc). I would guess that bruto is in the latter

Hi, I'm raising the issue, again, of those pesky thin lines that show up with a pdf graphics device, otherwise known as FAQ 7.36. I have read the FAQ and understand that I can eliminate them by adjusting the anti-aliasing option for viewing them on the screen, but this doesn't necessarily eliminate them for a print device or for transferring them into a powerpoint or similar type of

Hi all, I have a rather naive question. I have the height of 100 individuals in a table and I want to assign the tallest 30% as Case=1 and the bottom 30% as Case=0. How do I do that? thanks. jiong The email message (and any attachments) is for the sole use of the intended recipient(s) and may contain confidential information. Any unauthorized review, use, disclosure or distribution is

Dear r-helpers, Could you please help me solve the following problem: When I run require(AlgDesign) trt <- LETTERS[1:5] blk <- 10 trtblk <- 3 BIB <- optBlock(~., withinData = trt, blocksizes = rep(trtblk, blk)) In response to the last command, R complains: Error in optBlock(~., withinData = trt, blocksizes = rep(trtblk, blk)) : object ".Random.seed" not found The

Dear r-helpers, When I do an xYplot and display the result in a JavaGD() window, the font is sans-serif (presumably Helvetica). When I send the figure to a pdf, I get a serif font (presumably times). How do I insure that the font in the pdf is indeed the default sans serif? > sessionInfo() R version 2.4.1 (2006-12-18) i386-apple-darwin8.8.1 locale: C attached base packages: [1]

In all likelihood, this doesn''t have anything to do with backgroundrb, but I''m stuck so I thought someone might be able to help. I have a pretty simple C extension that I built into a gem and have installed on my system. It works fine when required from rails (1.2.1). But when I have a Backgroundrb worker class use the extension (under os x), I get: