similar to: Avoiding loops to spare time and memory

The update of VR by downloading VR_7.1-6.zip and using install.packages (from local zip files) fails with the following error message: Error in file(file, "r") : unable to open connection In addition: Warning message: cannot open file `VR/DESCRIPTION' Other packages can be installed without problems, except of dse_2003.4-1 with a similar error message. Why? Operating System:

Can someone show me how to specify zero inflated poisson and zero inflated negative poisson models in R? I would like to replicate the example given in Long (1997: Regression Models for Categorical and Limited Dependent Variables) in Chapter 8.5 (pp. 242-247). TIA Dirk ************************************************* Dr. Dirk Enzmann Criminological Research Institute of Lower Saxony

I wonder whether the functions for skewness and kurtosis in the e1071 package are based on correct formulas. The functions in the package e1071 are: # -------------------------------------------- skewness <- function (x, na.rm = FALSE) { if (na.rm) x <- x[!is.na(x)] sum((x - mean(x))^3)/(length(x) * sd(x)^3) } # -------------------------------------------- and #

Just out of curiosity (and I hope that my question does not confuse those who intuitively understood the par(new=TRUE) statement right in the beginning): I would like to know whether beginners (learning R) had the same difficulty as I had: Intuitively I thought that par(new=TRUE) would draw a new plot (and not into an already existing plot) and that par(new=FALSE) would not draw a new plot

Trying to select a subset of cases (rows of data) I encountered several problems: Firstly, because I did not read the help to read.spss() thoroughly enough, I treated the data read as a data frame. For example, dr2000 <- read.spss('myfile.sav') d <- subset(dr2000,RBINZ99 > 0) and thus received an error message (Object "RBINZ99" not found), because dr2000 is not a

I receive R-help messages in digest form that makes it difficult to answer a post. I noticed that my answer is not added to the thread (instead, a new thread is started) although I use the same subject line (starting with "Re: ") as the original post. Is there a solution (I prefer the digest to separate messages for several reasons and don't want to change my email reader)? The

In trying to streamline various optimization functions, I would like to have a scratch pad of working data that is shared across a number of functions. These can be called from different levels within some wrapper functions for maximum likelihood and other such computations. I'm sure there are other applications that could benefit from this. Below are two approaches. One uses the <<-

Using ranef() (package nlme, version 3.1-75) with an 'lme' object I can obtain random effects for intercept and slope of a certain level (say: 1) - this corresponds to (say level 1) "residuals" in MLWin. Maybe I'm mistaken here, but the results are identical. However, if I try to get the standardized random effects adding the paramter "standard=T" to the

After loading the package "Matrix" (version 0.995-12), using the summary function with an lme (package nlme version 3.1-75) object results in an error message saying Fehler in dim(x) : kein Slot des Namens "Dim" f?r dieses Objekt der Klasse "correlation" (translated: 'Error in dim(x) : no slot of the name "Dim" for this object of class

Although I tried to find an answer in the manuals and archives, I cannot solve this (please excuse that my English and/or R programming skills are not good enough to state my problem more clearly): I want to write a function with an indeterminate (not pre-defined) number of arguments and think that I should use the "..." construct and the match.call() function. The goal is to write

I hesitate to post this question in the light of recent threads, indeed I have hesitated for several weeks, however I have come to a full stop and really need some help if I am going to progress. I am a new user of R for medical statistics. I have attempted to read all the relevant documents, but would welcome any suggestions as to what I have missed. I am trying to contruct "table 1"

Hello, I just set up a samba server (version 2:4.1.13+dfsg-4) on Debian 7.7, being a member server of a Win2k8-Domain (before that, that server was an old SuSE (10.4)-Samba with own user-management (standalone-server). I would like to use winbind with the idmap backend "ad" or "rfc2307" instead. When using rfc2307 (like in my conf specified), I can do successfully: wbinfo

Dear All, I have a codes which calculates the result of Ripley's K function of my data. I want to repeat this process 999 times. However, i am getting an error when i use the "for i in" function. Is there any way to repeat this analysis 999 times. Here are the codes i used ; data4 <- matrix(c(sample(id),data1),203,3) a <- data4[,1] random.case=data4[a==0,]

Hi! Since we upgraded to 1.6.1.11, asterisk is only outputing monitored files -in and -out. It is not mixing them in the end. queues.conf has monitor-type=MixMonitor... Would somebody help me debug why it doesn't mix the sounds?? Thanks -------------- next part -------------- An HTML attachment was scrubbed... URL:

Dear friends, I am studying the spatial distribution of trees in a tropical swamp forest in Brazil, and the spatial association between the trees and other forest and non forest elements. These are point pattern analysis of mapped data, preferentially done with the L function of Ripley. I have recently discovererd the R-project and would like to know whether it is possible to do these

One Question: if I make a lattice-plot (beware: under development!) with a special layout e.g. c(2,2,4): p1 <- xyplot(y~x|id,data=data,layout=c(2,2,3)) is there a way to print.trellis e.g. the second page of the trellis object "p1"? print.trellis(p1,page=2) #!wrong! I want this to automatically dev.print the pages (no, I don't want multiple postscript pages). If this

One of our students did something like: x[a < b < c] instead of x[a < b & b < c] But why is 3 < 2 < 1 # [1] TRUE ??? Is there any reason? Or wouldn't it be better to get a warning / error? Uwe Ligges -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html

Hi LLVM community, I'd like to ask a question regarding the behavior of llvm-link: My code contains Classes which are structurally equivalent but they are totally unrelated and distinct on a c++ point of view. However, if the compiled IR gets processed by llvm-link, these types are merged together. My question is: Is this expected behavior or a bug? To explain it more in detail, a reduced

I have to admit that I'm a bit ignorant about regular expressions...I have a problem with "gsub": under unix I can : > cat tmp.txt F.123 F.123 F123 F123 sed 's/F\./d/g' tmp.txt d123 d123 F123 F123 that is: replacing "F.", however under R: gsub("F\.","d",c("F123","F.123"),extended=F) [1] "d23"

Found this. http://www.tnt.uni-hannover.de/~baumgart/peri_ear_v0.1.tgz The peripheral ear model is based on the structure of Eberhard Zwicker's "Analogmodell". That model consists of analog electrical elements. ftp://ftp.tnt.uni-hannover.de/pub/papers/2000/Diss-FB.pdf Input: Raw 100KHz monoural 16bit audio Output: hair cell excitation signal at 251 equidistant locations in