similar to: A question on type="h" plot lines

Dear R users, I need to to the the following. Let a= 1 2 3 4 5 6 and b= -1 -2 -3 be (2x3) matrices. -4 -5 -6 I need to combine the two matrices into a new (2x6) matrix like this: ab = ( 1 -1 2 -2 3 -3 ) 4 -4 5 -5 6 -6 How can this be done in R? ----------------------------------------------------------------- ?????? ???

Hello, Can anyone suggest a simple way to generate a Kaplan-Meier plot with 2 survfit objects, just like this one:? https://drive.google.com/file/d/1fEcpdIdE2xYtA6LBQN9ck3JkL6-goabX/view?usp=sharing Suppose I have 2 survfit objects: fit1 is for the curve on the left (survtime has been truncated to the cutoff line: year 5), fit2 is for the curve on the right (minimum survival time is at the

Hello! After I type "upsmon", in syslog the following message appears: Aug 16 14:11:09 marto upsmon[7731]: Master privileges unavailable on UPS [Inform@localhost] Aug 16 14:11:09 marto upsmon[7731]: Reason: Access denied Can you tell me where my error is, my configuration files are: ups.conf: [Inform] driver = powercom port = /dev/ttyS0 desc = "inform guard" linevoltage = 220

Dear R users, I am using the vars package to calculate the impulse response functions and the forecast error variance decomposition of a VAR model. Unfortunately I do not know whether these functions assume unit or one standard deviation shocks. I tried to look into the code of these functions, but in vain: neither irf, nor vars::irf, nor vars:::irf output the code of the functions. Does someone

Dear R users, I need perform structural analysis on a no intercept VAR model. Unfortunately the functions irf.VAR and dfev that come with the MSBVAR package only work with objects output by the reduced.form.var function, which seems to only evaluate VAR models with intercept. Is there a way to suppress the estimation of intercept term in reduced.form.var? Do I need to modify the code, and if I

Dear R users, I need to plot to perpendicular straight lines. However, although I set the coefficients so that the lines are perpendicular, they do not look to be so in the plot. Here is a minimal working example: plot(x=c(-1, 1), y=c(-1, 1)); abline(a=0, b=1/sqrt(2)); abline(a=0, b=-1/sqrt(2)) Please tell me if the same problem is valid by you. I am running R-2.10.1 on Linux. Is there a way

Dear R users, I have a lot of experience with traditional R graphics, but I decided to turn to trellis as it was recommended for spatial graphs by the sp package. In traditional R graphics I always first set the size of the device region absolute units (e.g. mm) and then I firmly fix the inner margins with mai and the outer margins with oma also in absolute units. What is left from the device

Thank you for your response. The point is that I have configured my system to shut down at the event "upsc Inform@localhost ups.status"!="OL", not at "LB". So there will be no /etc/killpower file created. At an ups.status != "OL" event, the following script named "notifyme" is run: echo `date` >> /root/upsLog echo " No power, shutting

Dear nut users, I am running nut-2.6.0 on Slackware-Linux-13.37.0 with kernel 2.6.37.6. I am trying to get the software work for a repotec UPS with model name: RPT-1003AU. The UPS communicates with the computer via USB. I know that the model is not officially supported but I want to try out whether it will work with some of the repotec drivers. Here is the result with the genericups upstype=13

Hello! I need to perform weighted multidimensional scaling analysis(WMDS). I did rsitesearch, googled, but I could find no info on how to perform WMDS using R. In several places they say it is possible with the ALSCAL algorithm, but I could not find the relevant function to carry it out. ----------------------------------------------------------------- ???????? ?? ??????? ?? ??? ???????????!

Dear forum, I need to create biplot graphic. I took a look at the code of the biplot.default function from the stats package and I found that the arrows for the variables are multiplied by a coefficient of 0.8 before plotting: ... if (var.axes) arrows(0, 0, y[, 1L] * 0.8, y[, 2L] * 0.8, col = col[2L], length = arrow.len) ... I have no idea why this is done. Could somebody

Dear R users, I am running R-2.15.1 in Linux Slackware64-14.0. Here is my minimal working example: testfun <- function (x) { a <- 0; sapply(X="a", FUN=get, envir=sys.frame(which=x)); } Inside R, that is R called from within a Linux terminal, the following code works: testfun(x=5) print(testfun(x=6)) But within rkward the above code fails and the following works: testfun(x=1)

hi R gurus I saw some graphs with vertical band like this one: http://pragcap.com/wp-content/uploads/2010/04/GS.png how to draw the blue band in R, can't find any clue to do this,any ideas? thanks in advance [[alternative HTML version deleted]]

Dear All, I am running Slackware64-14.1(Linux marto 3.10.17 #1 SMP Fri Feb 14 16:39:21 CST 2014 x86_64 Intel(R) Pentium(R) CPU G3250 @ 3.20GHz GenuineIntel GNU/Linux) with kernel 3.10.17 and nut package nut-2.6.5-x86_64-1_SBo, buit from source using the build script http://slackbuilds.org/repository/14.1/system/nut/ In the README.Slackware it is suggested that the commands for killing the

Hello list, I've been working on this problem for a while and I haven't been able to come up with a solution. I have a couple of functions that plot a bunch of data, then a single point on top of it. What I want is to be able to change the plot of the point without replotting all the data. Consider the following example: x = rnorm(100,1,0.5) y = rnorm(100,1,0.5) plot(x,y,pch=16)

Hi all, I would like to ask your help to add a color band (Ι am not sure regarding the right term, this color band at the right of the plot "describing" values with their corresponding color. For now I have only this code test<-matrix(data=runif(10000),nrow=100) plot(test,axes="FALSE") axis(1,at=c(0,1),labels=c("a","b")) # but  I would like to add

Hi. I'm trying to add some vertical lines to an xyplot in which x axis is a temporal variable (class=dates) and y axis a factor variable. This is my code: xyplot(factor(Abacus$Emisor)~Abacus$Dia,xlab="Date", ylab="Fish",pch=124, scales=list(x=list(format="%b %Y",tick.number=30)) , panel=function(...){ panel.xyplot(...) panel.abline(v=08/08/2008) }) I'd

Hi, I try to create a vertical line in my plot, which has a xaxis comprising time formated data. This is what I tried: ---------------- y<-152833 x<-strptime(y, format="%H%M%S") abline(v=x, col="red") ---------------- for some reason, it doesn't work and no error msg is displayed... I hope someone can tell me what I do wrong, or maybe an alternative. Thanks

Hi, I have a question for the group, perhaps someone can help me figure this out. I've already looked in the help files and they were no help to me. I have a vector of values and I am plotting an ecdf graph. 1. How can i draw a continuous line through the ecdf points? (lines and type for the plot with an ecdf object does not work) 2. Supposing I have this line drawn. I can add a

I have a problem with the ylabs in barplot. When I draw the bars horizontal the text for the y axis is not drawn horizontal too. The text remains vertical. How can I change that. Thank's for your advice. There is an example below. ##Example plot.new() jib <- data.frame(c(1:15),c(11:25),c(15:1),c(25:11)) colnames(jib) <-