search for: huntsinger

Duncan Murdoch wrote: > Huntsinger, Reid wrote: > >>Sorry, I was looking at 2.0.1 when I meant to be looking at 2.1.0. The line >>numbers for the latter are 161-179 and line 164 is the one with what I think >>is a typo. >> >>Reid Huntsinger >> >>-----Original Message----- >>From: r...

I am currently the only use-R at my company, but they are considering buy a more powerful server and letting multiple people use it. They asked me if R supports client-server setups. I know S+ has a server version that does that. I didn't find anything about that on CRAN, but hopefully someone can correct me. I did see some stuff about R web servers

I would like to call C routines from Fortran under linux as suggested in section 5.6 of the "Writing R extensions" documentation. I'm familiar with Fortran but not with C. I understand the example provided in Fortran: subroutine testit() double precision normrnd, x call rndstart() x = normrnd() call dblepr("X was", 5, x, 1) call rndend() end but I don't understand

Hi, Does anyone know if R has the functionality to calculate a simple moving average. I cant seem to find it in the help menu. thanks, Wayne Dr Wayne R. Jones Statistician / Research Analyst KSS Group plc St James''s Buildings 79 Oxford Street Manchester M1 6SS Tel: +44(0)161 609 4084 Mob: +44(0)7810 523 713 KSS Ltd A division of Knowledge Support Systems Group plc Seventh

Thanks very much, Professor Ripley. Reid Huntsinger -----Original Message----- From: ripley at stats.ox.ac.uk [mailto:ripley at stats.ox.ac.uk] Sent: Tuesday, August 06, 2002 3:00 AM To: Huntsinger, Reid Cc: r-help at stat.math.ethz.ch Subject: Re: [R] pipe and binary i/o (on Linux) On Mon, 5 Aug 2002 ripley at stats.ox.ac.uk wrote: > pip...

Hi, It is my understanding that the eigenvectors of a circulant matrix are given as follows: 1,omega,omega^2,....,omega^{p-1} where the matrix has dimension given by p x p and omega is one of p complex roots of unity. (See Bellman for an excellent discussion on this). The matrix created by the attached row and obtained using the following commands indicates no imaginary parts for the

Assume we have a function like: foo <- function(x, y) how is it possible to define a binary indexing operator, denoted by $, so that x$y functions the same as foo(x, y)

Hi guys im having a problem getting R to numerically integrate for some function, say f(bhat)*optimize(G(bhat)), over bhat. Where id like to integrate this over some finite range, so that here as we integrate over bhat optimize would return a different optimum. For instance consider this simple example for which I cannot get R to return the desired result: f <- function(bhat) exp(bhat) g

Dear helpers in this forum, This is a clarified version of my previous questions in this forum. I really need your generous help on this issue. > Suppose I have the following data set: > > id x y > 023 1 2 > 023 2 5 > 023 4 6 > 023 5 7 > 412 2 5 > 412 3 4 > 412 4 6 > 412 7 9 > 220 5 7 > 220 4 8 > 220 9 8 > ...... > Now I want to compute the

...om the node they've migrated to you also need DFSA enabled. That sounds like a great idea but it seems to cause problems in our setup--that is, with DFSA off we go a lot longer without oMFS misbehavior or other odd phenomena than when DFSA is on. How many processes do you run at a time? Reid Huntsinger -----Original Message----- From: Jim Thomas [mailto:james at staarfunds.com] Sent: Wednesday, April 21, 2004 3:14 PM To: Huntsinger, Reid Cc: r-help at stat.math.ethz.ch; openmosix-general at lists.sourceforge.net Subject: Re: [openMosix-general] openMosix and R: File I/O issues? Memory could b...

...y contains your problem (section 3, Testing the Rasch model) and gives a very simple Markov chain for sampling from the uniform distribution on these matrices. If you need other than the uniform distribution, see the modifications Diaconis and Sturmfels make (the "Metropolis" step). Reid Huntsinger -----Original Message----- From: r-help-bounces at stat.math.ethz.ch [mailto:r-help-bounces at stat.math.ethz.ch] On Behalf Of Huntsinger, Reid Sent: Thursday, January 27, 2005 10:50 AM To: 'Michela Marignani'; r-help at stat.math.ethz.ch Subject: RE: [R] Request for help Persi Diaconis...

Hi, I am trying to run a very computationally expensive procedure in R-2.0.0. and the process always gets killed after approx 8 minutes. This procedure calls some of my own C++ code - in case it was this code causing a memory leak I unload and then reload the .so file every time, however I still get the same problem. The procedure is run 16000 times and always calls the lm() function. My

Could anyone help me with the following issue. Using the GSL library in R I define the following code: ######### library(gsl); S<-function(n) { r<-0:n; ans<-sum(gsl_sf_choose(n,r)*(-1)^r*2^(2*r)*gamma_inc(6-2*r,2)) ans } ######### >SS(10) yields 34.91868 >SS(40) yields 5.340422 >SS(60) yields 180.3162 Doing the same computations in maple I get

Good afternoon, do you know if R provides a special package to make networks ? Thank you, Angela

Hello! Does anyone now if there exist a function that would compute a "null space" (or "kernel" - "Ker") of a matrix and maybe also one that would compute an "image" ("Im") of a matrix. I tried R-site search and google, However I found notnihg useful! Thanks for any sugestions! I am also not sure what an "image" of a matrix is, so

...of the mode. This seems to work (fingers crossed). Then I noticed gzfile has the same problem, but bzfile doesn't, and uses the same solution as I did for pipe. Can I assume that this is a correct fix? Or were there reasons that pipe and gzfile wouldn't allow "rb"? Thanks, Reid Huntsinger -----Original Message----- From: Thomas Lumley [mailto:tlumley at u.washington.edu] Sent: Monday, August 05, 2002 12:17 PM To: Jim Lemon Cc: r-help at stat.math.ethz.ch Subject: Re: [R] filename of source file On Mon, 5 Aug 2002, Jim Lemon wrote: > David Richmond wrote: > > Yes, sourc...

...ss[i],1,word] + 1 } and the "absent" (==2) parts are then easy: tab[1,2,] <- sum(class == 1) - tab[1,1,] tab[2,2,] <- sum(class == 2) - tab[2,1,] so now you can use chisq.test on each of the 2 x 2 tables tab[,,i] for i a word index, all at once using apply() if convenient. Reid Huntsinger -----Original Message----- From: r-help-bounces at stat.math.ethz.ch [mailto:r-help-bounces at stat.math.ethz.ch] On Behalf Of Weiwei Shi Sent: Wednesday, June 15, 2005 5:10 PM To: R-help at stat.math.ethz.ch Subject: [R] coding to generate a matrix to prepare for chi-sqr test for text mining Hi...

...Win32 system x86, Win32 status major 1 minor 4.1 year 2002 month 01 day 30 language R > Following is a diagnosis of the problem by one of my colleagues: > From: Huntsinger, Reid > Sent: Monday, April 22, 2002 5:15 PM > To: Raubertas, Richard > Subject: RE: R problem > > [ ... ] This happens > because of the attempt to convert. The function type.convert > is called on columns with unspecified types; this tries to > convert to "logica...

...icant compared to 1, say, then you get essentially the "theoretical quantiles" ie the "x" component of the list back from l$x - l$y. l$x is basically a sample from a normal(0,1) distribution so they do line up perfectly in the second qqnorm(). Is that what's happening? Reid Huntsinger -----Original Message----- From: r-help-bounces at stat.math.ethz.ch [mailto:r-help-bounces at stat.math.ethz.ch] On Behalf Of WeiWei Shi Sent: Thursday, April 28, 2005 1:38 PM To: Vincent ZOONEKYND Cc: R-help at stat.math.ethz.ch Subject: [R] have to point it out again: a distribution question...

Hi - I am doing a monte carlo experiment that requires to do a linear regression of a matrix of vectors of dependent variables on a fixed set of covariates (one regression per vector). I am wondering if anyone has any idea of how to speed up the computations in R. The code follows: #regression function #Linear regression code qreg <- function(y,x) { X=cbind(1,x) m<-lm.fit(y=y,x=X)