similar to: Save results

Hi, I have finally installed TeachingDemos trough zip file and installed on lbrary. One new question: I want to draw a line (like abline) on diagonal or the result of running a ols. How can I do it. On the other hand, I want to mark the length on each sub plot y axis. It is possible? My function is sca<-function(){ z2<-read.table("Ase.txt",header=T)

Sorry all for these two easy questions: First, I have a matrix with trhee columns: A= Name Eight Puntuation Pepe 1,85 10 Paco 1,7 7 Pablo 1,82 6 I want to scatter the two columns (I could use pairs or other functions...) but the only doubt is that I cannot find the way to add in the associated intersection (the point scattered) the label "Pepe" "Paco" and so

Hi I have an output data Data.csv this style: ,"V1","V2","V3" 1,"-9552","9552","C" 2,"0","9653","0" 3,"9614","9614","V" 4,"0","9527","0" 5,"-9752","9752","C" 6,"0","9883","0"

Hi all, I want to plot trough pairs() plot a matrix with 4 columns. I want to make a trhee plot in a graph. Plotting pairs colum 2,3,4 on y axis and 1 on X axis. You mean (a plot with three graphs) ommitting the first pair with itself. And only the pairs with colum 1 with the other not all pairs. I. e. this matrix 4177 289390 8740 17220 3907 301510 8530 17550 3975 316970 8640 17650 3651 364220

Hi all, I´m using titt<- expression(paste("R con ventanas de 200 ", italic(datos))) And it works properly on a plot(...,main=titt,..) But if I want to add an improvement to avoid typing n on the plot so that n<-200 titt<- expression(paste("R con ventanas de ",n, italic(datos))) It doesn´t recognize that n is a "variable" and adds "n" letter

Hi all I want to apply different colors on a simple plot: If I type par(br="gray") before a plot it puts all the image in gray but (imagine I run a simple plot) want to let the centrall box (where the dots are plotted) in white or image in lightblue. Can anyone guide me to apply this second step (make the box where the series are plotted in different colours). Thanks in advance.

Hi I wrotte this function but when I get the "tmp.xls" file it shows data in a rare way. I mean not appears a matrix with 2000 rows and 100 columns. Can anyone help me, guide me? kim<-function(){ tmp<-randz<-matrix(rnorm(200000, mean=0,sd=0.01 ),2000,100); dim(tmp) write(tmp,file="tmp.xls") Thanks in advance. On the other hand, everytime I execute a function

Hi all, I have two vectors (columns) called "efinal" and "efinal 2". I want to plot them on the same plot and "draw" a shaded area beween the two lines using function polygon I have tried all but I don ?t understand the polygon area, can you help me with examples? plot(efinal,type="l",ylim=range(min(efinal2),

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20091124/39eccb61/attachment-0001.pl>

Dear List, I just googled to find out if notched box plots are possible with ggplot2, but couldn't find a answer to it. boxplot() has the option: notch = TRUE, e.g.: boxplot(mpg$hwy, notch=TRUE) My example code (taken from the net) is: require(ggplot2) qplot(class, hwy, fill=factor(year), data=mpg, geom="boxplot", position="dodge")+theme_bw() Thank you for you help!

Hi I want to simulate random numbers normal distributed with this size (2000,10000). I tried this but my computer exhaust there is a fast way to make it? randz<-matrix(rnorm(2000000),2000,10000) [[alternative HTML version deleted]]

Hi all, Imagine I have a matrix and the first colum is a list that repeats the same names, I want to sum the second column on each unique name on first column. Imagine this: Pepe 2 Pepe 3 Pepe 4 Jose 2 Jose 5 Manuel 4 Manuel 2 I want to make a new matrix that calculates and recognizes that there are 3 different names ans sum second column. But a priori I don´t know the list of the different

Dear Colleagues, I hope to know how ordinal logistic regression with a mixed model is made in R. We (My colleague and I) are studying the behavior of a beetle. The attraction of beetles to a stimulus are recorded: the response is Slow, Mid, or Fast. They are based on the time after the presentation of the stimulus to the beetles. Because we do not observe the behavior continuously but do

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20091125/654989c6/attachment-0001.pl>

Hi all,,Im using function arima() I.e series is my data model<-arima(series,c(1,0,0)) forecast<-predict(model,80) I want to create a variable: b1<-forecast$pred - 1.96*forecast$se and save in a txt file but using this: save(b1, file= "b1.txt") creates afile butwith this inside: ]{HQ ~|LJIiW*-l)% )m#), RIiSfdbGur9p94; H"L#Rez 1y3pN8{,I6W!6= {6l?OMw_5KZ+ =

Hi all, Imagine I have a matrix G with N rows and M columns So L=NxM is the number of different cells in my matrix. I want to create a column vector F whose size will be F(L,1) So the fisrt row in F is G(1,1) Second row in F is G(1,2) When we arrive to a point M the element M+1 will be G(2, 1) Element M+2 will be G(2,2) and so on. I´m trying but allways error.... Easy Example: G= 2 3 4

Hi all I have installed quantmod package but when I try to obtain GOOG data appers this message: Can anyone inform why itappears? I type getSymbols("GOOG",src="google") Thanks and Best Regards for all Error en download.file(paste(google.URL, "q=", Symbols.name, "&startdate=", : no fue posible abrir la URL '

Dear List, I'm trying to solve this simple optimization problem in R. The parameters are the exponents to the matrix mm. The constraints specify that each row of the parameter matrix should sum to 1 and their product to 0. I don't understand why the constraints are not satisfied at the solution. I must be misinterpreting how to specify the constrains somehow. library(alabama) ff <-

Hi all, I have a matrix dd with dates and data. i WANT to create a grpah with shade between the line and the zero axis. I?m trying to use polygon but something doesn? t work Imagine. I want to plot g<-read.table("dd.txt", col.names=c("fecha","DP")) g$fecha <- as.Date(g$fecha, format="%d/%m/%Y") t<-g$fecha st<-length(g$DP)

Hello everyone, I want to calculate type I error rate for modified F statistic for one way robust anova. I need to find the j group trimmed mean and winsorized sum of squared deviations. Here I attached my code for j=2 to make it simple. Originally I have j=4. Hope you can help. I need to run it for 1000 times My problem is: i) the value of F-test obtain from my simulation below is in negative