similar to: inter-specific competition - community matrices and two species models using Lotka-Volterra

I have a simple (!) problem. What I want to do is to create a legend with Greek letters, and substituting numbers into the legend. Like this: mu=1:3 Mean=rep(mu, each=20) plot(runif(60), rnorm(60,Mean,0.1), pch=Mean) # create a plot legend(0.6,1.7, paste("mu =", 1:3), pch=mu) but with a Greek letter mu. I think the solution has something to do with substitute() possibly along with

Does anyone know of some code, and examples that implement game theory/Nash equilibrium hypothesis testing using existing packages like igraph/statnet or similar? Perhaps along the lines of this article: Zhang, Y., Aziz-Alaoui, M. A., Bertelle, C., & Guan, J. (2014). Local Nash Equilibrium in Social Networks, *4*, 6224. Best, Chris Buddenhagen cbuddenhagen at gmail.com [[alternative HTML

I am executing the following loops in R 15.1.0 and the first one works and produces results. However, the second one is not with initiating the object as x[i]. Further, I was able to get the second loop work in R 15.0. Am I missing something here? for(i in 1:length(b)) {y[i]<-paste(a,b[i],c,sep="") print(y[i]) print(fromJSON(file=y[i],method="C")) } for(i in 1:length(b))

(Using BRugs 0.7-5, R 2.14.2 32-bit on 64-bit Windows 7, OpenBUGS 3.2.1) 1. BRugs crashes R for me as follows. Sorry about the lack of detail; please let me know if / how to supply a more useful bug report on this issue. fit <- BRugsFit(...) # BRugs and OpenBUGS runs fine, the parameter estimates are reasonable # across 3 chains samplesBgr("beta") # crash

Dear users, I'm trying to get a dot plot but always end up with a boxplot. Can someone please tell me what I am doing wrong? df <- structure(list(FACETTE = structure(c(1L, 1L, 1L, 1L, 2L, 2L, + 2L, 2L), .Label = c("base", "tip"), class = "factor"), Sq = c(274836, + 0.74182, 0.709205, 0.984552, 279869, 255712, 26566, 301464)), .Names = c("FACETTE",

I try to do a test for dirichlet process for Multivariate normal, but Winbugs always says "expected multivariate node", does that mean I miss something at initialization? I will really appreciate the help to solve this problem Here is the R code, and Winbugs code. model { for(i in 1:N){ y[i,1:2] ~ dmnorm(mu[i,],tau[i,,]) S[i] ~ dcat(pi[]) mu[i,1:2] <- mu.star[S[i],]

I responded to the unhelpful suggestion "Why don't you implement and uplad the package to CRAN?" No mention of a search engine. Is this what you are commenting on Jeff? > On Jun 28, 2017, at 5:41 AM, Jeff Newmiller <jdnewmil at dcn.davis.ca.us> wrote: > > In what way does reminding people that packages exist because others just like them contributed something

I don't think OP asked an unreasonable question at all. Civility! > On Jun 27, 2017, at 2:00 PM, Suzen, Mehmet <mehmet.suzen at gmail.com> wrote: > > Why don't you implement and uplad the package to CRAN? > > On 27 Jun 2017 17:45, "Chris Buddenhagen" <cbuddenhagen at gmail.com> wrote: > > Does anyone know of some code, and examples that

Thanks I too wondered about the tone. The first suggestion was that I should "google it" and the second, write my own code. I think if I did I'd be reinventing the wheel, (and it'd be a big challenge for me). Also, I have been searching and not found such code, despite evidence that it has been coded (just not sure if it was in R). BTW I did write to authors of the article I

Hi. I am new here, and I am writing this Winbugs code with BRugs. n=length(bi.bmi) Lagegp=13 Lgen=2 Lrace=5 Lstra=15 Lpsu=2 #model gen x race bi.bmi.model=function(){ # likelihood for (i in 1:n){ bi.bmi[i]~ dbern(p[i]) logit(p[i])<- a0 + a1[agegp[i]]+a2[gen[i]]+a3[race[i]] + a12[agegp[i], gen[i]] + gam[stra[i]]+ u[psu[i],stra[i]] } # constraints for a1, a2, a3, a12 a1[1]<-0.0

Why don't you implement and uplad the package to CRAN? On 27 Jun 2017 17:45, "Chris Buddenhagen" <cbuddenhagen at gmail.com> wrote: Does anyone know of some code, and examples that implement game theory/Nash equilibrium hypothesis testing using existing packages like igraph/statnet or similar? Perhaps along the lines of this article: Zhang, Y., Aziz-Alaoui, M. A., Bertelle,

Hello Chris, I was implying you are capable enough to implement it, while you have already identify a research paper. If there is no package out there, uploading to CRAN would help future user too. I am more than happy to help if you want to implement from scratch. Best, Mehmet On 27 June 2017 at 17:45, Chris Buddenhagen <cbuddenhagen at gmail.com> wrote: > Does anyone know of some

In what way does reminding people that packages exist because others just like them contributed something count as being uncivil? Terse, perhaps, since it bypassed the obvious suggestion to use a search engine, but not rude. -- Sent from my phone. Please excuse my brevity. On June 28, 2017 5:08:16 AM EDT, Boris Steipe <boris.steipe at utoronto.ca> wrote: >I don't think OP asked an

Hi, I have a simple task I can't figure out. I'd like to take some measurements I made, e.g.: year (y-axis) 1 2 3 4 5 6 counts (x-axis) 10 10 20 30 40 50 And then, make a barplot with the x-axis ticks (representing the borders between years) between the bars. However, barplot seems to force you to make the x-axis arbitrary categories. I want it to be continuous (as in a

To account for likely differences between families in naturalization rates, we fitted a generalized linear mixed model, using PROC NLMIXED in SAS10, with the naturalization rate per genus (that is, the number of naturalized species in a genus as a proportion of the total number of introduced species in a genus) as the response variable, a variable coding genera as containing at least one native

Hi, I am trying to covert a Winbugs code into R code. Here is the winbugs code model{# model’s likelihoodfor (i in 1:n){time[i] ~ dnorm( mu[i], tau ) # stochastic componenent# link and linear predictormu[i] <- beta0 + beta1 * cases[i] + beta2 * distance[i]}# prior distributionstau ~ dgamma( 0.01, 0.01 )beta0 ~ dnorm( 0.0, 1.0E-4)beta1 ~ dnorm( 0.0, 1.0E-4)beta2 ~ dnorm( 0.0, 1.0E-4)#

Hello All, I am trying to compile Samba on POSIX under OpenVMS on a DEC 2000 AXP, however I am having difficulty defining the right FLAGSM and LIBSSM in the Makefile. So thus I'm having trouble getting the correct #includes, can anyone help? Regards, Gavin Butler -------------- next part -------------- A non-text attachment was scrubbed... Name: not available Type: application/ms-tnef

Hello, I am an R newbie, trying to use lsoda to solve standard Lotka-Volterra competition equations. My question is: how do I pass a parameter that varies with time, like say, phix <- 0.7 + runif(tmax) in the example below. # defining function lotvol <- function(t,n,p){ x <- n[1]; y <- n[2] rx <- p["rx"]; ry <- p["ry"] Kx <-

Hello, I have a syntactic problem with ode. How do I specify vectors of equations in ordinary differential equation systems. (i.e. in my case I want to simulate an a priory undefined number of species that have different parameters but the same behaviour) I demonstrate this using the Lotka Volterra example. The code below does not work and I have not a good idea how to specify this right. ##

hello, I want to iterate through a list of names and use each element as an argument in a function. For instance: > a = c('one','two','three') > data= c() > for(elem in a){data=cbind(elem = 2,data)} > data elem elem elem [1,] 2 2 2 instead I want 'elem' to be substituted by the string in the list. Doing it by hand would be: > data =