similar to: mean reverting model

Is there a package that will allow me to fit Brownian motion and Ornstein?Uhlenbeck models of evolution for discrete traits? I know that geiger and ouch have commands for fitting these models for continuous traits, but these aren't suitable for discrete trait evolution, correct? -- View this message in context:

Dear friends i am trying to fit an Ornstein-Uhlenbeck process by MAXIMUM LIKELYHOOD method. i found these formulas on http://www.sitmo.com/article/calibrating-the-ornstein-uhlenbeck-model/ this is the mean-reverting process http://r.789695.n4.nabble.com/file/n4637271/process.txt process.txt and this is the script that i am using....... ouFit.ML=function(spread) { n=length(spread)

hello everyone, i'm trying to simulate a 2-parameter integrated ornstein-uhlenbeck (IOU) process, but i'm not sure exactly where to start (which package, which function). the motivation is the paper by taylor et. al. (JASA 1994) "a stochastic model for the analysis of longitudinal aids data." the model they suggest consists of a combination of fixed and random effects, a

dear r-community, currently i have a workflow as following: data in online-MYSQL-database via .csv into R for analysis, results back to mysql. generating graphics for reports also. that's ok, but i need some process-optimization. i'd like to run r directly on the webserver. analysis should work automatically, after the appropriate r-script was selected in a browser-gui for R. until

Dear all, I wanted to let you know about four phylogenetic comparative methods (PCM) packages that have become available on (3 on CRAN and 1 on GitHub) recently that hopefully will be interesting to somebody. Three of them go significantly beyond the Brownian motion (BM) and Ornstein-Uhlenbeck (OU) processes. 1) There is a new version of mvSLOUCH available. The most important change is that the

Dear all, I wanted to let you know about four phylogenetic comparative methods (PCM) packages that have become available on (3 on CRAN and 1 on GitHub) recently that hopefully will be interesting to somebody. Three of them go significantly beyond the Brownian motion (BM) and Ornstein-Uhlenbeck (OU) processes. 1) There is a new version of mvSLOUCH available. The most important change is that the

Hello. I try to compute a "White-corrected" ANOVA. Therefore I use the command "Anova" implemented in the car()-package): Anova(modelname, white.adjust="hc3") and receive the error message: "Fehler in SS[i] <- SS.term(names[i]) : nichts zu ersetzen" (German) "Error in SS[i] <- SS.term(names[i]) : nothing to replace" (probably being the

Hi R Users, I have produced a simulated scatter plot of y versus x tightly clustered around the 45 degree line through the origin with the following code: x <- seq(1,100) y <- x+rnorm(100,0,10) plot(x,y,col="blue") abline(0,1) Is there some way to generate a 3-dimensional analogue of this? Can I get a similar simulated scatter plot of points in 3 dimensions where the points

Hi, I have a matrix of results of multiple 2x2 chi^2 tests, non- significant tests are marked as TRUE. Is there a function for grouping those variables in a similar way LSD.test from agricolae library does? I reviewed LSD.test's source but it's not helpful for me. This is my matrix: [,1] [,2] [,3] [,4] [,5] [,6] [,7] 1 TRUE FALSE TRUE TRUE TRUE FALSE TRUE 2 FALSE TRUE

Dear all, I have just uploaded to CRAN the version 1.6 of the package 'adehabitat'. Significant changes are listed below: * The package has been reorganized into four parts (see ?adehabitat-package for a description): (i) management of raster maps, (ii) habitat selection / ecological niche analysis, (iii) home range analysis, and (iv) analysis of animals trajects. The package contains

Dear all, I have just uploaded to CRAN the version 1.6 of the package 'adehabitat'. Significant changes are listed below: * The package has been reorganized into four parts (see ?adehabitat-package for a description): (i) management of raster maps, (ii) habitat selection / ecological niche analysis, (iii) home range analysis, and (iv) analysis of animals trajects. The package contains

Hi R users, I am getting the following error when using the splsda function in R v2.12.1: "Error in switch(classifier, logistic = { : EXPR must be a length 1 vector" What does this mean and how do I fix this? Thank you in advance! Best, Savi

Dear R-users, I have installed the latest version of Tinn R (Version 2.2.0.2). With the older versions it was possible to run Tinn R next to R and then select code that one could send to R. Now it seems that one has to start the R gui in Tinn R and then it is possible to select en run code of R. But if I start R within Tinn R it does NOT start the R-Gui but only the R console. When I start R

Dear list, I have a question and I'm going to give an example of my problem f<- function(d1, d2, d3) { d<- d1*d2/d3 return(d) } v1<- 1 v2<- 2 If I try f(v1, v2, v3) Error in f(v1, v2, v3) : object "v3" not found I obviously got the above error message. I would like to add something to my function to allow me to get a certain value (say zero) if one of the arguments is

Hi, when I have made a decision tree with rpart, is it possible to "apply" this tree to a new set of data in order to find out the distribution of observations? Ideally I would like to plot my original tree, with the counts (at each node) of the new data. Reagards, Jay

In a previous message I was not clear enough in my querry. I have the following program: tst<- function() { x <- c(32.7,32.3,31.5,32.1,29.7,29.1,35.7,35.9,33.1, 36.0,34.2,31.2,31.8,28.0,29.2,38.2,37.8,31.9, 32.5,31.1,29.7) g <- rep(1:7,rep(3,7)) s <- rep(1:3,7) cat(" Only x and g \n") aov1(x,g) cat("\n\n Now x, g and s \n") aov1(x,g,s=s) }

Hello ,everyone. I am a fresh user of R. I wrote several several R script in R with the version 2.9.1 in Windows XP. But I only can run these scrpit via R-GUI. Can anyone tell me how to run the scripts in background without openning R-gui.In another word, run the R-script and get the result via some batch proceesing files or shell scripts directly? Thanks. . -- TANG Jie Email:

Hi, I am looking for a R package to draw 3d plot. But not in a static way like scatterplot3d or stuff like that. I would like to make rotate the plot, to zoom/unzoom etc..., actually to render the graph 'dynamic'... (for the biologist and bioinformatician, a R package which permitts to reproduce what pdb viewer does). Does a Way exist to do something like that in R ? Thanks,

Hi R helpers, If I have a dataset looks like: ID record 1 20 . 30 . 25 2 26 . 15 3 21 4..................... And I want it becomes ID record 1 20 1 30 1 25 2 26 2 15 3 21 4..................... That is, I have to duplicate IDs for those with multiple records. I am wondering it is possible to be

An indicator matrix is a binary matrix with orthogonal columns whose rows sum to 1. A row of this matrix could be [0 1 0 0]. My problem is to group the similar rows (profiles) so that to create a compact form of the matrix. Is there an R function that deals with this problem or do I have to write it from scratch? Thanks, Angelos Markos Dr. Applied Informatics, University of Macedonia, Greece