similar to: Matrix in R

Hi everybody. I'm trying to use R with Sweave but I have a problem perhaps with the directory path of sweave in R. The windows path is this: C:\Program Files (x86)\R\R-2.9.2\share\texmf\Sweave When I run the latex file with R, the program works well, without any errors, but when I create the pdf file, the Scode doesn't work. I think the problem is about the path. Can anybody suggest to

Hi R users. I've got a problem in producing the pdf file from Latex with R code. When I run the code Sweave("example.Rtex") in R it seems working, but when I run the Latex file it doesn't. The code error shown to me is below: *Runaway argument? {echo=FALSE} data<- read.csv("C:\\Users\\Daniele\\Desktop\\dati\\dati_england ! File ended while scanning use of \FV at

Hi R-users, I'm trying to built a plot of two series of data, but thees series result "superimposed". The R-code is like this: goal <- c(125, 143, 81, 26, 2, 3) numgoal <- 0:5 lambda <- sum(goal*numgoal)/sum(goal) plot(numgoal, goal, type="h") x <- 0:5 y<-dpois(x, lambda) att<-y*380 lines(x,att, type="h", col="red")

HI, I am using lmer() for a simple mixed effects model. The model is of the form logit(y)~ x + (1|z), where x is an indicator variable and z a multi-level factor. I would like an estimate of the response variable (either y or logit y) with an associated confidence interval for a given value of x. There does not appear to be a predict function written for lmer(). The output for the fixed

In trying to run John Fox's new Rcmdr package on Redhat Linux 9, I got this error message. The problem is in loading the tclck library. library(tcltk) produces the same error. It has been discussed before, e.g., http://finzi.psych.upenn.edu/R/Rhelp02/archive/11898.html but with respect to earlier versions of Linux and R. I tried updating to the latest RPMs of tcl and tk (8.3.5-89), but that

Hi everyone, I have a very quick question: Is there a ready-made function in R or any R packages to find the prime factorization of an integer? -- View this message in context: http://r.789695.n4.nabble.com/Prime-Factorization-tp2548877p2548877.html Sent from the R help mailing list archive at Nabble.com.

Hi, I have two vector with the marginal distribution like this: > a [1] -0.419 -0.364 -0.159 -0.046 -0.010 -0.002 0.000 0.000 0.000 > b [1] 0.125 0.260 0.270 0.187 0.097 0.041 0.014 0.004 0.001 How can I calculate the joint distribution with R? Thank you to all Dan -- View this message in context: http://r.789695.n4.nabble.com/Joiny-probability-with-R-tp3314059p3314059.html Sent

Dear all, I am looking at the R-help entry below: http://finzi.psych.upenn.edu/R/Rhelp02/archive/26640.html I have a more complicatedt problem. I have a zoo time series frame with 100+ sequences. I want to cycle through them back and forth and compare them to the 1st column at any time. I need also a button to click when I need the viewed-selected sequence (that is being compared to the

With the approach of R 2.12.0: with mild apologies for re-opening this perennial issue: is there any hope, if appropriate patches are submitted, of adding a drop.levels argument (with default equal to FALSE to preserve backward compatibility/efficiency) to the subset function ... ? If not, would a patch to the documentation and/or the R FAQ be accepted? This does seem to be a continuing

hi,R users, Now I have some scatter figures, is there some method can plot equalprobable error ellipse ? -- TANG Jie Email: totangjie@gmail.com Tel: 0086-2154896104 Shanghai Typhoon Institute,China [[alternative HTML version deleted]]

Hi, How can I generate data from multivariate gamma distribution & multivariate beta distribution? I only found command for multivariate normal only. Many thanks in advance :) Regards Rofizah [[alternative HTML version deleted]]

Does anyone know which package include the computation of Cochran’s Q statistic in R? jlfmssm [[alternative HTML version deleted]]

Hi, I am aware that this question might be a basic one. I did browse the help and archives but I still haven't understood how to do the following: I run R locally and would like to read in data from a server for which there is a username and password. That is, how do I open a connection to a server with a password? Thanks, Hadassa -- Hadassa Brunschwig PhD Student Department of Statistics

Dear All, Can you help me, with the code below how do I obtain the fisher information from it. Is my q<-replicate(1000,x) the right way to do simulation. thank you. x<-rweibull(100,0.8,1.5) q<-replicate(1000,x) z<-function(p){ beta<-p[1] eta<-p[2] log1<-(n*log(beta)-n*beta*log(eta)+(beta-1)*sum(log(x))-sum((x/eta)^beta)) return(-log1) } zz<-optim(c(0.5,0.5),z) zz Chris

Dear Group: I want to use lattice with a formula such as y ~ x | v to plot a data frame in which v varies to indicate which "x" is really being plotted. I know how to make the x-axis scales vary with the panel but is it possible to let the x-axis label vary, i.e., to let the user specify a vector of x-axis labels? Thanks Frank -- Frank E Harrell Jr Professor and Chair

Dear all! I am using wireframe function from lattice package. Is it possible to remove box around the plot but to keep axis (x, y, z)? Thanks in advance! [[alternative HTML version deleted]]

I am interested in studying the binary representation of numerics (doubles) in R, so am looking for possibilities of output of the internal binary representations. sprintf() with format "a" or "A" is halfway there: sprintf("%A",pi) # [1] "0X1.921FB54442D18P+1" but it is in hex. The following illustrate the sort of thing I want: 1.1001 0010 0001 1111

A common point made in discussion of contrasts, type I, II, III SS etc is that for sensible comparisons one should use contrasts that are 'orthogonal in the row-basis of the model matrix' (to quote from http://finzi.psych.upenn.edu/R/Rhelp02/archive/111550.html) Question: How would one check, in R, that this is so for a particular fitted linear model object? Steve Ellison

Hello UseRs, I need help on how to plot several raster images (such as those obtained from a kernel-smoothed intensity function) in a layout such as that obtained from the lattice package. I would like to obtain something such as obtained from using the "levelplot" or "xyplot" in lattice. I currently use: >par(mfrow=c(3,3) to set the workspace, but the resulting plots

Dear R-helpers, sorry if this is obvious, but I can't find this in the documentation. I am using Sweave, and have some code that does not actually work - but I want to include it anyway, including the error message that R produces. But on running Sweave() on my .Rnw file, it simply halts when it gets to the code chunk with the bad code. I can use \verbatim{} instead, but then I don't