search for: simu

Dear all, I need to read in 4 sets of tab delimited data in a loop. The 4 data sets are called "simu1.dat", "simu2.dat" and so on. I know what I need on the righthand side of the read.table expression but I can't the left hand side of it to work (see the line in bold below). Can you kindly help? Many thanks. simu1 <- matrix(0,30,3) simu2 <- matrix(0,30,3) simu3 <- mat...

It runs more than twice as slowly using 8 core than using a single core in inverting large matrix. Tested on 8 core Windows XP 64 machine. > n = 1000 > n.simu = 100 > func1 = function() + { + x = rnorm(n*n) + dim(x)=c(n,n) + y = solve(x) + } > > setMKLthreads(1) > system.time(for(i in 1:n.simu) func1()) user system elapsed 69.48 2.42 71.91 > > setMKLthreads(8) > system.time(for(i in 1:n.simu) func1...

...nd saved it in c:\Perl directory. My path is set as: C:\bin;C:\MinGW\bin;C:\Perl\bin;%SystemRoot%\system32;%SystemRoot%;%SystemRo ot%\System32\Wbem;C:\Program Files\Common Files\Adaptec Shared\System;c:\matlab6p5\bin\win32; When I input the following command: c:\rw1071\bin> Rcmd SHLIB -o simu.dll simu.c I met the following problem: 'make' is not recognized as an internal or external command, operable program or batch file. Could any one tell me the problem here? which steps I have lost or the path has not been correctly set??? Thank you for your help in advance! Ping

...eck: PAM: Account Validation Failed - Rejecting User ws13! Nothing interesting in auth.log and the same message in syslog (where slapd logs to). I don't know if this is a samba issue or ldap or network... Any suggestions are highly welcome. We've got 20+ clients and users typically log in simultaneously. timbo smb.conf: panic action = /usr/share/samba/panic-action %d dos charset = 850 unix charset = ISO-8859-15 display charset = ISO-8859-15 netbios name = KRAKE workgroup = GHSWA hosts allow = 192.168. inherit acls = yes update encrypted = yes obey pa...

Does anybody have any suggestions on how I might simulate from fitted GAM model? I am using the gam function in the mgcv package to fit a variable coefficient model like the following from the examples. I would like simulate based on the fitted model like the simulate function in the stats package does for lm models. library(mgcv) set.seed(10) ##...

...u.mean.initial); b <- b - solve(obj$Hessian, obj$score); D.u <- obj$uu + TT; u.mean.initial <- obj$u.mean.initial; print(i); print(date()); print("inside REML"); print(TT); if(i==n.iter) write(TT, file="c:/liao/reml/simu50_8.dat", append=T); print(D.u); } list(b=b, D.u=D.u); } bias <- function(b, D.u, x, z) { yy <- glmm.sample.y(b, D.u, x, z); #this is the sampling stage obj1 <- mle(b, D.u, x, yy, z, F, F, n.iter=50); obj2 <- mle(b, D.u, x,...

hi. i did some research first, but didn't find what i was looking for... the thing is: i generated data with correlated errors and simulated the power with using aov(). what i wanna do now is something similar while using lme(), so that the corr-structure will be paid attention to. i'm not quite sure how to do that... i hope someone can help me. thanks in advance. mel code: ee<-mvrnorm(n=100, mu=rep(0,N),Omega.block) #...

...e function following the instruction in the Pinheiro and Bates (2000), but I could not make it work. Then I used varComb in gls with weights=varComb(varExp(form=~W), varExp(form=~I(W^2))))). But the estimated variance parameters seems to have a large discrepancy from the true values (I used the simulated data). This makes me wonder if it is a right way to model variance function exp(c0+nu0*W +nu1*W^2) using "varComb". The codes and outputs are copied below. Any suggestions and help are very apprecited library(nlme) simulate.pilot = function(m, mn, sigma, alpha0, al...

...ce function following the instruction in the Pinheiro and Bates (2000), but I could not make it work. Then I used varComb in gls with weights=varComb(varExp(form=~W), varExp(form=~I(W^2))))). But the estimated variance parameters seems to have a large discrepancy from the true values (I used the simulated data). This makes me wonder if it is a right way to model variance function exp(c0+nu0*W +nu1*W^2) using "varComb". The codes and outputs are copied below. Any suggestions and help are very apprecited library(nlme) simulate.pilot = function(m, mn, sigma, alpha0, alpha1, c0...

Hi there, I have some question about how to use this package: My code structure is as follows: rm(list=ls()) set.seed(2014) library(snowfall) ############# #part 1: generate the initial value for simulation #part 2: some functions will be used in the simulation #simulation part, start from "while loop" sink('SnowFallExample.Rout', split=TRUE) .Platform .Machine R.version Sys.info() sfInit(parallel=TRUE,cpus=4,type="SOCK") simu<-1 while(simu<=1000){ function...

...haracteristics: Intel Core 2 duo 2.33Ghz and 3Gb of RAM memory (with Windows XP Pro SP2) and when I am running the same R syntax I receive the following error message: "Error: cannot allocate vector of size 216.0 Mb". This is an extract of my R syntax: > memory.limit(4000) NULL > simu.BYM<-bugs(data, inits, parameters, model.file="Modelo.txt", n.chains=3, n.iter=200000, n.burnin=20000, n.thin=180, debug=FALSE, DIC=FALSE, digits=5, codaPkg=FALSE, bugs.directory="c:/Archivos de programa/WinBUGS14", working.directory=NULL) Error: cannot allocate vector of siz...

Hello, everyone! I picked up R several months ago and have adopted it as my choice for statistical programming. Coming from a Java background, I can honestly say that R is not only free, it is better tha S-plus: the lexical scope in R makes it very simple to simulate Java's object model. For this, I encourage everyone to read the artcle: Robert Gentleman and Ross Ihaka (2000) "Lexical Scope and Statistical Computing", Journal of Computational and Graphical Statistics, 9, 491-508. I have coded a simulation program which runs at 1/5 the speed...

Anyone using the rrdgraph report on the puppetmaster on Debian stable? I''m trying to get this working and not getting far. I even took rrdtool and librrd0-dev from backports and I still get the same problem. It seems that the rrdfiles are being somewhat created, but nevertheless, the puppetmaster isn''t happy and keeps spewing: err: Report rrdgraph failed: Could not create

Dear all, I am doing 200 times simulation. For each time, I generate a matrix and define some function on this matrix to get a 6 dimension vector as my results. As the loop should be slow, I generate 200 matrice first, and save them into a list named ma, then I define zz<-sapply(ma, myfunction) To my surprise, It almost costs...

Dear all, I am doing 200 times simulation. For each time, I generate a matrix and define some function on this matrix to get a 6 dimension vector as my results. As the loop should be slow, I generate 200 matrice first, and save them into a list named ma, then I define zz<-sapply(ma, myfunction) To my surprise, It almost costs...

Hi there, me again! I´m looking for a solution for the following: I have two (a bunch really..) of services that work together and have to exchange some configinformation between each other. I can get this to work in only one of three cases: 1. Works # class serviceB { file {"serviceBconfig": content => $serviceA::serviceAvar } } # # class serviceA { $serviceAvar

Hola, ¿qué tal? Tengo una pregunta de esta

...g Splus to run a multiphase mixed-effects model. The quations of the models are as below: gf[ij]=b0[i]+b1[i]*age[ij]+b2[i]*max(0,(age[ij]-tau[i]))^2+e[ij] b0[i]=b00+e[i0] b1[i]=b10+e[i1] b2[i]=b20+e[i2] tau[i]=tau+e[i3] i: 1,2,...,100 subjects j: 1,2,...,6 occasions The main scripts of Splus is: simu1<-groupedData(gf~age|id) simu.nlme<-nlme(gf~(b0 + b1 * age + b2 * max(0,(age-tau))^2),data=simu1,fixed=list(b0~1,b1~1,b2~1,tau~1), + random=list(b0~1,b1~1,b2~1,tau~1),start=c(b0=4,b1=5.32,b2=-5.29,tau=14.8)) But it seems like that there are some errors in it. Error in .C("mixed_EM&qu...

Hi, Sorry to bother those who are not interested in the question. I am using R to run a multiphase mixed-effects model. I simulated a data set by using the following model. And now I want to use R to estimate the parameters and compare the results with the true values. The equations of the models are as below: gf[ij]=b0[i]+b1[i]*age[ij]+b2[i]*max(0,(age[ij]-tau[i]))+e[ij] b0[i]=b00+e[i0] b1[i]=b10+e[i1] b2[i]=b20+e[i2]...

Hi everyone, I have posted a similar question to this list, but I don't get a reply. I really want to solve this problem, so I post it again... I am trying to use R to fit some mixed-effects models for a nested data. The data is a simulated data with 111 subjects. Each subject has 6 waves' data. Below are the first two subjects' data : > simu1[1:12,] Grouped Data: gf ~ age | id id occasion age gf tau b0 b1 b2 1 1 1 4.056714 35.08943 14.8 14.502078 5.209333 -5.199...