search for: elmandjra

...pesetting R code. I encountered one problem: \begin{lstlisting} X %*% V \end{lstlisting} in the output the * in %*% disappears! same with %/%, etc, the / disappears. Any ideas? Kjetil -- Kjetil Halvorsen. Peace is the most effective weapon of mass construction. -- Mahdi Elmandjra

Suppose I have x<-data.frame(v1=1:4, v2=c(2,4,NA,7), v3=rep(1,4), v4=LETTERS[1:4],v5=rep('Z',4)) or a much larger frame, and I wish to test for and remove the constant numeric columns. I made: is.constant<-function(x){identical(min(x),max(x))} and apply(x,2,is.constant) # Works for numerics x[,-which(apply(x,2,is.constant))] I'd really like to be able to

...y.as.list.matrix(na.omit(cbind(Tmed, resid(VRS.mod1,type="response")))), xlab="", ylab=""))) there shoud be an easier way for such a daily task? Kjetil -- Kjetil Halvorsen. Peace is the most effective weapon of mass construction. -- Mahdi Elmandjra

...eck on a package with an R subdirectory empty, reports: * checking R files for syntax errors ... ERROR Syntax error in file On removing the empty subdirectory, the error disappears. Kjetil -- Kjetil Halvorsen. Peace is the most effective weapon of mass construction. -- Mahdi Elmandjra

...n halted na.omit of course is in package stats, and that is listed in the Depends field in the DESCRIPTION file. What is happening? (Windows XP with all development tools installed) Kjetil -- Kjetil Halvorsen. Peace is the most effective weapon of mass construction. -- Mahdi Elmandjra

Hello R'ers In previous versions of R (I now use 2.0.1) there was a package "bootstrap". I wrote some programs that depended heavily on this package but can unfortunately not find it on Cran today. Is there any way to find an older version of "bootstrap" and use it with the new version - 2.0.1? Sincerely Fredrik Lundgren Norrk??ping Sweden

...rom stats. Then VGAM:::family returns an object which doesn't have a $family component, (it has a component $vfamily). I thought namespaces should protect us from this happening? Kjetil -- Kjetil Halvorsen. Peace is the most effective weapon of mass construction. -- Mahdi Elmandjra -- No virus found in this outgoing message. Checked by AVG Anti-Virus.

...-1.077 ... $ y : num 0.453 1.227 -1.514 -0.904 -0.129 ... $ z : int 10 4 15 19 14 3 9 17 18 5 ... > summary(testframe) Error: protect(): protection stack overflow Kjetil -- Kjetil Halvorsen. Peace is the most effective weapon of mass construction. -- Mahdi Elmandjra -- No virus found in this outgoing message. Checked by AVG Anti-Virus.

...ame. data(vowel) loads the two first in the list above. Could it be this (which should be allowed, is mentioned in "writing R extensions") or is it something else, or a bug? Kjetil -- Kjetil Halvorsen. Peace is the most effective weapon of mass construction. -- Mahdi Elmandjra -- No virus found in this outgoing message. Checked by AVG Anti-Virus.

...ame. data(vowel) loads the two first in the list above. Could it be this (which should be allowed, is mentioned in "writing R extensions") or is it something else, or a bug? Kjetil -- Kjetil Halvorsen. Peace is the most effective weapon of mass construction. -- Mahdi Elmandjra -- No virus found in this outgoing message. Checked by AVG Anti-Virus.

Hi, Would you help me solve the following question? Thanks. Question: If I try to set the probability=0.05 and find the approximate x. (The answer should be somewhere between 2.1782 and 2.1783.) I write about this R program as follows but I don¡¦t know how to get the value of x which is between 2.1782 and 2.1783. library(mvtnorm) value<-array(1000) a<-array(1000)

Hello, I am hoping someone can help me with the following multivariate issue: I have a model consisting of about 50 covariates. I would like to reduce this to about 5 covariate for the reduced model by combining cofactors that are strongly correlated. Is there a package or function that would help me with this in R? I appreciate any suggestions. Thanks, Ian

Hello! I was reading R-intro and I have some suggestions: R-intro.html#A-sample-session rm(fm, fm1, lrf, x, dummy) suggestion rm(fm, fm1, lrf, x, y, w, dummy) The next section will look at data from the classical experiment of Michaelson and Morley to measure the speed of light. file.show("morley.tab") mm <- read.table("morley.tab") suggestion mm <- data(morley)

Dear List, We are planning a genotyping study to be analyzed using false discovery rates (FDRs) (See Storey and Tibshirani PNAS 2003; 100:9440-5). I am interested in learning if there is any consensus as to how many features (ie. how many P values) need to be studied before reasonably reliable FDRs can be derived. Does anyone know of a citation where this is discussed? Bill Dupont William D.

Dear R-users, I have an outcome variable and I'm unsure about how to treat it. Any advice? I have spending data for each county in the state of California (N=58). Each county has been allocated money to spend on any one of the following four categories: A, B, C, and D. Each county may spend the money in any way they see fit. This also means that the county need not spend all the money that

Dear R list members, Please excuse my ignorance but as a new comer to R I was wondering if anyone knows of any functions in R or Splus that can test a time-series for non-stationarity such as the Pettitt or the Mann-Kendall tests. Kind regards, Jon Nott [[alternative HTML version deleted]]

Dear R-Help readers, I suspect that this question must be a FAQ, but my investigation of the archives has not been very revealing. Is there an R function for calculating moving averages of time series objects? Thank you for your time and patience. -Paul Paul Schwarz, Ph.D. Associate Director of Methodology Gartner Vendor Marketing Solutions, Custom Research +1 503 241 8036 x186

What is the difference between arima and arima0?

Hallo, I would like to add the marginal distributions along the X and the Y axis to a scatter plot. Can anybody help me, please? Thank you, Paolo -- Paolo Bulla Istituto di Metodi Quantitativi Universit?? "L. Bocconi" viale Isonzo 25 20136 Milano paolo.bulla at unibocconi.it

Dear Sir/Madam: I'm currently running R 1.9.1 on a Mac OS 10.3. I have been able to use it fine for a while now. But for some reason, it now crashes each time I open it with the error: Fatal error: unable to restore saved data in .Rdata I read on the online help forum that I can just delete the bad .Rdata file. But I cannot find it. I even tried to delete everything related to R and