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

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 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).

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 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),

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

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

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

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