"Deepayan Sarkar" <[EMAIL PROTECTED]> writes:
> On 7/30/07, Patrick Drechsler <[EMAIL PROTECTED]> wrote:
>> what is the recommended way of adding a strip to a lattice plot?
>
> See ?strip.default.
>
>> In the example below I would like to add the value of mean(y) to a new
>> strip.:
>>
>> --8<---
I have to compare four different grape varieties proteome in two different
years. I don't know what test would be more suitable for my data. I think that
an anova two way can be usefull also if someone suggested me to perform a
manova. In addiction, I can perform each test on a single protein a
I have to compare four different grape varieties proteome in two different
years. I don't know what test would be more suitable for my data. I think that
an anova two way can be usefull also if someone suggested me to perform a
manova. In addiction, I can perform each test on a single protein a
p <- seq(0.001,0.999,,1000)
x <- qt(p,df=9)
y <- dt(x,df=9)
plot(x,y,type="l")
polygon(x=c(x,rev(x)),y=c(y,rep(0,length(y))),col="gray90")
Hope this helps.
ST
--- "Nair, Murlidharan T" <[EMAIL PROTECTED]> wrote:
> Indeed, this is what I wanted, I figured it from the function you and
> Mark poi
Sehr geehrte Damen und Herren.
Air Turkey ist in Istanbul-Izmir-Antalya ein
Spezialveranstalter.www.air-turkey.com
Air Turkey spezialisiert für ihre
Fluege-Chartermiete,Spezialreisen,Hotels-Gruppenreisen von Kultur und
Bildungsreisen.Fuer Vereinen,Sportgruppen,Kirchen.
Wir sind spezialisiert a
Hi Conor,
I hope I interpreted your question correctly. I think for the first one you
are looking for a conditioning plot? I am going to create and use some
nonsensical data - 'iris' comes with R so this should be reproducible on your
machine:
library(lattice)
data(iris)
x <- iris
# make some fac
Hello,
Can anyone explain the following behaviour? To me it seems a bug, but maybe
it is intentional.
It seems that a diff on a zooreg class that is not _strictly_ regular only
considers those entries that are 'deltat' apart.
In the following, diff on the zooreg class only returns values where
Hi Montserrat,
What exactly would you like to plot?
Your differential equation can be easily integrated so
that you can get an implicit expression for F(x), i.e.
expression like G(c,x,F(x)) = 0 where G is a known
function and c is an arbitrary constant. For every
value of c and each value of x su
Hi,
I would like to plot the following equation:
dF(x)/dx=(k1+k2F(x))(1-F(x))
where k1 and k2 are parameters that I have estimated already.
How can I plot the curve in R?
Thanks!
Montserrat Rue
Universitat de Lleida (Spain)
[[alternative HTML version deleted]]
__
Perhaps more fun is
> library(XML)
> res = htmlTreeParse("http://www.omegahat.org/RSXML/";, useInternalNodes=TRUE)
> xpathApply(res, "//h1", xmlValue)
[[1]]
[1] "An XML package for the S language"
Martin
Quoting Steven McKinney <[EMAIL PROTECTED]>:
>
>
> >-Original Message-
> >From: [
>-Original Message-
>From: [EMAIL PROTECTED] on behalf of Am Stat
>Sent: Wed 8/1/2007 2:19 PM
>To: r-help@stat.math.ethz.ch
>Subject: [R] Extracting a website text content using R
>Dear useR,
>Just wandering whether it is possible that there is any function in R could
>let me get the t
Walter,
Here is what I do in similar situations. I am on WinXP but this should be
similar on other systems (I hope).
1. I start R with the new .RData workspace (usually by double-clicking on
it).
2. I go to File < Change Dir menu item.
3. I change the directory to where the old .RData is.
4.
work with it as text. for text mining use:
1- http://wwwpeople.unil.ch/jean-pierre.mueller/
2- tm by Ingo F.
Am Stat wrote:
> Dear useR,
>
> Just wandering whether it is possible that there is any function in R could
> let me get the text contents for a certain website.
>
> Thanks a lot!
>
> Be
If you are feeling altruistic you could write a predict method for
slm objects, it wouldn't be much work to adapt what is already
available and follow the predict.lm prototype. On the other
hand if you are looking for something quick and dirty you can
always resort to
newX %*% coef(slmo
First thing to do is to use Rprof to determine where the time is being
spent and then you can pinpoint what section of code is taking the
time. A quick look says to do all your subsetting at once. You might
look into using 'split' to create the subsets and then access the
subsets with "[[...]]".
All right, my question is, if there is(are) such function(s), what is(are)
it(they) ?
Best,
Leon
2007/8/1, Bert Gunter <[EMAIL PROTECTED]>:
>
> Yes, there are.
>
> (Please see and follow the posting guide if you wish to obtain something
> more specific)
>
>
> Bert Gunter
> Genetech Nonclinic
On 8/1/07, Jonathan Williams <[EMAIL PROTECTED]> wrote:
> I thought one nice addition to a splom figure would be to have the
> scatterplots in the upper triangle and a color-coordinated
> correlation matrix on the bottom. So I tried my hand at customizing
> panel.pairs(), and was rebuffed. Many t
Hi,
I am trying out the SparseM package and had the a
question. The following piece of code works fine:
...
fit = slm(model, data = trainData, weights = weight)
...
But how do I use the fit object to predict the values
on say a reserved testDataSet? In the regular lm
function I would do so
Hi,
Suppose I have a dataframe in one workspace (a .RData file) dedicated
to one project. I then create a new workspace (a new .RData file) for
another project but I want to move or copy the dataframe to the new
workspace. How can I do this efficiently? Just do a dump and th
Dear R-users,
I have written the following code to generate some trellis plots. It
works perfectly fine except that it is quite slow when it is apply to my
typical datasets (over several thousands of lines). I believe the
problem comes from the loops I am using to subset my data.frame. I read
Yes, there are.
(Please see and follow the posting guide if you wish to obtain something
more specific)
Bert Gunter
Genetech Nonclinical Statistics
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Am Stat
Sent: Wednesday, August 01, 2007 2:19 PM
To: r-h
Dear useR,
Just wandering whether it is possible that there is any function in R could
let me get the text contents for a certain website.
Thanks a lot!
Best,
Leon
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
I'm coming from the scipy community and have been using R on and for
the past week or so. I'm still feeling out the language structure,
but so far so good. I apologize in advance if I pose any obvious
questions, due to my current lack of diction when searching for my
issue, or recognizing it if I
I am not exactly sure what you are asking for. I am assuming that you
want a vector that represent the combinations that are given
combinations that are present:
> N
[1] 11 22 31 42 51 12 21 32 41 52
> table(i,j)
j
i 1 2
1 1 1
2 1 1
3 1 1
4 1 1
5 1 1
> z <- table(i,j)
> which(z==1
If want you want to do is to delete the diagonal values in your
matrix, here is one way of doing it assuming that your matrix is 'x':
# ignore the diagonal value in each column
z <- lapply(1:ncol(x), function(a) x[-a, a])
# put back into a matrix
z <- do.call('cbind', z)
On 8/1/07, Urmi Trivedi
see ?polygon function.
--
Henrique Dallazuanna
Curitiba-Paraná-Brasil
25° 25' 40" S 49° 16' 22" O
On 01/08/07, Nair, Murlidharan T <[EMAIL PROTECTED]> wrote:
>
> Indeed, this is what I wanted, I figured it from the function you and
> Mark pointed me. Thank you both.
>
> I am trying to plot it to
On 01-Aug-07 19:18:05, [EMAIL PROTECTED] wrote:
> Well, is "t = 1.11" all that accurate in the first place? :-)
>
> In fact, reading beween the lines of the original enquiry, what the
> person probably wanted was something like
>
> ta <- pt(-1.11, 9) + pt(1.11, 9, lower.tail = FALSE)
>
> which
On 7/30/07, Patrick Drechsler <[EMAIL PROTECTED]> wrote:
> Hi,
>
> what is the recommended way of adding a strip to a lattice plot?
See ?strip.default.
> In the example below I would like to add the value of mean(y) to a new
> strip.:
>
> --8<---cut here---start---
On 7/30/07, Patrick Drechsler <[EMAIL PROTECTED]> wrote:
> "Deepayan Sarkar" <[EMAIL PROTECTED]> writes:
> > On 7/30/07, Patrick Drechsler <[EMAIL PROTECTED]> wrote:
> >> The Gmane interface seems to have some lag at the moment...
> >>
> >> "Deepayan Sarkar" <[EMAIL PROTECTED]> writes:
> >>
> >> >
Indeed, this is what I wanted, I figured it from the function you and
Mark pointed me. Thank you both.
I am trying to plot it to illustrate the point and I tried this
plot(function(x) dt(x, df = 9), -5, 5, ylim = c(0, 0.5), main="t -
Density", yaxs="i")
Is there an easy way to shade the area un
Hello,
My name is IKJIN LEE and I am studying in University of Iowa.
Recently, I have been trying to prove the fact d F(x, y; rho) / d rho =
f(x, y; rho),
but I could not. In the mean time, I found that Professor S. Le Cessie
from Netherlands kindly gave you an elegant proof of
d F(x, y; rho) /
Well, is "t = 1.11" all that accurate in the first place? :-)
In fact, reading beween the lines of the original enquiry, what the
person probably wanted was something like
ta <- pt(-1.11, 9) + pt(1.11, 9, lower.tail = FALSE)
which is the two-sided t-test tail area.
The teller of the parable wi
csiro.au> writes:
>
> for the upper tail:
>
> > 1-pt(1.11, 9)
> [1] 0.1478873
>
wouldn't
pt(1.11, 9, lower.tail=FALSE)
be more accurate?
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
?pt is what you want.
Hope this is helpful,
Dan
Daniel Nordlund
Bothell, WA
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
> On Behalf Of Nair, Murlidharan T
> Sent: Wednesday, August 01, 2007 11:43 AM
> To: r-help@stat.math.ethz.ch
> Subject: [R] t-distributi
for the upper tail:
> 1-pt(1.11, 9)
[1] 0.1478873
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Nair, Murlidharan
T
Sent: Thursday, 2 August 2007 4:43 AM
To: r-help@stat.math.ethz.ch
Subject: [R] t-distribution
If I have a calculated t can I get the p
if you mean the area to the left of the 1.11 point on the x axis of a t
dist with 9 degrees of freedom,
Then you need to use pt(1.11,9). See ?pt for more info.
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Nair, Murlidharan
T
Sent: Wednesday, August 01,
If I have a calculated t can I get the probability associated with it
using an R function by giving it the df and t? I know I can do the whole
calculation using t.test() or get the t-distribution using qt(). If
t=1.11 and df =9 can I get the probability?
Thanks../Murli
[[alternativ
Hi,
I'm having trouble using the constrOptim function to generate the
9-component vector argmin of the function ELfsds:
ELfsds <- function(pvechat){
LG=0
for(i in 1:9){
LG=LG+log(pvechat[i])
}
return(-LG)
}
with accompanying gradient function:
see ?rect, or, for more general shapes, ?polygon
## EXAMPLES
plot(c(0,500),c(0,500),type="n",las=1)
rect(par("usr")[1],200,par("usr")[2],300,col="grey90")
points(seq(0,500,length=3),seq(0,500,length=3))
plot(c(0,500),c(0,500),type="n",las=1)
polygon((par("usr")[1:2])[c(1,1,2,2)],
(c(200,3
I thought one nice addition to a splom figure would be to have the
scatterplots in the upper triangle and a color-coordinated
correlation matrix on the bottom. So I tried my hand at customizing
panel.pairs(), and was rebuffed. Many times. Four hours of
fruitless debugging later, I turn to yo
08/01/2007 18:46:14 [GMT+0100]
For security reasons certain items found in an email with your address as the
sender have not been accepted.
File name: TRANSCRIPT.SCR
Filtered by: Malformed messages
Sender: r-help@stat.math.ethz.ch
Recipients: [EMAIL PROTECTED]
CC:
It works after rm(list=ls())
I had a function called 'fn'; if I had paid close attention when I
loaded gsubfn, I would have seen the warning. My fault.
Thanks again for a most useful package!
David L. Reiner
Rho Trading Securities, LLC
-Original Message-
From: Gabor Grothendieck [mail
Emilio-
One possible way to do this is to create a new factor and put the levels
in the order you want them in.
For example:
related.differences$header2 <- factor(related.differences$header,
levels=c("spontaneous recovery", "negative reinforcer", ... etc. )
where you want spontaneous recovery
I forget to specify that the data in the table are random. and I want to know
how it is possible to carry out all the two way ANOVA test in a single run for
all the columns. thank you
-- Initial Header ---
>From : [EMAIL PROTECTED]
To :
Cc : "r-help" r-h
Dear R users,
I used the following code to draw a scatter plot.
plot(x,y,type="n")
points(x,y,pch=1)
And then I used the abline functions to draw two lines. I want to add
the shadow between those two lines.
abline(h=200)
abline(h=300)
Any suggestions?
Thanks
Rebecca
-
Dear R-users,
I would like a suggestion.
I have two time series covering the period 1000-2000 a.D. and I would like to
understand if there are some time significant correlations between them.
Samples in the first time series are quite compact because about a sample per
year is available.
I am using R 2.5.0 and lattice 0.15-5. plot1, plot2, and plot3 all
work fine. plot4 gives the error.
x<-1:12
f<-gl(3,4)
g<-gl(4,3)
plot1<-xyplot(y~x|f*g)
plot2<-xyplot(y~x|f*g,scales=list(relation="free"))
plot3<-densityplot(~x|f*g)
plot4<-densityplot(~x|f*g,scales=list(relation="free"))
Than
On 8/1/07, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote:
> Extremely cool and useful
>
> I immediately saw uses for it in some of the work I do.
> So went to look at it. I especially liked the examples with joins!
>
> I observed a few small bugs in the documentation:
>
> - The help references
>
I've got this dataframe
X12 X14 X17 X19 vitigni years
11 4 78 54 rie 2005
21 4 7 4 rie 2005
31 4 75 5 rie 2005
42 5 66 5croa 2005
51 4 4 46croa 2005
62 5 7 6croa 2005
73 2 56 5 rie 2006
83 6
Hello and sorry to bother.
Please help.
I searched the archives but could not find out why --args is being ignored
on Windows 2000.
I try
R CMD BATCH --slave 11.R 11.Rout --args "12"
and 11.R has
x=commandArgs(trail=T)
print(x)
a=x[length(x)]
write.csv(a,file="13.out")
q("no")
the argument i
Dear all,
I have been successful so far in plotting matrices and getting the regression
line. But, as the matrices contains values like 1 and 0 (diagonal line) which
is actually not needed as they are for the same gene, is creating bias in my
regression line. I wish to neglect that part of the
?matrix
On 8/1/07, Dong Guo <[EMAIL PROTECTED]> wrote:
> Thanks, Gabor.
>
> My matrix is from a big array(year, regions,variables). so, matrix does not
> have row names or col names, how could i add the col names or row names??
>
> Thanks again.
> Dong
>
>
> On 8/1/07, Gabor Grothendieck <[EMAIL P
Thanks, Gabor.
My matrix is from a big array(year, regions,variables). so, matrix does not
have row names or col names, how could i add the col names or row names??
Thanks again.
Dong
On 8/1/07, Gabor Grothendieck <[EMAIL PROTECTED]> wrote:
>
> Try this:
>
> mat <- matrix(1:24, 6, dimnames = lis
Try this:
mat <- matrix(1:24, 6, dimnames = list(year = 2001:2006, region = letters[1:4]))
library(lattice)
xyplot(Freq ~ year | region, as.data.frame.table(mat))
On 8/1/07, Dong Guo <[EMAIL PROTECTED]> wrote:
> Dear all,
>
> I have a matrix, dim = (years, regions)
>
> I would like to plot th
Dear all,
I have a matrix, dim = (years, regions)
I would like to plot the data in a lattice so that each panel is region's
plot with y-axis based on values, x-axis based on year.
how can I do that?
Many thanks in advance.
Dong
[[alternative HTML version deleted]]
__
--- David Lloyd <[EMAIL PROTECTED]> wrote:
> Hi,
>
> I've plotted a Kaplan-Meier curve but the curves
> only range from 0.7 to
> 1 on the y-axis. Therefore I have used: -
>
> ylim=c(0.7,1)
>
> [although I think convention dictates that you plot
> 0.5 to 1 to show the
> median? A few papers I'v
Use wtd.mean from Hmisc and by:
Lines <- "
11:00:0134 1000
11:00:0135 500
11:00:0135 1000
11:00:0234 500
11:00:0235 500
"
library(Hmisc) # for wtd.mean
library(zoo)
library(chron)
# replace with DF <- read.table("mytable.dat")
On Wed, 2007-08-01 at 14:06 +0200, Dong GUO 郭东 wrote:
> Dear all,
>
> How can I clear the workspace, as we do in Matlab "clear all"??
>
> Many thanks in advance.
>
> Dong
?rm
E.g.:
rm(list = ls())
will remove everything shown by ls(). Look at ?ls to see possible
arguments to that function t
Dear all,
How can I clear the workspace, as we do in Matlab "clear all"??
Many thanks in advance.
Dong
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read t
Hallo,
Im trying to find out how to tabulate frequencies
of factors when the data have a frequency variable.
e,g:
i<-rep(1:5,2)
j<-rep(1:2,5)
N<-10*i+j
table(i,j) gives a table of ones
as each combination occurs only once.
How does one get a table with the corresponding N's?
Thanks!
Gerrit.
-
Thanks again, Petr. Following the reference, that would be true that "="
only assign values to the top level...So apparently using '<-' is the safe
all the time to assign values.
Dong
On 8/1/07, Petr PIKAL <[EMAIL PROTECTED]> wrote:
>
> Hi
>
> "Dong GUO ¹ù¶«" <[EMAIL PROTECTED]> napsal dne 31.07
David Lloyd wrote:
> Hi,
>
> I've plotted a Kaplan-Meier curve but the curves only range from 0.7 to
> 1 on the y-axis. Therefore I have used: -
>
> ylim=c(0.7,1)
>
> [although I think convention dictates that you plot 0.5 to 1 to show the
> median? A few papers I've read have done this]??
>
>
Hi,
I've plotted a Kaplan-Meier curve but the curves only range from 0.7 to
1 on the y-axis. Therefore I have used: -
ylim=c(0.7,1)
[although I think convention dictates that you plot 0.5 to 1 to show the
median? A few papers I've read have done this]??
BUT, I would like a symbol like // (but r
Hi,
I've plotted a Kaplan-Meier curve but the curves only range from 0.7 to
1 on the y-axis. Therefore I have used: -
ylim=c(0.7,1)
[although I think convention dictates that you plot 0.5 to 1 to show the
median? A few papers I've read have done this]??
BUT, I would like a symbol like // (but r
What about the other way round?
x <- 1:10
y <- rnorm(10,10,1)
x2 <- 3*x + 2
plot(y ~ x2, xaxt = "n")
axis(side=1,at = x2, labels = x)
Florent Bresson
Maybe I do not explain well what I would like to do. I do not want to change
the labels of the axis, but the scale. What I want is a general procedure for
changing the scale. Its like using a logarithmic scale on a plot. Labels are
the same, but the increases of x along the x-axis are defined by
x <- 1:10
y <- rnorm(10,10,1)
x2 <- 3*x + 2
plot(y ~ x, xaxt = "n")
axis(side=1,at = x, labels = x2)
Joris
Florent Bresson
<[EMAIL PROTECTE
Dear R users,
I would like to draw a plot with a custom scale for the axis. More precisely,
instead of plotting y on x, I want to plot y on a monotone function of x (for
instance a*x+b). Which command and/or package should I use in order to get this
result?
Thanks
Florent Bresson
If you put in penalties for breaking constraints,
it is generally better to make the penalty depend
on the size of the violation. This keeps the function
continuous (though usually not differentiable at the
boundary), and gives the optimizer a hint about
which way to go.
Patrick Burns
[EMAIL PRO
Hello,
I'm trying to plot a set of barplots like a matrix (2 rows, 10 columns
from"reduced_mat") to a pdf. It works with the following parameters:
pdf("test.pdf",width=ncol(reduced_mat)*2, height=nrow(reduced_mat)*2, pointsize
= 12)
par(mfcol = c(nrow(reduced_mat),ncol(reduced_mat)), oma = c(0
Hello,
I have two questions concerning the RWeka package:
1.) First question:
How can one perform a cross validation, -say 10fold- for a given data
set and given model ?
2.) Second question
What is the correct syntax for the parametrization of e.g. Kernel
Make a function and add the constraints to the function eg. if (Ai1 > x )
Inf.
By making the result Infinite, you have added the constraints manually.
Vu Nguyen-4 wrote:
>
> Hello R community,
>
> I am sorry for the previous accidental posting as I pressed a wrong key.
>
> I am using R for
> "UweL" == Uwe Ligges <[EMAIL PROTECTED]>
> on Tue, 31 Jul 2007 12:13:45 +0200 writes:
UweL> francogrex wrote:
>> Hi, 2 questions:
[]
>> Question 2:
>>
>>> deriv(~gamma(x),"x")
>>
>> expression({
>> .expr1 <- gamma(x)
>> .value <- .expr1
Hi
"Dong GUO 郭东" <[EMAIL PROTECTED]> napsal dne 31.07.2007 15:27:35:
> Thanks, Petr.
>
> I changed the equation mark from "=" to "<-", then, it works fine. Dont
know
> what difference it has made between the "=" and "<-"..
from help page
The operators <- and = assign into the environment in
On Tue, 2007-07-31 at 15:10 -0700, Jennings, Eric wrote:
> I'm working on a project involving reliability values (known failure
> rates) for a system with approximately 700 components with a set
> cconfiguration.
>
> I'm looking to compute a "parts-count" MTBF (mean time between failures)
> for th
75 matches
Mail list logo