[R] labels outliers in boxplot

2007-01-12 Thread antoniababe
Dear R-users,

Following is part of my data, where slide has 36
levels and block 48 levels. I have done boxplot for
each slide on the same graph. There are outliers for
each slide and I tried to use indentify functtion to
identify outliers in such a way that when I click on
an outlier or point, the points will be labelled by
either their block or ID or by both but without
success. How can I make it work or are there other
ways to do it than using identify function?


Thanks in advance, 

Jenny,

dat1[1:10,]
 y Slide Block  ID Control
1   0.03147823 1 1   IgG-human   5
2  -0.23815974 1 1 LPPAANDVSVLTAAR   0
3  -0.71926359 1 1 HTKHYRVVSKPAALV   0
4  -0.14607826 1 1 FVALPAATADAYATT   0
5   0.89553073 1 1 NYPAMMAHAGDMAGY   0
6  -0.67587100 1 1 RRALRQIGVLERPVG   0
7   0.32636034 1 1 DCGTIRVGSFRGRWL   0
8  -1.44057259 1 1 MAKLSTDELLDAFKE   0
9  -0.37064338 1 1 LELSDFVKKFEETFE   0
10 -0.20387233 1 1 VSRRAKVDVLIVHTT   0


 tb_ncs-subset(dat1,dat1$Control==1) ### this
data contains only negative controls

   par(las=2,mar=c(10.1,4.1,4.1,2.1)) 
  
boxplot(split(tb_ncs$y,tb_ncs$Slide),col=orange,
cex=.65,
   outline=TRUE,main=Negative control
response of each patient, cex.main=1, font.main=1,
col.main=blue,
names=c(1:35,B))
   grid(nx=NA, ny=NULL)   
  ### grid over boxplot 
   legend(bottomright, B = Buffer +
sec,text.col=blue)
   

out.block-
identify(tb_ncs$y,tb_ncs$Slide) 






_
Flyger tiden iväg? Fånga dagen med Yahoo! Mails inbyggda

__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] labels outliers in boxplot

2007-01-12 Thread antoniababe
Dear talepande,

Thanks for your suggestion,  
I have already tried to use it, but the identify
function gave me only  the observation number
everytime I clicked on any point.What I want is
instead of obervation numbers it would be block and/or
slide numbers.

Any other idea how I can make it works ?

Thanks 
--- talepanda [EMAIL PROTECTED] skrev:

 because given data is a part of your data, I cannot
 examine,
 however, try:
 
 ##out.block-identify(tb_ncs$y,tb_ncs$Slide)
 out.block-identify(tb_ncs$Slide,tb_ncs$y)
 
 
 On 1/11/07, [EMAIL PROTECTED]
 [EMAIL PROTECTED] wrote:
  Dear R-users,
 
  Following is part of my data, where slide has 36
  levels and block 48 levels. I have done boxplot
 for
  each slide on the same graph. There are outliers
 for
  each slide and I tried to use indentify functtion
 to
  identify outliers in such a way that when I click
 on
  an outlier or point, the points will be labelled
 by
  either their block or ID or by both but without
  success. How can I make it work or are there other
  ways to do it than using identify function?
 
 
  Thanks in advance,
 
  
 
  dat1[1:10,]
  y Slide Block  ID Control
  1   0.03147823 1 1   IgG-human   5
  2  -0.23815974 1 1 LPPAANDVSVLTAAR   0
  3  -0.71926359 1 1 HTKHYRVVSKPAALV   0
  4  -0.14607826 1 1 FVALPAATADAYATT   0
  5   0.89553073 1 1 NYPAMMAHAGDMAGY   0
  6  -0.67587100 1 1 RRALRQIGVLERPVG   0
  7   0.32636034 1 1 DCGTIRVGSFRGRWL   0
  8  -1.44057259 1 1 MAKLSTDELLDAFKE   0
  9  -0.37064338 1 1 LELSDFVKKFEETFE   0
  10 -0.20387233 1 1 VSRRAKVDVLIVHTT   0
 
 
   tb_ncs-subset(dat1,dat1$Control==1) ### this
  data contains only negative controls
 
par(las=2,mar=c(10.1,4.1,4.1,2.1))
 
  boxplot(split(tb_ncs$y,tb_ncs$Slide),col=orange,
  cex=.65,
outline=TRUE,main=Negative
 control
  response of each patient, cex.main=1,
 font.main=1,
 col.main=blue,
  names=c(1:35,B))
grid(nx=NA, ny=NULL)
   ### grid over boxplot
legend(bottomright, B = Buffer +
  sec,text.col=blue)
 
 
 out.block-
  identify(tb_ncs$y,tb_ncs$Slide)
 
 
 
 
 
 
 

_
  Flyger tiden iväg? Fånga dagen med Yahoo! Mails
 inbyggda
 
  __
  R-help@stat.math.ethz.ch mailing list
  https://stat.ethz.ch/mailman/listinfo/r-help
  PLEASE do read the posting guide
 http://www.R-project.org/posting-guide.html
  and provide commented, minimal, self-contained,
 reproducible code.
 
 






_
Flyger tiden iväg? Fånga dagen med Yahoo! Mails inbyggda

__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] Boxplot issue

2007-01-08 Thread antoniababe
Dear R-users,

I have a data frame containing 2 colums: column 1 is
the patient numbers (totally 36 patients), column 2 is
patient's response values (each patient has 100
response values). If I produce a boxplot for each
patient on the same graph in order to compare them
against each other then the boxplots are very small. 

How can I instead of creating one graph containing 36
boxplots, create four different graphs where three of
them have 10 boxplots each representing data of 10
patients and the fourth graph has boxplots of
remaining patients ?

Thanks alot for any suggestion,
Greetings,
Antonia

__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] different points and lines on the same plot

2007-01-07 Thread antoniababe
Dear all,

I have following data called paitent

daypatient1patient4patient5patient6   
0   -0.27842688 -0.04080808 -0.41948398 -0.04508318 
56  -0.22275425 -0.01767067 -0.30977249 -0.03168185  
112 -0.08217659 -0.26209243 -0.29141451 -0.09876170 
252  0.08044537 -0.26701769  0.05727087 -0.09663701 

where each patient have response values at four time
points. I want to plot each patient's values over time
on the same plot where the value points are connected
by line. That is, the graph will have four lines for
the four patients. I tried the program below but
couldn't make it work correctly. I'm new beginner and
haven't yet learned how functions line and points work
together. Hope you can help me out.

Thanks for your help,

Antonia

par(mfrow=c(1,1))
plot(patient[,1],patient[,2], pch=1,
type=l,col=1,cex=1,lwd=2,
 xlab=Days, ylab=Patient
response,cex.main =1,font.main= 1,
 main=NULL)
   
 
points(patient[,1],patient[,3],col=2,pch=1,cex=1)
lines(patient[,1],patient[,3],col=2,lty=1,cex=1)
   
points(patient[,1],patient[,4],col=3,pch=1,cex=1)
lines(patient[,1],patient[,4],col=2,lty=2,cex=1)
   
points(patient[,1],patient[,5],col=4,pch=1,cex=1)
lines(patient[,1],patient[,5],col=2,lty=1,cex=1)
  points(patient[,1],patient[,6],col=5,pch=1,cex=1)
 lines(patient[,1],patient[,6],col=2,lty=1,cex=1)

__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] new R-user needs help

2006-10-18 Thread antoniababe
Dear all,

I have a dataset 

Exp name id x1 x2 x3 x4 x5 
1   a1  23 24 23 22 30
1   b2  25 26 27 23 24
1   c3  32 19 23 25 28
2   a4  28 32 24 26 27
2   b5  23 24 25 26 28
2   c6  23 31 30 38 23
3   a7  24 25 31 27 29
3   b8  28 25 26 32 28 
3   c9  21 31 28 23 29

where each name a,b and c is repeated three times.
Now I want to calculate the the averages of a,b and c
over the three experiments for each x. For example for
x1 I want to calculate the average of a like

[a(from exp1) + a(exp2)+ a(exp3)]/3= (23+28+24)/3

How can I do that ?

Thanks for your help,
Regards,

Antonia

__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.