similar to: Remove

Hi, I'm writing a short course tutorial to Browninan Motion/ Langevin Equation. At the end of the theory section I wanted to add a short GNU R example, so the students can play a little around. I already looked in the MASS book (by Venables and Ripley) but I couldn't find any Brownian Motion/ Langevin Equation package. Are there any good packages or tutorials available which cover R and

Dear all, I have been looking for a while for ways to simulate from Langevin distributions and I thought I would ask here. I am ok with finding an algorithmic reference, though of course, a R package would be stupendous! Btw, just to clarify, the Langevin distribution with (mu, K), where mu is a vector and K>0 the concentration parameter is defined to be: f(x) = exp(K*mu'x) / const where

Hello, We use the mix package and we have a problem with the DA function. We aren't sure, but it's maybbe a memory problem. We have done: > Ent<--read.table("C:/.../File.txt") > attach(Ent) > Ent V1 V2 V3 V4 ... V16 V17 1 1 1 2 6 18 18 2 1 1 1 NA 14 17 3 1 1 2 1 16 14 .... 199 2 1 NA 7 19 18 200

I'm using R version 1.1.0 with Redhat 6.2 in a Dual-Pentium When I type postscript() I have the next error: -- Error in old$command == "default" : comparison (1) is possible only for vector types I have updated R next week using the rpm package http://www.stat.math.ethz.ch/CRAN/bin/linux/redhat/6.x/i386/R-base-1.0 .0-2.i386.rpm I know that postscript() has new features, but I

Hi, I''m having trouble connecting to a remote database (mysql). I know that my connection settings are correct, since I have a php app that is connecting to the same database. Are there any settings besides the following that I need in my database config file? Any other suggestions? development: adapter: mysql database: db_name username: username password: password host:

Hello, I've just installed R (Linux) and I have some problems with graphical viewing.Could you, please, explain me why I can't see graph when I just want to make a plot ? Many thanks ! St?phanie Langevin Helios BioSciences ---------------------------------------------------------------- Ce service de mailing vous est offert par http://www.freesurf.fr. FreeSurf, votre acces ADSL a

Hi! I am using R1.5.0. When using functions form lattice the background is dark grey. I have tried to change it using par() but it does not work. By the way, if I use par and the standard functions like hist() there is no problem with the par function. Any hints? Thank in advance, Ignacio Perez Escuela Colombiana de Ingenieria

Dear useRs, I'm having troubles with using ylim in barplot(): even though I reduce the y-scale using ylim, the bars still extend down to 0into the x-labels. The sample data is below, and here is the code. #This works fine but I would like to plot only from 50 to 70: barplot(t(mydata), beside=TRUE, col=c("orange", "green", "yellow", "purple"),

Hi, I've used various patches with asterisk 1.4 to have support for call pickup and notification with good results. Now I'm trying vanilla 1.6.2 with its official support for "dialog-info +xml" notifications with no success. This is what i'm doing: - Phone A has a key configured as type "extension" pointing to Phone B. - In sip.conf I added

Hi, I have two questions. First, I'd like to compare a ZINB model to a negativ binomial model with the Vuong test, but I can't find how to performe it from the zicount package. Does a programm exist to do it ? Second, I'd like to know in which cases we have to use a double hurdle model instead of a zero inflated model. Many thanks, St??phanie Payet REES France R??seau

Thank you David, That's a nice workaround using plotrix::barp(), but that doesn't explain why ylim doesn't work as intended (or at least, as I expect it to work), or why xpd has no influence when using devEMF::emf()... The problem with saving directly in RStudio is that it requires to manually save the plot, and this becomes troublesome when there are a lot of plot commands in a

Hi - I would like to create a scatterplot with the symbols coded by a factor in each record. Is their a better way than what I have used below? I attempted to use the plot.factor, but cannot force it to make a scatterplot instead of boxplots. I have considerd making a function with a 'for' statement for each level of the factor. Is there a better way for when I am dealing with more

Hola! Necesito evaluar la normalidad de una variable para la que tengo más de 5000 observaciones. Shapiro-Wilks no funciona para muestras mayores 5000 observaciones. AAlshap <- lapply(AAdf, shapiro.test) Error in FUN(X[[i]], ...) : sample size must be between 3 and 5000 Alguna sugerencia? Gracias. [[alternative HTML version deleted]]

Hi everyone I'm trying to use sigma2.mat<-apply(ri.mat,2,'-n/2*log(sum(ri.mat^2)/n)') but it's not working. Is it wrong to use that kind of FUN in apply ? or am I doing something wrong ?! Thank's ernesto -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send

If you are using a Windows system, you can Export the plot from RStudio and save it as a metafile without using package devEMF and it will crop the bars with xpd=FALSE. When I used devEMF on a Windows machine, the bars were not cropped with barplot() as you indicated, but when I switched to plotrix::barp() they were cropped. The arguments are a bit different, but I did not need xpd=FALSE:

Dear R-users, Is there a function in R that classifies data in more than 2 groups using logistic regression/classification? I want to compare the c-indices of earlier research (lrm, binary response variables) with new c-indices obtained from 'multiple' (more response variables) logistic regression. Best regards, Stephanie Delalieux Department Biosystems M?-BIORES Group of Geomatics

Hola a todos. Soy un poco torpe manejando cadenas de texto, así que os pido ayuda. Tengo un vector de texto de este tipo datos$tipo [1] m.1.p.Álava m.1.p.Albacete [3] m.2.p.Alicante m.1.p.Almería [5] m.3.p.Asturias m.1.p.Ávila [7] m.1.p.Badajoz m.1.p.Baleares (Illes) [9] m.1.p.Barcelona m.1.p.Burgos [11] m.1.p.Cáceres m.1.p.Cádiz Y quiero extraer el

Hi guys, I resently used xapian to sort some documents by distance between 2 points. I implemented a MatchDecider which work well. I now tried to implement a Weight scheme to put my document in ascending order depending on the distance... My information to calcul distance is in values in the document. How I can access document values from Weight to be able to add some sum_extra weight ??

Hello, I have a general formatting question. I have two columns of data: ColA <- c("m", "m", "m", "m") ColB<- c("d","d","d","d") And I would like to reorder them into a new column that looks like this: ColC<-

When I use the command for PACF, lags 5, 10, 15, and 20 are labeled. I would like to label lag 1. I would greatly appreciate if someone could tell me how to do this. Below is the command that I am using: pacf(data$R1,main="Series R1 Residuals") [[alternative HTML version deleted]]