similar to: Empirica Copula

Dear all, I would like to plot under the same plot many curves like mycdf<-ecdf(runif(100)) plot(mycdf) curve(mycdf,from=-1,to=1) mycdf<-ecdf(runif(200)+3) plot(mycdf) curve(mycdf,from=-1,to=1) as you can see each new call to curve redraws (erases) the window. How I can have an effect like plot(mycdf) lines(mycdf2) lines(mycdf3) I would like to thank you in advance for your help

Hi r-users,   I really need help in fitting the t-copula.  I try to reproduce the example given by Jun Yan in “Enjoy the joy of copula” but I’m not sure how to correct the error based on the error message.  I tried so many ways but still could not get it working.   loglik.marg <- function(b, x) sum(dgamma(x, shape = b[1], scale = b[2], log = TRUE))   ctrl <- list(fnscale = -1)   #dat <-

Hi r-users, I really need help in fitting the t-copula. I try to reproduce the example given by Jun Yan in "Enjoy the joy of copula" but I'm not sure how to correct the error based on the error message. I tried so many ways but still could not get it working. loglik.marg <- function(b, x) sum(dgamma(x, shape = b[1], scale = b[2], log = TRUE)) ctrl <- list(fnscale

Hi r-users, I try to fit the t copula using the gamma marginals.  But I got error message which I don't really understand. Thank you for any help given. myCop.t <- ellipCopula(family = "t", dim = 2, dispstr = "toep", param = 0.5, df = 8) myCop.t myMvd <- mvdc(copula = myCop.t, margins = c("gamma", "gamma"), paramMargins = list(list(mean = 0, sd

Hello, I am trying to do a fit a loglikelihood function with Multivariate distribution via copulas with fitMdvc. The problem is that it doesn't recognize that my beta is a vector of km parameter and when I try to run it it say that the length of my initial values is not the same as the parameter. Can somebody guide me where my mistake is. Thanks, Elisa. #################################

Hello, Im trying to combine 3 affybatches (1x hgu133+2 array and 2x hgu133a array) Im useing this script: library(matchprobes) library(affy) library(AnnotationDbi) library(hgu133plus2probe) library(hgu133aprobe) library(hgu133a.db) u133p2 = ReadAffy() # reading hgu133 +2 cel file into affybatch u133a1 = ReadAffy() # reading hgu133a cel file into affybatch u133a2 = ReadAffy() # reading hgu133a

I am trying to make a loglikelihood function using copulas. I am trying to use mvdc to find the density function. When I run this I got the error that the pdf and cdf of my function tobit doesn't exist. Can somebody guide me where my mistake is? dtobit <- function(beta,sigma, x, y) {ifelse(y>0, dnorm(y,x%*%beta, sigma),(1-pnorm((x%*%beta)/sigma)))} ptobit <- function(beta,sigma, x,

Hi, I'd like to know why using the program "R" I can't add a number of margins> 3, I have a problem with the graphics. Post here my script: > myCop.norm <- ellipCopula(family = "normal", dim = 3, param = 0.4) > myMvd <- mvdc(copula = myCop.norm, margins = c("norm", "norm","norm"), > paramMargins = list(list(mean = 0, sd

ok, try it like this.   hosts allow = 127.0.0.1 192.168.1. EXCEPT 192.168.1.251   Greetz,   Louis   Van: Bernadette.McCluskey at gmb.org.uk [mailto:Bernadette.McCluskey at gmb.org.uk] Verzonden: donderdag 27 augustus 2015 14:36 Aan: L.P.H. van Belle Onderwerp: Re: [Samba] smb.conf I added a host ip to hosts deny in smb.conf and saved & restarted services. I then added that ip's

Hi! I want to do a meta-analysis with count data for treatement/control cases. Mi problem is that I need to use zero values (an informative value) for the mean and standard deviation for one of the treatement, but R has a problem: "Studies with zero values for sd.e or sd.c get no weight in meta-analysis". I can agroup the case by Family (byvar=Family). ¿anybody help me? Thanks! >

Hi, I want to do a meta-analysis with count data for treatement/control cases. Mi problem is that I need to use zero values (an informative zero value) for the mean and standard deviation for one of the treatement, but R has a problem: "Studies with zero values for sd.e or sd.c get no weight in meta-analysis". I can agroup the case by Family (byvar=Family). ¿Can you help me? Thanks!

Hola Comunidad, Tenia al siguiente duda ¿hay alguna forma de que R me de la Esperanza y Varianza de una distribución según sus parámetros? es decir, si tengo una weibull de parámetros "X" y "Y" para la esperanza hay una formula, pero me gustaría saber si hay una función en R que le diga algo como Esperanza(weibull, shape = "X", scale = "Y") y ya devuelva

I added a hosts deny line to above and denied one specific ip address. I then deleted that line but added hosts allow and made this equal to any ip address within previous "denied ip address's subnet. However, the afore-mentioned ip was still denied and I had to allow it individually. Can anyone tell me why. Thanks Bernadette McCluskey Ops/Comms Administrator GMB National

Hi all, I am going to sample two variables from two exponential distributions, but I want to specify a covariance structure between these two variables. Is there any way to do it in R? Or is there a "Multivariate Exponential" thing corresponding to the multivariate normal? Thanks in advance. Sincerely, Yanwei Zhang Department of Actuarial Research and Modeling Munich Re America Tel:

On Tue, Nov 23, 2010 at 10:15 AM, Bo Persson <bop at gmb.dk> wrote: >  > > > ----- Original Message ----- > From: Chandler Carruth > To: Zhanyong Wan (λx.x x) > Cc: Argyrios Kyrtzidis ; clang-dev Developers > Sent: Tuesday, November 23, 2010 6:40 AM > Subject: Re: [cfe-dev] draft rule for naming types/functions/variables > On Mon, Nov 22, 2010 at 9:36 PM,

Dear Group, if anyone has experience with the RMySQL Package maybe this behaviour is know: Reading data from a table into R the fields with datatype bigint(20) unsigned are transformed in some way: e.g. the query "select * from orders where userid = 14929859848712890325" selects the correct case but in R the userid is changed to 14929859848712890368. What happened here? This

Hi We want optimize my extensions file conf on asterisk 11.4.0 : We have a big quantity of extensions, all are same "design": ; Destination: Gambia Type: Fixe exten => _00220X.,1,Set(CDR(CodeCom)=BUS-GMB) exten => _00220X.,2,AGI(Caller-ID.agi,${IAXVAR(ACCOUNTID)}) exten => _00220X.,3,Set(CALLERID(all)=${NUMID}) exten =>

Greetings, I'm interested in generating data from various bivariate or mulitivariate distributions (e.g. gamma, t, etc), where I can specify the parameter values, including the correlations among the variables. I haven't been able to dig anything up on the faq, but I probably missed something. A nudge in the right direction would be appreciated. David --

hi, first I want to say that I'm new here, and new with copula and R. That is the reason why I'm writing, if somebody can help me. I have to make an example of Copula. On internet I've found this forum and that copula can calculate with R. Can somebody help me with the thing how can I start and where can read about these stuffs. Thank to all who can help! -- View this message

Hello, I am very new in R so I am so sorry for this question. I have the Barro-Lee data set which contains 98 countries and I want to run the regressions only for the Latin America countries, so what do you recomend? How can I delete all the other countries or how can I select the countries of Lat. Am. thank you this is the list of countries SHCODE COUNTRY NAME WBCTRY