search for: simul

...x2,x3,x4),ncol=5) Y <- matrix(y,ncol=1) Design.data <- cbind(X, Y) for (j in 1:nrow(X)) { result <- vector("list", ) for( i in 1: 3100) { data <- Design.data[sample(50,50,replace=TRUE),] ##### and unique(Design.data.....) dataX <- data[,1:5] dataY <- data[,6] B.cap.simulation <- solve(crossprod(dataX)) %*% crossprod(dataX, dataY) P.simulation <- dataX %*% solve(crossprod(dataX)) %*% t(dataX) Y.cap.simulation <- P.simulation %*% dataY e.simulation <- dataY - Y.cap.simulation dX.simulation <- nrow(dataX) - ncol(dataX) var.cap.simulation <- crossprod...

Um texto embutido e sem conjunto de caracteres especificado associado... Nome: n?o dispon?vel Url: https://stat.ethz.ch/pipermail/r-help/attachments/20070508/1c7def0a/attachment.pl

Hello, All: ????? The default "simulate" method for lm and glm seems to ignore the sampling variance of the parameter estimates;? see the trivial lm and glm examples below.? Both these examples estimate a mean with formula = x~1.? In both cases, the variance of the estimated mean is 1. ??? ??????? * In the lm example with x...

I have found the "simulate" method (incorporated in some packages) very handy. As far as I can tell the only class for which simulate is actually implemented in base R is lm ... this is actually a little dangerous for a naive user who might be tempted to try simulate(X) where X is a glm fit instead, because it defaul...

On 2019-12-27 04:34, Duncan Murdoch wrote: > On 26/12/2019 11:14 p.m., Spencer Graves wrote: >> Hello, All: >> >> >> ? ????? The default "simulate" method for lm and glm seems to ignore the >> sampling variance of the parameter estimates;? see the trivial lm and >> glm examples below.? Both these examples estimate a mean with formula = >> x~1.? In both cases, the variance of the estimated mean is 1. > > That...

I have been talking with Gabriel Wainer about the possibility of some community building with R and simulation workers. The conference announced below is just before UseR (but in that capital to the north of Washington). If anyone on this list is participating or thinking of participating, perhaps they could let me know off-list. Ottawa-Gatineau R User Group would likely try to set up some informal...

Dear R helpers Suppose I have a portfolio of securities with exposure to Equity, Bonds and Forex (say $ 1000000 each). Is there any fucntion in R that will help me calculate Value at Risk (VaR) using Monte carlo Simulation , Historic simulation and Variance - Covariance Simulation. With regards Maithili

...ll, I have the following R script : #!/usr/bin/Rscript out_file = "hybrid.pdf" pdf(out_file, height=8.5, width=11) myvalues_1M <- read.csv("hybrid_sims_1M.csv",head=TRUE,sep=",") plot(myvalues_1M$num_sims_per_thread,myvalues_1M$time_per_sim,xlab="Number of Simulations per Thread",ylab="Time per 1 million Simulations(in milliseconds)",col="red",main="For 1 million simulations") lines(myvalues_1M$num_sims_per_thread,myvalues_1M$time_per_sim,col="red") myvalues_5M <- read.csv("hybrid_sims_5M.csv",head...

Full_Name: Reginaldo Constantino Version: 2.8.0 OS: Ubuntu Hardy (32 bit, kernel 2.6.24) Submission from: (NULL) (189.61.88.2) For many tables, chisq.test with simulate.p.value=TRUE gives a p value that is obviously incorrect and inversely proportional to the number of replicates: > data(HairEyeColor) > x <- margin.table(HairEyeColor, c(1, 2)) > chisq.test(x,simulate.p.value=TRUE,B=2000) Pearson's Chi-squared test with simulated p-value...

On 26/12/2019 11:14 p.m., Spencer Graves wrote: > Hello, All: > > > ????? The default "simulate" method for lm and glm seems to ignore the > sampling variance of the parameter estimates;? see the trivial lm and > glm examples below.? Both these examples estimate a mean with formula = > x~1.? In both cases, the variance of the estimated mean is 1. That's how it's docu...

Dear R users, I need to run 1000 simulations to find maximum likelihood estimates. I print my output as a vector. However, it is taking too long. I am running 50 simulations at a time and it is taking me 30 minutes. Once I tried to run 200 simulations at once, after 2 hours I stopped it and saw that only about 40 of them are simulated...

Hi All, I,m trying to do a GARCH simulation in R 2.3.0 release in Windows XP. I've seen garchsim function but that is for garch (1,1) and ?garch gives an example for ARCH simulation. Can anyone help me how can i extend the help shown in ?garch to GARCH simulation? Please help me in this regard. Thanks, Sumanta Basak.

...of experimental data (with replicae in time) >experiment<-data.frame(times=c(0,0,10,10,20,20,30,30),expval=c(1,1,2,2,3,3,4,4)) > experiment times expval 1 0 1 2 0 1 3 10 2 4 10 2 5 20 3 6 20 3 7 30 4 8 30 4 and a set of simulated values > simul<-data.frame(times=c(0,10,20,30),simul=c(3,4,5,6)) > simul times simul 1 0 3 2 10 4 3 20 5 4 30 6 >From this simulated data, I would like to obtain a data frame like this times expval simul 1 0 1 3 2 0 1 3 3...

Hello, can anyone tell me if R has any special function for simulating the structure of human populations? Something like the genetic algorithm? I need to simulate a sample of a population with a specific structure. Is there something on R that can help me? Thanks to everyone. __________________________________________________________

I'd like to carry out a Monte Carlo simulation test where given data is a contingency table. I think this is something to do with using rmultinonom(), but I'm not sure how to code this, to simulate contingency tables. Could anyone please help with how to use R to simulate contingency tables like this? -- View this message in context:...

Dear All, I have been trying to simulate data from a fitted glm using the simulate() function (version details at the bottom). This works for lm() fits and even for lmer() fits (in lme4). However, for glm() fits its output does not make sense to me -- am I missing something or is this a bug? Consider the following count data, model...

Hello, I created an array to hold the results of a series of simulations I'm running: d.eta <- array(0,dim=c(3,3,200)) <simulation goes here and populates the array but it's not important> Then I tried to save the results using this: save(d.eta,file="D:/Simulation Results/sim 9-23-11 deta") When I later tried to reload them using th...

Hi, sorry for bothering your guys again. I want to simulate 100 AR(1) data with cor(x_t, x_t-1)=rho=0.3. The mean of the first 70 data (x_1 to x_70) is 0 and the mean of the last 30 data (x_71 to x_100) is 2. Can I do it in the following way? x <- arima.sim(list=(ar=0.3), 100) mean <- c(rep(0, 70), rep(2, 30)) xnew <- x+mean If the above code...

Dear Friends, I want to know if any one could give me help on this. Suppose you want to simulate the p-value of a test statistic but you do not know how the distribution of the test statistic, how can one simulate the p-value. It will be helpful to get an example code in R on this Respectfully Submitted, Clement

Dear R-users, I'm trying to do some spatial simulation. I have two covariates, Z and C. I want to examine their relationship under different spatial distribution. I have no problem simulating completely spatial random process but I'm totally stuck on poisson (cluster) pattern. I already have a dataset with Z and C (obs=575) and I know the re...