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
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 lucaserv...@gmail.com
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 question?
David Winsemius dwinsem...@comcast.net
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
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
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 the
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
detebe...@hotmail.com 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
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 this
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 are
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 I
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
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
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
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(
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
On Fri, Sep 4, 2009 at 11:25 AM, Markku
Karhunenmarkku.karhu...@helsinki.fi 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 -
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
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() == 1L)
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
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
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...
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
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
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,
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]
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
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.
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 pulling teeth,
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
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 smaller
than 0
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
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 helpful (and don't I feel
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
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):
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
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
Try this:
lapply(grep(Name, x$A), function(idx)x[idx + 1:3,])
On Fri, Sep 4, 2009 at 9:57 AM, Dimitri Liakhovitski ld7...@gmail.comwrote:
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
Try this:
x2 - subset(x, A != )
split(x2, cumsum(x2$B == ))
On Fri, Sep 4, 2009 at 8:57 AM, Dimitri Liakhovitskild7...@gmail.com 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
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 can not
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,
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. I
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 ba...@psych.upenn.edu:
A couple of other ideas about embedding fonts and setting bounding
boxes. These all work on Linux, so in
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
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
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
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.
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.
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 ted.hard...@manchester.ac.uk:
On 04-Sep-09 14:01:44,
-- Forwarded message --
From: Jonathan Baron ba...@psych.upenn.edu
Date: 2009/9/4
Subject: Re: [R] eps file with embedded font
To: Simone Gabbriellini simone.gabbriell...@gmail.com
On 09/04/09 17:16, Simone Gabbriellini wrote:
thanks Jonathan,
I was wondering about the
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 Sweetingm...@paulsweeting.co.uk wrote:
Hi
Well, I think the title says it all! I've looked through the documentation
but I can't
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))
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 separate
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 Milin
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. the
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
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
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,
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't
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
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
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
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 list
Hi,
Try this,
x = 10
noquote(dQuote(x))
noquote(sQuote(x))
HTH,
baptiste
2009/9/4 sailu Yellaboina bio.sa...@gmail.com
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
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 convert
On Fri, Sep 4, 2009 at 8:28 AM, Paul Sweetingm...@paulsweeting.co.uk 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
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'10'
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
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
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
Hi Gene,
Yes, ggplot2 has replaced ggplot.
Hadley
On Fri, Sep 4, 2009 at 12:58 PM, Gene Leynesgleyne...@gmail.com 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 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
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 solution is
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
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
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
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).
I
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,]24
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
Hi Peng,
Try
rbind(x,y)
cbind(x,y)
HTH,
Jorge
On Fri, Sep 4, 2009 at 6:26 PM, Peng Yu pengyu...@gmail.com 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
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árdicsa...@rmki.kfki.hu wrote:
Dear All,
I would like to do run R without having write permissions to any
directory on the system. It
-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 directory
Hi Bert,
Thank you for your note! I tried changing the REML default, and it still
produces the same result (see below). Is that what you meant for me to try?
Incidentally, I am using lmer() not lme()
### ORIGINAL ###
f1 - (lmer(outcome ~ predictor.1 + (1 | person), data=i))
f2 -
1 - 100 of 113 matches
Mail list logo
