similar to: Simulation Function - Save results

hi R-help! I would like to know a simple and easy way (if posible) to do the following data manipulation. I have a matrix 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

I am trying to run this loop but it keeps coming up with the following error: 'Error: no function to return from, jumping to top level' n.simul <- 1000 A <- {for(i in 1:12){{ nT <- rpois(n.simul, T[i]) nP <- rpois(n.simul, P[i]) nD <- rpois(n.simul, D[i]) nC <- rbinom(n.simul, nT, C[i]) Score <- (5*nT + 2*nC + 3*nP + 3*nD) } return(Score)}} My aim is to calculate

Hello, I have a question on export and import data. Thank you for any suggestions. data 'simul' is generated as follows: N <- 20 n <- N/2 nsets <- 10 simul <- matrix(0,nsets,N) th <- c(0,1, 1) for(i in 1:nsets){ simul[i,] <- rnorm(N,mean= rep(th[1:2],N/2),sd=th[3]) } I exported data as follows: write.table(simul,

Hi List, I would like to plot multiple curves (parametric density curves) in one plot. For example: # parameters for three normal density curves parms = data.frame(ID=c(1,2,3),mu=c(50,55,60),sigma=c(10,12,15)) # I can easily draw three normal density curves using curve(): curve(dnorm(x,mean=parms$mu[1],sd=parms$sigma[1]),from=0, to=150, ylab="density", col="red")

I have a question about the lme function in R. My question is: After I got the object from function lme, why the numIter value of the object is always NULL? Following is my code: jjww<-lme(y~x*zz,data=simul,random=~x|group, control=lmeControl(returnObject=TRUE)) attributes(jjww) jjww$numIter the first 20 observation of data simul are: > simul y

Hi listers, I've been looking for a procedure, but I am not succeding... I have a function that give multiple results... Then, I would like to simulate this function n times, so I need to save/keep the n multiple results, in order to calculate my desired statistics... I have already tried with the RETURN and LIST FUNCTION, but I am not getting it right... An example of what I am looking for

Dear R-users, I would like to ask whether it's possible (for sure it would be), to plot each rows (or columns) in different graphs and in the same figure region without using the function "par" and then struggling around with "axes" and labels etc. Luckily, I would always have "rows + columns = even number" and the same "ylim". The next one could be a

Hello, I would like to add to axis labels special caracters. Instead of writing :plot(simul, xlab="beta", ylab="sigma11") It would be great if I could write something as in LaTeX :plot(simul, xlab="\beta", ylab="\sigma_{11}") Is there a way to do that ? Thank you _________________________________________________________________ ues clics pour retrouver tout

Hey gang, I've been able to use sipp to produce some call volume on our asterisk server. The server has no problems handling 50 simul calls. But then again, no RTP is being done. I tried to use the rtp echo ability of sipp but that doesn't seem to work right. I also setup a fake number in asterisk that when called by sipp, would dial another number via PRI, hoping that some 729

Dear useRs, I have written a function that implements a Bayesian method to compare a patient's score on two tasks with that of a small control group, as described in Crawford, J. and Garthwaite, P. (2007). Comparison of a single case to a control or normative sample in neuropsychology: Development of a bayesian approach. Cognitive Neuropsychology, 24(4):343?372. The function (see

Hello, i only got a small problem. i try to create automatic new dataframes, or png?s. the main problem i got is: how can i create automatic a new name for a file (read out by simply "for") - i tried to use "(paste...) but theres an errormessage, about a wrong declination. R told it is as.character, but need as.Real. Should i use another method than "paste"? i tried as

Hi I've two servers : 1. UK 2. Pakistan Pakistan * server has ISDN30. Pakistan(ISDN30) <====> UK ===> User Im planning to setup an IAX2 trunk between these two server ? so , how much bandwidth I need for 30 simul. calls ? Im planning to use G729 on both my server ? to support 30 calls over IAX2 do I've to change some setting during compile time or not ? pls suggest.

> Date: Thu, 10 Aug 2000 15:33:39 -0700 > From: Jesus Maria Frias Celayeta <iosu at bureau.ucc.ie> > > hi R-help! > > I would like to know a simple and easy way (if posible) to do the following > data manipulation. I have a matrix of experimental data (with replicae in time) actually, they seem to be data frames. > >

I was asked by my boss to do an analysis on a large data set, and I am trying to convince him to let me use R rather than SPSS. I think Sweave could make my life much much easier. To get me a little closer to this goal, I ran my analysis through R and SPSS and compared the resulting values. In all but one case, they were the same. Given the matrix [,1] [,2] [1,] 110 358 [2,] 71 312 [3,]

Hi Violeta, I compiled with clang -Oz and clang -Oz -mno-red-zone for comparisons against Oz and clang -O0 and clang -O0 -mno-red-zone for comparisons against a default clang. I unfortunately don’t have the clang commit I worked with on my home laptop, and don’t have access to the computer I was using at Apple, so I can’t help you there. Jessica > On Sep 28, 2016, at 11:22 AM, Violeta

Hi, I have a question about the lme function in R. My question is: After I got the object from function lme, why the numIter value of the object is always NULL? Following is my code: jjww<-lme(y~x*zz,data=simul,random=~x|group, control=lmeControl(returnObject=TRUE)) attributes(jjww) jjww$numIter the first 20 observation of data simul are: > simul y

Dear useRs, I'm happy to announce a substantial update of package actuar that bumps the version number to 2.0-0. This release focuses on additional support for continuous and discrete distributions, new functions to simulate data from compound models and mixtures, and revised and improved documentation. A slightly shortened version of the NEWS file follows: NEW FEATURES ? Support for the

Dear useRs, I'm happy to announce a substantial update of package actuar that bumps the version number to 2.0-0. This release focuses on additional support for continuous and discrete distributions, new functions to simulate data from compound models and mixtures, and revised and improved documentation. A slightly shortened version of the NEWS file follows: NEW FEATURES ? Support for the

i think your email at fisher is stripping the images On Nov 30, 2004, at 1:00 PM, asterisk-users-request@lists.digium.com wrote: > Send Asterisk-Users mailing list submissions to > asterisk-users@lists.digium.com > > To subscribe or unsubscribe via the World Wide Web, visit > http://lists.digium.com/mailman/listinfo/asterisk-users > or, via email, send a message with subject

Hey guys, I'm looking for some realworld specs on somebodys machine that will work with the Digium 4-port T1/PRI card and that will support 96 simultaneous calls. Dell is soon to release the PowerEdge 1850: 2U, Dual 3.6Ghz Xenon, 1Gb DDR2 RAM, Dual 36GB Ultra320 SCSI RAID, Hot swap Powersupply, one 64bit 133Mhz PCI and one 64bit 100Mhz PCI for about $3,000. Tack on a 4 port Digium card and