Hi,
Except the function postscript ('file_name.eps'), are there functions to
generate an image file in eps format?
thanks
Carol
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/li
Using tail() for the selection is more elegant and slightly faster.
> N<- 100
> x <- runif(N)
> system.time(x[order(x)[c(N-1,N)]])
user system elapsed
1.080.011.10
> system.time(sort(x)[c(N-1,N)])
user system elapsed
0.360.000.35
> system.time(tail(sort(x), 2)
Hi there.
I am not sure what it is that you are trying to get.
If you want to get the anova results for the interaction, why not just use:
summary(aov(x~gen*con))
?
Are you trying to do something else with the result?
Tal
On Thu, Sep 3, 2009 at 1:32 PM, Luca Servolo wrote:
> When I perfo
On Fri, Sep 04, 2009 at 12:38:29AM -0700, carol white wrote:
> Hi,
> Except the function postscript ('file_name.eps'), are there functions to
> generate an image file in eps format?
Yes, there is cairo_ps and probably there are other ways. But more
importantly: what is the reason for your questio
> David Winsemius
> on Thu, 3 Sep 2009 22:34:05 -0400 writes:
> I'm having difficulties with plot.locfit.3d, at least I think that is
> the problem. I have a large dataframe (about 4 MM cases) and was
> hoping to see a non-parametric estimate of the hazard plotted agai
Hi everyone,
I want to write some data to file, for example
a = c(1,2,3,4,5)
b = c(5,6,7,8,9)
write.table(cbind(a,b), file = "R output.csv", sep = ",")
The result is this:
a b
1 1 5
2 2 6
3 3 7
4 4 8
5 5 9
It added an i
On Thu, 2009-09-03 at 04:49 -0700, swertie wrote:
> Hello, I read a lot about ordination, but I am still confused... I have data
> on species presence/absence for 8 different sites and I would like to
> represent my species and the sites on an ordination plot to see if some
> species are associated
Hi,
This may come close to what you want,
x <- data.frame(ID=rep(letters[1:5],2), A1=rep(10:14,2), A2=rep(2:6,2),
A3=c(101:105,95:99), A4=c(-60, rep(c(0, 3), length=9)))
# basic conditions
cond1 <- quote(ID == "a" & A2 >1)
cond2 <- quote(A1 < 10)
cond3 <- quote(A1 == 10) # note t
Hi everybody,
I am struggling to save txt. or csv. large sparse matrices (1024 X 768
image size), initialized at 0.1.
image1[0:10,0:10]
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,] 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1
[2,] 0.1 0.1 0.1 0.1 0.1 0.1 0.1
On Fri, 4 Sep 2009 01:07:49 -0700 (PDT) dennis11
wrote:
D> It added an index and shifted the column names by one, b is above a,
D> and a is above the added index.
Try the option row.names=FALSE in write.table then no index is written.
Stefan
__
R-hel
thomas.david...@sauder.ubc.ca wrote:
I did a speed test with a colleague. We basically have identical Lenovo
ThinkCentres. He has 8 vs my 4 G RAM, but I don't think that's the
issue.
code:
length <- 2*10^6
a <- runif(length)
b <- runif(length)
print(summary(lm(a~b)))
Running Arch x86_64 thi
Hello Rasmus,
> I am not sure whether this is a bug or lack of R experience.
It is not your lack of R experience, and whether it is a bug
I am not sure ;-)
> However, I am using your Sweavelistingutil package, which is very
> nice.
Thank you.
> Obviously I use it to create LaTeX files. These
AFAIK, R only has "pnorm" which computes the probability of getting a
value smaller or equal to "x" from
a normal distribution N[mean,stdev]. For example:
R> pnorm(0, 4, 10)
[1] 0.3446
means there is 34.46% chance of getting a value equal to or smaller
than 0 from a N(4, 10) distribution.
What
Moreno Ignazio Coco wrote:
Hi everybody,
I am struggling to save txt. or csv. large sparse matrices (1024 X 768
image size), initialized at 0.1.
image1[0:10,0:10]
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,] 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1
[2,] 0.1 0.1
Hi,
Does anybody know, what is going on here?
diag(sqrt(1))
[,1]
[1,]1
diag(sqrt(0.))
<0 x 0 matrix>
sqrt(1)
[1] 1
sqrt(0.)
[1] 0.5773214
BR, Markku Karhunen
researcher
University of Helsinki
__
R-help@r-project.org mailing l
Please tell us what you did exactly when you got the empty file:
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
well, that's exactly the point, I didn't do anything fancy or special,
just save th
On 09/04/2009 12:25 PM, Markku Karhunen wrote:
Hi,
Does anybody know, what is going on here?
diag(sqrt(1))
[,1]
[1,] 1
diag(sqrt(0.))
<0 x 0 matrix>
sqrt(1)
[1] 1
sqrt(0.)
[1] 0.5773214
BR, Markku Karhunen
researcher
University of Helsinki
Try this instead;
> diag( sqrt(.33
Markku Karhunen wrote:
Hi,
Does anybody know, what is going on here?
diag( x ) produces a round(x) x round(x) identity matrix when x is
length 1. (This is the third case listed on the man page ?diag). See
the note there about a safer form if you wanted a matrix with x on the
diagonal.
On Fri, Sep 4, 2009 at 11:25 AM, Markku
Karhunen wrote:
> Hi,
>
> Does anybody know, what is going on here?
>
>> diag(sqrt(1))
>
> [,1]
> [1,] 1
>>
>> diag(sqrt(0.))
>
> <0 x 0 matrix>
>>
>> sqrt(1)
>
> [1] 1
>>
>> sqrt(0.)
>
> [1] 0.5773214
>
Read the help for diag yet?
'diag
it's documented as "unexpected"
?diag
Note
Using diag(x) can have unexpected effects if x is a vector that could be of
length one. Use diag(x, nrow = length(x)) for consistent behaviour.
And the result follows from this part,
else if (length(x) == 1L && nargs() == 1L) {
n <- as.int
dear Claudia,
I was recently in touch with Vito Muggeo (the developer of the
segmented package) with a similar question. This is an adapted
version of his answer to your problem.
In fact, the essential aspect is that predict.segmented has not
(yet?) been implemented. Nevertheless, you could u
Hi
it is rather difficult to understand what you mean by your
questions/answers without real reproducible code.
r-help-boun...@r-project.org napsal dne 03.09.2009 13:41:11:
>
> I'm posting answers to my own Q's here - as far as I have answers -
first so
> that people don't spend time on them,
Dear R-help list,
Does anyone have a function that I could use to determine power for 2 way
Anova??
an A x B repeated measures study,power is 0.95, I'd like to draw separate lines
for three different combinations of A and B:
(2,2), (2,5), (2,8).
Thanks a lot.
Tammy
_
baptiste auguie wrote:
it's documented as "unexpected"
?diag
Note
Using diag(x) can have unexpected effects if x is a vector that could be of
length one. Use diag(x, nrow = length(x)) for consistent behaviour.
And the result follows from this part,
else if (length(x) == 1L && nargs() ==
Duncan Murdoch wrote:
baptiste auguie wrote:
it's documented as "unexpected"
?diag
Note
Using diag(x) can have unexpected effects if x is a vector that could be of
length one. Use diag(x, nrow = length(x)) for consistent behaviour.
And the result follows from this part,
else if (leng
Hi,
I'm trying to fit a smooth line in a plot(y ~ x) graph.
x is continuous variable
y is a proportion of success in sub-samples, 0 <= y <= 1, from a Monte
Carlo simulation.
For each x there may be several y-values from different runs. Each run
produces several sub-samples, where "0" mean no suc
Hi,
Is there a way I can specify linear contrasts in glm? I'm looking for
something equivalent to SAS' contrast statement.
I'd like to do the following, suppose I have a categorical input with 4
levels (a,b,c,d), I'd like to test something like: (i) a+b=c+d, (ii) a=b,
(iii) a=b+d, etc...
Thanks
You could use a glm with the binomial family to model that.
A solution with ggplot2
library(ggplot2)
ggplot(dataset, aes(x = x, y = y, weights = n)) +
geom_smooth(method = "glm", family = binomial)
geom_point()
Have a look at glht() from the multcomp package. That will allow you to
specify the contrasts for posthoc tests.
A more direct approach is the create the required dummy variables
yourself.
HTH,
Thierry
ir. Thierry Onk
Hi
image<-matrix(0.1, 768,1024)
> image[1:10, 1:10]
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,] 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1
[2,] 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1
[3,] 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1
[4,] 0.1 0.1 0
Hi,
I've been trying to pass a character vector from R to a FORTRAN subroutine.
There have been several posts discussing this issue (e.g.
http://tolstoy.newcastle.edu.au/R/help/98a/0547.html,
http://tolstoy.newcastle.edu.au/R/help/05/10/13558.html,
http://tolstoy.newcastle.edu.au/R/help/01a/25
True. Should have read ?diag.
However, this provokes a more general question: Is there some way I
can declare some scalar and _all its functions_ as matrices?
For instance, I would like to
A = as.matrix(0.98)
B = function(A)
C = diag(sqrt(B))
so that all scalars are explicitly [1,1] matric
Hi,
I'm trying to fit a smooth line in a plot(y ~ x) graph.
x is continuous variable
y is a proportion of success in sub-samples, 0 <= y <= 1, from a Monte
Carlo simulation.
For each x there may be several y-values from different runs. Each run
produces several sub-samples, where "0" mean no suc
Dear all,
I am new to R and would like to run a multinomial logistic regression on my
dataset (3 predictors for 1 dependent variables)
I have used the vglm function from the VGAM package and got some results. Using
the predict() function, I obtained the probability table I was looking for.
Ho
Nordlund, Dan (DSHS/RDA) wrote:
>
>
> You need to put the filename in quotes
>
> file = "C:/Documents and Settings/aelmore/Desktop/foo.pdf"
>
> Hope this is helpful,
>
> Dan
Oh yes, this was very helpful (and don't I feel just too silly...). I am,
however, still stuck. This is like pullin
Thank you very much. I am just concerned because I wonder if I used the best
method. I have presence/absence data. With isoMDS I can specify
"Bray-Curtis" distance method, which is adequate, but I was not sure of the
method used by metaMDS. I think that it is Euclidian distance. Can I use it
for p
Gundala Viswanath wrote:
>
> AFAIK, R only has "pnorm" which computes the probability of getting a
> value smaller or equal to "x" from
> a normal distribution N[mean,stdev]. For example:
>
> R> pnorm(0, 4, 10)
> [1] 0.3446
>
> means there is 34.46% chance of getting a value equal to or smal
Hello ALL!
Can anyone tell me how to specify qqmath() function with distribution
qgamma?
In help it is bit vague how to pass shape and scale parameters for gamma.
Thank you in advance. Best,
PM
__
R-help@r-project.org mailing list
https://stat.ethz.c
On 9/4/2009 8:36 AM, Petar Milin wrote:
Hello ALL!
Can anyone tell me how to specify qqmath() function with distribution
qgamma?
In help it is bit vague how to pass shape and scale parameters for gamma.
Take a look at the first example, which passes in the df parameter for a
t distribution.
Hi
r-help-boun...@r-project.org napsal dne 04.09.2009 13:53:24:
>
> Nordlund, Dan (DSHS/RDA) wrote:
> >
> >
> > You need to put the filename in quotes
> >
> > file = "C:/Documents and Settings/aelmore/Desktop/foo.pdf"
> >
> > Hope this is helpful,
> >
> > Dan
>
>
> Oh yes, this was very h
Hello !
I want to create a spatial stratified sampling scheme with the package
spsurvey. To do this with the function "grts" in spsurvey, I need to create a
list containing the specifications for each stratum. This specifications were
stored in a named list, where the name for each stratum is t
Dear R-ers!
I have a badly organized data base in Excel. Once I read it into R it
looks like this (all variables become factors because of many spaces
and other characters in Excel):
x<-data.frame(A=c("","Name1","text1","text2","text3","","","","Name2","text1","text2","text3","","","Name3","text1
I recently submitted a package on CRAN. When I went to update the package, I
tried to do the check/build procedures on a system with Windows XP, using the
recommended R-tools available from Dr. Ripley (Rtools29).
Everything works fine in the check and build process except under Windows XP
the p
On 9/4/2009 9:07 AM, Gates, Michael wrote:
I recently submitted a package on CRAN. When I went to update the package, I
tried to do the check/build procedures on a system with Windows XP, using the
recommended R-tools available from Dr. Ripley (Rtools29).
Everything works fine in the check and
Try this:
lapply(grep("Name", x$A), function(idx)x[idx + 1:3,])
On Fri, Sep 4, 2009 at 9:57 AM, Dimitri Liakhovitski wrote:
> Dear R-ers!
>
> I have a badly organized data base in Excel. Once I read it into R it
> looks like this (all variables become factors because of many spaces
> and other c
Try this:
x2 <- subset(x, A != "")
split(x2, cumsum(x2$B == ""))
On Fri, Sep 4, 2009 at 8:57 AM, Dimitri Liakhovitski wrote:
> Dear R-ers!
>
> I have a badly organized data base in Excel. Once I read it into R it
> looks like this (all variables become factors because of many spaces
> and other
"Quote from Petr PIKAL"
> No, no. Here you need to do your plotting routine.
> Something like
> plot(1:10, rnorm(10))
> Now you can open it by e.g. Acrobat Reader or any suitable programme. Do
> not forget to quit from Acrobat when trying to use the file again, if it
> is opened it ca
On 04-Sep-09 10:45:27, Markku Karhunen wrote:
> True. Should have read ?diag.
>
> However, this provokes a more general question: Is there some way I
> can declare some scalar and _all its functions_ as matrices?
>
> For instance, I would like to
>
> A = as.matrix(0.98)
> B = function(A)
> C =
Dear list,
I am trying to make eps file with embedded font.
I use:
postscript("ranking-exp-all.eps", horizontal=TRUE, onefile=FALSE,
paper="special", height=8, width=12, family="Helvetica")
# plot stuff
dev.off()
since R does not embed font, I then use:
embedFonts(file="indegdistr.eps", outfil
On Fri, 2009-09-04 at 03:36 -0700, swertie wrote:
> Thank you very much. I am just concerned because I wonder if I used the best
> method. I have presence/absence data. With isoMDS I can specify
> "Bray-Curtis" distance method, which is adequate, but I was not sure of the
> method used by metaMDS.
On 04-Sep-09 14:01:44, Simone Gabbriellini wrote:
> Dear list,
> I am trying to make eps file with embedded font.
> I use:
>
> postscript("ranking-exp-all.eps", horizontal=TRUE, onefile=FALSE,
> paper="special", height=8, width=12, family="Helvetica")
># plot stuff
> dev.off()
>
> since R does no
Rainer M Krug wrote:
Hi
is there a list of all roxygen tags which are available? I couldn't find them.
I am asking specifically towards the use of roxygen in documenting S4
classes - is that implemented yet (i am using roxygen 0.1 from CRAN at
the moment)?
Thanks
Rainer
I am using it do
Dear all,
I would like to calculate a cca (package vegan) with species and environmental
data. One of these environmental variables is cos(EXPOSURE).
The problem: for flat releves there is no exposure. The value is missing and I
can't call it 0 as 0 stands for east and west.
The cca does not run
thanks Jonathan,
I was wondering about the difference between your second option and
the Ted one: isn't it the same thing?
regards,
Simone
2009/9/4 Jonathan Baron :
> A couple of other ideas about embedding fonts and setting bounding
> boxes. These all work on Linux, so in theory they should al
Hi
Well, I think the title says it all! I've looked through the documentation but
I can't find a way of doing this. The situation is that I have 4 series, say
a, b, c and d. Series a and c are plotted on the lh y axis, series b and d are
plotted on the rh (secondary) y axis. I've worked out
Hi all,
I have got 2 function (see bellow) which are simplifications of what I need
to do. These functions are precisely the same, except for the last line.
My question is, why doesn't function testA work in the same way as function
testB.
Both functions produce two objects, "a" and "b" that
Ah! sorted!
it was NOT running the same code.
We're making a GUI using Perl (Tcl/Tkx) to facilitate a number of
analyses in our lab to other people who don't necessarily want to know
about R (their loss ;-)
I provided the R code to my colleague and he assured me he used it
without chang
Thank you very much for those useful informations. I've been reading some
papers and actually different people will use different ordination methods
also if the studies are very alike. So I will keep metaMDS for the moment
and see if my results are interpretable :)
--
View this message in contex
Hello,
I have I problem and I hope someone can help me.
I have a data set with 66 variables and 30 values for each variable. I need
to calculate the Pearson correlation and the p-values among variables. I
have the “large p, small n” problem so I used a shrinkage process to obtain
the r values. M
Hello,
I have I problem and I hope someone can help me.
I have a data set with 66 variables and 30 values for each variable. I need
to calculate the Pearson correlation and the p-values among variables. I
have the “large p, small n” problem so I used a shrinkage process to obtain
the r values. M
Thank for the help,
unfortunately it doesn't work... it looks impossible to embed fonts...
cairo should do this automatically but not for eps files at the
moment... really don't know how to solve it.
best regards,
Simone
2009/9/4 Ted Harding :
> On 04-Sep-09 14:01:44, Simone Gabbriellini wrote:
-- Forwarded message --
From: Jonathan Baron
Date: 2009/9/4
Subject: Re: [R] eps file with embedded font
To: Simone Gabbriellini
On 09/04/09 17:16, Simone Gabbriellini wrote:
> thanks Jonathan,
>
> I was wondering about the difference between your second option and
> the Ted one
Here is an example:
http://casoilresource.lawr.ucdavis.edu/drupal/node/510
make.groups() is your friend.
Cheers,
Dylan
On Fri, Sep 4, 2009 at 8:28 AM, Paul Sweeting wrote:
> Hi
>
> Well, I think the title says it all! I've looked through the documentation
> but I can't find a way of doing thi
Thank you very much for the answer!
However, instead of x being some uniform values, I have real values.
Hence, if I have dat$V1 of my interest, formula should be like this:
qqmath(~ V1, data=dat, distribution=function(V1) qgamma(V1,
shape=gammafit(dat$V1)$shape, scale=gammafit(dat$V1)$scale))
I
Hi,
I think you've got a problem with environments,
testA<-function(input=1)
{
dat <- data.frame(A=seq(input,5), B=seq(6,10))
vec.names<- c("a", "b")
env <- new.env()
for(i in 1:ncol(dat))
{
tab<- dat[,i]-1
assign(vec.names[i], tab, env=env)
}
I need to use a red-black tree, which package provides that data structure?
--
Best regards
Rune Schjellerup Philosof
Ph.d-stipendiat, Forskningsenheden for Biostatistik
Telefon: 6550 3607
E-mail: rphilo...@health.sdu.dk
Adresse: J.B. Winsløwsvej 9, 5000 Odense C
SYDDANSK UNIVERSITET
If you represent your series as zoo series then you can do this:
library(zoo)
z <- zoo(cbind(a = 1:3, b = 4:6, c = 3:1, d = 6:4))
xyplot(z, screen = 1, type = c("p", "p", "l", "l"))
See the three vignettes that come with zoo and also see ?xyplot.zoo
If you omit screen = 1 then the 4 will be on se
You should pay close attention to the definitions of the parameters
passed to the various gamma functions.
x <- rgamma(1000, 2, 2)
library(mhsmm)
qsh <- gammafit(x)$shape; qsc <- gammafit(x)$scale
> qsh
x
x 2.009260
> qsc
x
x 0.4827448
On Sep 4, 2009, at 11:59 AM, Petar M
To all (and Duncan in particular),
The problem has been solved. Apparently, you cannot have directories in the
path for files that contain "spaces" in their names. I noticed this problem
under Vista if the build directory was nested under several directories within
the working directory (i.e. th
On 9/4/2009 1:02 PM, Gates, Michael wrote:
To all (and Duncan in particular),
The problem has been solved. Apparently, you cannot have directories in the path for
files that contain "spaces" in their names. I noticed this problem under Vista
if the build directory was nested under several direc
This must be explained somewhere, but I've been searching for a couple of
hours and not found it.
What happened to ggplot? It appears to be missing on CRAN, except in the
archives.
http://cran.r-project.org/web/packages/ggplot/index.html
Has ggplot2 replaced ggplot?
I was trying to run some exam
Hi R people,
I have a very basic question to ask - I'm sorry if it's been asked before, but
I searched the archives and could not find an answer. All the examples I found
were much more complicated/nuanced versions of the problem - my question is
much more simple.
I have data with multiple, n
On Fri, 2009-09-04 at 17:15 +0200, Kim Vanselow wrote:
> Dear all,
> I would like to calculate a cca (package vegan) with species and
> environmental data. One of these environmental variables is
> cos(EXPOSURE).
> The problem: for flat releves there is no exposure. The value is
> missing and I can
A couple of other ideas about embedding fonts and setting bounding
boxes. These all work on Linux, so in theory they should also work on
OS X, although I have no idea how.
1. For setting bounding boxes, you can use gv, which is a PostScript
viewer. As you move the pointer around, you can see the
Hello,
I am using R to analyze a large multilevel data set, using
lmer() to model my data, and using anova() to compare the fit of various
models. When I run two models, the output of each model is generated
correctly as far as I can tell (e.g. summary(f1) and summary(f2) for the
multilevel mode
Hi Everyone,
I'm trying to write selfStart non-linear models for use with nls. In these
models some combinations of parameter values are illegal; the function value
is undefined.
That's OK when calling the function directly [e.g. SSmodel(x, pars...)]; I
return an appropriate non-value such as
I want to print a variable with in double quotes.
For example
x = 10 ;
x ;#prints 10
"x" ; #prints x
\"x\" ; # Error: unexpected input in "\"
I want to the out put as'10' or "10"
Thank you
Sailu Yellaboina
__
R-help@r-project.org mailing l
Hi,
Try this,
x = 10
noquote(dQuote(x))
noquote(sQuote(x))
HTH,
baptiste
2009/9/4 sailu Yellaboina
> I want to print a variable with in double quotes.
> For example
>
> x = 10 ;
> x ;#prints 10
> "x" ; #prints x
> \"x\" ; # Error: unexpected input in "\"
>
> I want to the out put
On 9/4/2009 12:18 PM, sailu Yellaboina wrote:
I want to print a variable with in double quotes.
For example
x = 10 ;
x ;#prints 10
"x" ; #prints x
\"x\" ; # Error: unexpected input in "\"
I want to the out put as'10' or "10"
There are lots of ways to do that. The simplest is to c
On Fri, Sep 4, 2009 at 8:28 AM, Paul Sweeting wrote:
> Hi
>
> Well, I think the title says it all! I've looked through the documentation
> but I can't find a way of doing this. The situation is that I have 4 series,
> say a, b, c and d. Series a and c are plotted on the lh y axis, series b and
do you mean something like this:
x <- 10
as.character(x)
I hope it helps.
Best,
Dimitris
sailu Yellaboina wrote:
I want to print a variable with in double quotes.
For example
x = 10 ;
x ;#prints 10
"x" ; #prints x
\"x\" ; # Error: unexpected input in "\"
I want to the out put as
Dear All,
I would like to do run R without having write permissions to any
directory on the system. It seems that I need to modify the R source
code for this, to make R start without creating a temporary directory.
So far, so good. But should I expect any more complications? Does R
really need th
Hi,
A error message arose while I was trying to fit a ordinal model with lrm() I am
using R 2.8 with Design package.
Here is a small set of mydata:
RC RS Sex CovACovBCovCCovDCovE
2 1 0 1 1 0 -0.0055752802
2 1 0
My guess would be:
"Likelihood comparisons are not meaningful for objects fit using restricted
maximum likelihood and with different fixed effects. "
(from ?anova.lme in the nlme package).
Are you using the REML = TRUE default?
Bert Gunter
Genentech Nonclinical Statistics
-Original Message
In R and experimental or mixed-model terminology, your lm model specifies
fixed effects. As long as each data row represents a unique subject, you are
fine with lm. If not, you have to account for the repeated measurement of
subjects and will need other methods (potentially involving random effects
Hi Gene,
Yes, ggplot2 has replaced ggplot.
Hadley
On Fri, Sep 4, 2009 at 12:58 PM, Gene Leynes wrote:
> This must be explained somewhere, but I've been searching for a couple of
> hours and not found it.
>
> What happened to ggplot? It appears to be missing on CRAN, except in the
> archives.
>
The idea of estimating 7 parameters from 11 cases using mostly binary
variables just seems begging for a singular matrix.
On Sep 4, 2009, at 3:59 PM, Ping-Hsun Hsieh wrote:
Hi,
A error message arose while I was trying to fit a ordinal model with
lrm() I am using R 2.8 with Design package.
Ping-Hsun Hsieh wrote:
Hi,
A error message arose while I was trying to fit a ordinal model with lrm() I am
using R 2.8 with Design package.
Here is a small set of mydata:
RC RS Sex CovACovBCovCCovDCovE
2 1 0 1 1 0 -0.005575280
Hi,
I don't see what the difference between read.delim and read.delim2
after reading the help. Can somebody let me know what it is?
Regards,
Peng
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posti
Dan,
> May I ask you to report this as an issue on the oce webpage,
> so that others can see the discussion? (The "R help" is
> perhaps not the right place to report bugs ... and, yes, this
> is a bug.)
> http://code.google.com/p/r-oce/issues/list
Yes, I will.
> A possible solu
Hello,
In limmaGUI, after running a linear model for a single gene & time point
(mutant vs control), the venn-diagram for differentially expressed genes
with a p-value cut-off gives a different number than the top-table with
the same cut-off (2600 vs. 791).
I have seen the same thing happen befor
It's pretty subtle, I will admit, but look more carefully at the dec=
parameter in the help page. Kind of like the historical battles
between between the Liebnizians and the Newtonians.
On Sep 4, 2009, at 5:58 PM, Peng Yu wrote:
Hi,
I don't see what the difference between read.delim and re
Thanks for posting your code, Scott. I'll have a look at the windrose
in OCE sometime next week, to see if there is anything I can do that
would have made your task easier.
I've never used windroses in my own work, so that may explain why the
interface is clunky! Also, it's not completely
d.le...@warwick.ac.uk wrote:
> Hello,
>
> In limmaGUI, after running a linear model for a single gene & time point
> (mutant vs control), the venn-diagram for differentially expressed genes
> with a p-value cut-off gives a different number than the top-table with
> the same cut-off (2600 vs. 791).
Hi,
Suppose I have the following two matrices.
> x=matrix(1:4,nr=2)
> y=matrix(5:8,nr=2)
I want to join the two matrices to get the following two new matrices.
I did a search but I don't find any relevant webpage. Can somebody let
me know how to do it?
[,1] [,2]
[1,]13
[2,]2
On Sep 4, 2009, at 6:26 PM, Peng Yu wrote:
Hi,
Suppose I have the following two matrices.
x=matrix(1:4,nr=2)
y=matrix(5:8,nr=2)
I want to join the two matrices to get the following two new matrices.
I did a search but I don't find any relevant webpage. Can somebody let
me know how to do it?
Hi Peng,
Try
rbind(x,y)
cbind(x,y)
HTH,
Jorge
On Fri, Sep 4, 2009 at 6:26 PM, Peng Yu wrote:
> Hi,
>
> Suppose I have the following two matrices.
> > x=matrix(1:4,nr=2)
> > y=matrix(5:8,nr=2)
>
> I want to join the two matrices to get the following two new matrices.
> I did a search but I don'
On 04/09/2009 3:42 PM, Gábor Csárdi wrote:
Dear All,
I would like to do run R without having write permissions to any
directory on the system. It seems that I need to modify the R source
code for this, to make R start without creating a temporary directory.
So far, so good. But should I expect
If you are willing to do a bit of work you could look into linking R
to an embedded file system in RAM -- maybe whefs.
On Fri, Sep 4, 2009 at 3:42 PM, Gábor Csárdi wrote:
> Dear All,
>
> I would like to do run R without having write permissions to any
> directory on the system. It seems that I nee
> -Original Message-
> From: r-help-boun...@r-project.org
> [mailto:r-help-boun...@r-project.org] On Behalf Of Duncan Murdoch
> Sent: Friday, September 04, 2009 5:05 PM
> To: Gábor Csárdi
> Cc: R mailing list
> Subject: Re: [R] Running R on read-only file system,without
> temporary direct
1 - 100 of 119 matches
Mail list logo