similar to: Question about using “par”

Dear R users, I have various vector geometry operations to perform on 3-D coordinate data located on multiple (500+) csv files. The code I have written for the calculations works just fine. I have written a 'for' loop to automate the task of extracting the coordinates from the files and perform the analyses. The loop works reasonable well, but if the number of csv files is greater than

Dear All, I try to invoke a second program (called "gencont") within a R script. However this second program needs some time to finish. After finishing the output should be read into R and used for further calculation. Under windows, it perfectly works: result<-shell("gencont < yn.txt",intern=TRUE,wait=TRUE) # R waits untill the program finishes (apporx 3 min) and

I think the function you need is 'help.search'; try: help.search("binomial") and look for something obvious in the 'stats' package. A good deal quicker and easier than posting to an internet forum! Cheers, Mike. cathelf wrote: > > Dear all, > > I am trying a find the value "p" in binomial. > > X ~ Bin(n,p) > > I want to find the

Thank you for your guys reply for my previous question. But I got one more question about the boxplot. With the code in the R-help: boxplot(len ~ dose, data = ToothGrowth, boxwex = 0.25, at = 1:3 - 0.2, subset = supp == "VC", col = "yellow", main = "Guinea Pigs' Tooth Growth", xlab = "Vitamin C dose mg", ylab =

Hi, sorry for bothering your guys. I will trying to make some nice graph using boxplot. when I check the help file of boxplot, there is a sample code as: boxplot(len ~ dose, data = ToothGrowth, add = TRUE, boxwex = 0.25, at = 1:3 + 0.2, subset = supp == "OJ", col = "orange") legend(2, 9, c("Ascorbic acid", "Orange juice"),

I explain better: I have this percentage z<-1/(n*365) where n=20 so z is 0,000136986 Imagine: a matrix with length(n*365) so that matrix result is resultc( 1, 1-z result previous row-z result previous row-z .... 0) If I plot this is a linear decreasing line I want to plot an upper line wich is 10% higher tan resultc on the first values and then decreasing so that last value will be 0.05.

Hi, Dear R-users: Sorry for bothering your guys again. I think I should rewrite my question. I know how to extend whisker by using range. The question is that I will set the range=1.5, and at the same time, I only want to show the extreme outlier, like 0.01% and 99.99% percentile, so what should I do? Thank you very much! Catherine -- View this message in context:

Hi, Dear R-users, I have a problem when I drawing a boxplot. I want to extend the whisker to the 5% and the 95% quantiles and only show the most extreme outlier, like 0.01% and 99.99% percentiles. What should I do? I saw something on boxplot.stat, but even I define the parameter in boxplot.stat, what I should do next? how to connect it to boxplot? Thank you very much! Catherine -- View

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 to simulate 100 AR(1) data is right, what

Dear all, I have just sent an email with my problem, but I think no one can see the red part, beacuse it is black. So, i am writing again the codes: rm(list=ls()) #remove almost everything in the memory set.seed(180185) nsim <- 10 mresultx <- matrix(-99, nrow=1000, ncol=nsim) mresultb <- matrix(-99, nrow=1000, ncol=nsim) N <- 200 I <- 5 taus <- c(0.480:0.520) h <-

Hi, Can you tell me what is the meaning for "tail, 1" in "aggregate"? I also want to get some similar graph, but the data is not time series data. Suppose here is my data one, I want a graph with x-axis is just the index(1:9). The graph plot all the variable A, B,C,D. So there should be 4 lines for each graph. For the A line, at each time point, the letter A should be on

Hello. I am still a newbie in R. Excuse me if I am asking something obvious. My efforts to get an answer through browsing the mailing archives failed. I want to perform an augmented Dickey-Fuller test and to obtain AIC and BIC and to be able to impose some linear restrictions on the ADF regression so as to decide the correct order of autoregression. However I could find no obvious way to impose

Hi all I want to add a band of fluctuaci?n (exponential decreading) to a linear deacrecing values Imagine: I have a matrix like c(10,9,8,7,6,5,4,3,2,1) The thing is I want two new lines so that the m?ximum value of the new colum on the m?ximum is from 10% to 5% higher and the same lower for the m?nimum the final two matix will be something like c(10+0.10*10,9+0.089*9,8+0.075*8.....,1+0.05*1)

Is there a way to estimate a nominal response model? To be more specific let's say I want to calibrate: \pi_{v}(\theta_j)=\frac{e^{\xi_{v}+\lambda_{v}\theta_j}}{\sum_{h=1}^m e^{\xi_{h}+\lambda_{h}\theta_j}} Where $\theta_j$ is a the dependent variable and I need to estimate $\xi_{h}$ and $\lambda_{h}$ for $h \in {1...,m}$. Thank you, Mauricio Romero Quantil S.A.S. Cel: 3112231150

Hi, Sorry for bothering your guys. I am still working on boxplot. How can I extend the whisker to the 5% and the 95% quantiles? What I should do if I only want to show the most extreme outlier, like 0.01% and 99.99% percentiles? I saw something on boxplot.stat, but I even donot know how to connect it to boxplot? Thank you very much! suppose the follow is the data to make the boxplot. a <-