similar to: scatterplot in CAR

Hi All. I'd like to change the default plotting colors used to construct a scatterplot with regression line in the CAR package. So, scatterplot(y~pred,smooth=FALSE, xlab="X", ylab="Y", lwd=2) If I change the palette (e.g., palette(ranbow(6)), I can change the color of the lines and points. However, the axes and labels remain in black (i.e., the first color in the

Network Blitz I'm trying to generate a graph to summarize Interest Rate Parity. This involves a scatterplot of x against y where the x and y limits are set so to center the graph on 0,0 and then adding each axis line and a 45 degree line. Using CAR's scatterplot (sample code below) everything except the axes plot fine: scatterplot( Interest.Rate.Dif ~ ForPrm| profit,

Hello all! System: R2.0.1, W2k. All packages apdated with update.packages(). I use the default graphic device for my system (no active choice). I?m trying to fine-tune scatterplot() graphs (package car), for publication. I want to control plotting symbols and colours to match plots from other, less flexible, systems. I also need to make all text and symbols bigger, to enable printing of the

I am getting an error message from scatterplot: > library(car) > scatterplot(Prestige$income~Prestige$type) Error in Summary.factor(c(2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, : range not meaningful for factors In addition: Warning message: In Ops.factor(x[floor(d)], x[ceiling(d)]) : + not meaningful for factors > The command does output the kind of graph that I want (boxplots).

Hi Folks, Please, when I ask the option reg.line at the scatterplot in package car, the OLS models includes a constant? If not how can I do it sing the following code: scatterplot(lfirms ~ lscale, data=dataset, reg.line=lm, smooth=FALSE, labels=FALSE, span=0.5, xlab="Relative Plant Fixed Cost", ylab="Relative Number of Firms", pch=c(18),

Hello R folks, I am encountering a problem with the following scatterplot function from the car package: > scatterplot(y~x|z) where y and x are continuous (interval) random variables and z is a categorical variable. When z is a categorical variable coded 1 or 2, I (appropriately) get a scatterplot of y by x, coded by z. Similarly, when z is a categorical variable coded 1, 2, or 3, there is

I am trying to fit a logistic regression model for this data set. Firstly, I want to plot P(Y=1) vs As and P(Y=1) vs Aa. Does any body know how to do these in R. Thanks, Aimin > p5 <- read.csv("http://www.public.iastate.edu/~aiminy/data/p_5_2.csv") > str(p5) 'data.frame': 1030 obs. of 6 variables: $ P : Factor w/ 5 levels "821p","8ABP",..: 1

I have a dataset like this: p aa index x y z sdx sdy sdz delta as ms cur sc 1 821p MET 1 -5.09688 32.8830 -5.857620 1.478200 1.73998 0.825778 13.7883 126.91 92.37 -0.1320180 111.0990 2 821p THR 2 -4.07357 28.6881 -4.838430 0.597674 1.37860 1.165780 13.7207 64.09 50.72 -0.0977129 98.5319 3 821p GLU 3 -5.86733 30.4759

scatterplot() is autogenerating a regression and loess line. I want to remove the loess line, but can only find samples for how to add it. I definitely do not have code that adds it (see below), so I am puzzling... scatterplot(Aroused ~ anxious | group, data=data, xlab="Aroused", ylab="Anxious", main="Arousal and anxiety ratings", xlim=c(0.5, 9.5),

hi, new to R and using the car package to do some scatterplots with ellipses hoping to add the area and center points of each ellipse to the legend? looking for some direction / ideas.... here is the script, the data is where golf shots end up by club. x (dispersion), y (distance), group (Club) thanks, sam ## usage: Rscript shotScatter.R infile outfile level1 level2 minX maxX minY maxY

Hi, using scatterplot in 'car' package. There are 4 plots to be pictured together (ignore data used): library(car) par(mfrow=c(3,1)) ts.plot(rnorm(100)) ts.plot(rnorm(100)) scatterplot(rnorm(100),rnorm(100),boxplot="",smooth=T) # scatter with non-parametric fit problem is that the last one comes 'always' separate. Is there any way to solve it?? of course can use

Hello, I use nnet to do prediction for a continuous variable. after that, I calculate correlation coefficient between predicted value and real observation. I run my code(see following) several time, but I get different correlation coefficient each time. Anyone know why? In addition, How to calculate prediction accuracy for prediction of continuous variable? Aimin thanks, > m.nn.omega

I have a data(ABC) like this: x y A 3 4 A 1 3 B 2 6 B 4 8 C 5 4 C 6 7 I want to split this data into A: x y A 3 4 A 1 3 B B 2 6 B 4 8 C C 5 4 C 6 7 anyone knows how to do that? thanks, Aimin Yan

I have a set data={A,B,C,D,E,F,G} I want to choose 2 letter from 8 letters, i.e. generate the combination set for choose 2 letters from 8 letters. I want to get the liking: combination set={AB,AC,AD,....} Does anyone konw how to do in R. thanks, Aimin

I have a question in R. In directory H:/Delta_angle I have 19 files like this: ALA.delta ASN.delta ASP.delta CYS.delta GLN.delta GLU.delta HIS.delta ILE.delta LEU.delta LYS.delta MET.delta PHE.delta PRO.delta SER.delta THR.delta TRP.delta TYR.delta VAL.delta I want to read these files to 19 data sets in R. All these data sets have "P","AA",index","delta"

I try Sweave and get Sweave-test-1.tex but hot to run LaTeX on 'Sweave-test-1.tex'? I am using WinEdt. thanks, Aimin > Sweave(testfile) Writing to file Sweave-test-1.tex Processing code chunks ... 1 : print term verbatim 2 : term hide 3 : echo print term verbatim 4 : term verbatim 5 : echo term verbatim 6 : echo term verbatim eps pdf 7 : echo term verbatim eps pdf

I have a data set like this I want to assign "outward" to Y if sc <90 and assign "inward" to Y if sc>=90. then cbind(p1982,Y) to get like these p aa as ms cur sc Y 1 154l_aa ARG 152.04 108.83 -0.1020140 92.10410 inward 2 154l_aa THR 15.86 28.32 0.2563560 103.67100 inward 3 154l_aa ASP 65.13 59.16 0.0312137 7.27311 outward 4 154l_aa CYS 57.20 49.85

Dear R list, I have a question in R, it could be very simple, but I don't know how to do it? for example: I assign 6 to x in beginning of of my R script code > x<-6 ...... After many line code, I forget using x variable before, I use x again, and do assignment like this > x<-45 > x [1] 45 then value 6 of previous x is replaced by 45. I am wondering if there are some way R

R version 2.11.1 on WinXP How do I get 3 scatterplots with marginal boxplots (from the car package) onto a single plot? I have a data frame called bank > dim(bank) [1] 46 5 head(bank) x1 x2 x3 x4 pop 1 -0.45 -0.41 1.09 0.45 0 2 -0.56 -0.31 1.51 0.16 0 3 0.06 0.02 1.01 0.40 0 4 -0.07 -0.09 1.45 0.26 0 5 -0.10 -0.09 1.56 0.67 0 6 -0.14 -0.07 0.71 0.28 0

How to add mean,sd, number of observation in each panel for lattice histogram? Aimin