similar to: FracSim set.seed

Hi all, I met with a problem about the weighted least square regression. 1. I simulated a Normal vector (sim1) with mean 425906 and standard deviation 40000. 2. I simulated a second Normal vector with conditional mean b1*sim1, where b1 is just a number I specified, and variance proportional to sim1. Precisely, the standard deviation is sqrt(sim1)*50. 3. Then I run a WLS regression without the

I want to load many files in the R. The names of the files are "Sim1.txt", " Sim2.txt", "Sim3.txt", "Sim4.txt", "Sim5.txt" and so on. Can I read them at one time? What should I do? I can give the same names in R. Thanks. For example: > tst=paste("Sim",1:20,".txt",sep="") # the file names > tst [1]

Hi, Can anyone help to write a more elegant version of my code? I am sure this can be put into a loop but I am having trouble creating the objects b1,b2,b3,...,etc. b1 <- rigamma(50,1,1) theta1 <- rgamma(50,0.5,(1/b1)) sim1 <- rpois(50,theta1) b2 <- rigamma(50,1,1) theta2 <- rgamma(50,0.5,(1/b2)) sim2 <- rpois(50,theta2) b3 <- rigamma(50,1,1) theta3 <-

Dear useRs, I want to simulate a time series (stationary; the distribution of values is skewed to the right; quite a few ARMA absolute standardized residuals above 2 - about 8% of them). Is this the right way to do it? #-------------------------------- load("rdtb") #the time series > summary(rdtb) Min. 1st Qu. Median Mean 3rd Qu. Max. -1.11800 -0.65010 -0.09091

2015-02-27 10:25 GMT-06:00 A J Stiles <asterisk_list at earthshod.co.uk>: > O.K. So what does your existing Dial() statement in extensions.conf look > like? > apology, put the gateway was sangoma but is a openvox , all my outgoing calls out for this context: [my-mobile-out] exten => _NXXXXXXX,n,Dial(SIP/1003/${EXTEN},55,rT) exten =>

On Friday 27 Feb 2015, ricky gutierrez wrote: > the problem is that my pbx all incoming calls using only the channel > gsm 1 , the idea is that an incoming call to channel 1 is passed to > channel 2 Ah. *Incoming* calls are not something that is within your control; they have already been routed onto a line by your telco. So you will need to speak to someone at your telco about doing

I don't understand what is happening! I have a (large) object sim1, an matrix list with dim c(101,101) where each element is an 3*3 matrix. I am subsetting that with a matrix coo, of dim c(100,2), of unique indices, but the resulting object has length 99, not 100 as expected. Code reproducing the problem follows: library(RandomFields) set.seed(123) sim0 <- GaussRF(x=seq(0, 100, by=1),

Dear all, My double for loop as follows, but it is little efficient, I hope all friends can give me a "vectorized" program to replace my code. thanks x: is a matrix 202*263, that is 202 samples, and 263 independent variables num.compd<-nrow(x); # number of compounds diss.all<-0 for( i in 1:num.compd) for (j in 1:num.compd) if (i!=j) { S1<-sum(x[i,]*x[j,])

To whoever this may concern I am trying to write a repeat loop and can't make out from the documentation on the website how exactly to construct the repeat, break structure of the loop. Below is the function sim2.dat that I am trying to create, in which firstly I create n random uniform(0,1) numbers. From there I assign values to k, P and F. From then on for each random number u[j] of u, I

Are there any packages that help with statistical analysis in situations where fractal geometry is relevant? Perhaps something that supports computation of fractal related statistics, such as the Hurst exponent or fractal dimension? Or perhaps which support obvious tasks such as taking samples from a dataset at different levels of granularity (such as sampling spatial data at cm, m, km

If the observations are normally distributed and the 2xk design is balanced, theory requires that the tests for interaction and row effects be independent. In my program, appended below, this would translate to cntT (approx)= cntR*cntI/N if all R routines were functioning correctly. They aren't. sim2=function(size,N,p){ cntR=0 cntC=0 cntI=0 cntT=0 cntP=0 for(i in 1:N){

hi R-listers, I would like some help recoding a variable. I have a dataframe 'cause' that translates between a set of codes: acc nds - - 1 2 3 4 5 8 ... ... the desired result for dataframe 'p': a - 1 5 5 would be: a b - - 1 2 5 8 5 9 I have tried: transform(p, b=cause$nds[cause$acc==p$a]) but for some reason it complains about the difference in length between the

Hi, Am using R 2.13.1 on Linux (Fedora). Is anybody else having problems with dev.copy2pdf xyplot output with the pch=1 (open circle) symbol? The symbols come out as "q" in the PDF. dev.copy2eps produces the correct results as does cairo_pdf. Other symbols produced with dev.copy2pdf seem ok. Thanks, Selwyn

Dear R People: At one time, there were packages called fractal and Fractal, respectively, which had different functions. I can't seem to find Fractal any more. Does it still exist somewhere, please? Thanks in advance for any help! Sincerely, Erin -- Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown mailto: erinm.hodgess

# How do I make from a SpatialGridDataframe a normal data frame? If I do this in package gstat: library(gstat) g.dummy <- gstat(formula = z~1, locations = ~x+y, dummy = TRUE, beta = 0, model = vgm(1,"Exp",15), nmax = 20) set.seed(1) iso <- predict(g.dummy, newdata = xy, nsim = 10) gridded(iso) = ~x+y # spplot(iso,main="isotropic") # not essential for this question # I

Greetings of Peace! I am new in R software. I want to use the correlation dimension corrDim in computing a map. I have installed the packages(fractal) in R. Now my problem is this, when I tried your example in the following site http://rgm2.lab.nig.ac.jp/RGM2/func.php?rd_id=fractal:corrDim  > library(fractal) > png(filename="corrDim_%03d_med.png", width=480, height=480) >

Dear all, I use Trixbox as my PBX. Until a couple of days I've installed a GSM Gateway to communicate with our three cellular phones: 15 64227777 15 64228888 15 64229999 The GSM Gateway has just one SIM. I use the Free PBX web interface in order to set up the route and trunk parameters: Trunk: ******* Name: SIM1 Peer details: host=10.10.1.2 (IP from GSM Gateway) port=5060 type=peer

I have a dataframe with many rows like this: > df X1 X2 X3 X4 X5 X6 X7 week d sim1 FALSE TRUE TRUE TRUE TRUE TRUE TRUE 1 0.3064985 sim1 is the rowname, X1..X7,week,d are the column names. X1..X7 are factors, booleans in this case. I need to add another row, represented by the following list: list(rep(T,7),5,0.0) -- i.e, TRUE in all boolean columns,

Dear R-users, I would like to show you a simple example that gives an overview of one of my current issue. Although my working setting implies a different parametric model (which cannot be framed in the glm), I guess that what I'll get from the following example it would help for the next steps. Anyway here it is. Firstly I simulated from a series of exposures, a series of deaths (given a

Hi all, I tried to use rsync-2.6.2 to sync a set of subdirectories inside a directory. While most of the files are being sync'd correctly, rsync seems to wrongly consider one file as uptodate. Both the source and destination have this file, with the same size and same date (according ls says). But md5sum shows that they are different. One possible solution is to enable the checksum, but it