On Wed, Dec 01, 2010 at 07:12:37PM -0500, David Winsemius wrote:
...
> One of the workarounds... the one I remember anyway... has been to
> insert:
>
> if(FALSE) {
> stuff your don't want executed
> }
>
Switching a block of code off/on with editing a single character
may be done usin
Hi R users,
Thanks in advance.
I am using R 2.12.0 on Windows XP.
May I request you to assist me in the following please.
1. I am getting error while downloading quote data from yahoo finance.
The example code is below (taken from tseries help):
library(tseries)
con <- url("http://quo
Hi Ivan,
Thanks for your further advice.
RODBC works without problem. Following steps are correct.
> library(RODBC)
> data=odbcConnectExcel(file.choose())
> sqlTables(data)
TABLE_CAT TABLE_SCHEM
1 C:\\Users\\satimiswin764\\Documents\\copy_of_RDemoDe
Look at the code for plot.default to see how it handles this issue
usling locally defined versions of functions with (possibly) locally
set parameter values.
Cheers,
Bert
On Wed, Dec 1, 2010 at 6:39 PM, Gabor Grothendieck
wrote:
> On Wed, Dec 1, 2010 at 2:34 PM, wrote:
>> Dear R-users,
>>
>>
On Wed, Dec 1, 2010 at 2:34 PM, wrote:
> Dear R-users,
>
> I'm trying to work out a way to set default values for arguments in a
> function which includes the optional argument '...'.
> In practice, I have a 'plot' method for a function which specifies different
> types of plots. Every differen
On 2010-12-01 15:59, Jonathan Flowers wrote:
Thanks Sarah, you are right. In perl, "cut" serves the function of
eliminating a block of code from being executed in a script. When "=cut" is
placed above and below the code that you do not wish to execute then the
interpreter will skip over the co
On Dec 1, 2010, at 6:59 PM, Jonathan Flowers wrote:
Thanks Sarah, you are right. In perl, "cut" serves the function of
eliminating a block of code from being executed in a script. When
"=cut" is
placed above and below the code that you do not wish to execute then
the
interpreter will ski
On 2010-12-01 07:09, Robert Quinn wrote:
I am having problems trying to get R to graph data input that is log-normal
on the horizontal (x) axis.
The data is log (base 10), and I am more interested in viewing the tails of
the distribution. The closest I can get with this is log on the vertical
(
Hi all,
How can a multivariate Poisson time series be modeled? Aspects of glm,
forecast, dse and dynlm seem relevant but not quite complete--but hopefully
what I am missing is how to assemble them effectively. What I am looking to
do is model my dependent variable y_t as a Poisson family functi
Thanks Sarah, you are right. In perl, "cut" serves the function of
eliminating a block of code from being executed in a script. When "=cut" is
placed above and below the code that you do not wish to execute then the
interpreter will skip over the code. There are lots of ways to solve the
proble
It would help if you told us what you wanted this function to do,
and provided an example. Not everyone speaks Perl.
Sarah
On Wed, Dec 1, 2010 at 6:10 PM, Jonathan Flowers
wrote:
> Does anyone know of a command in R that is equivalent to the =cut function
> in Perl?
>
> Thanks.
>
> Jonathan
>
-
Does anyone know of a command in R that is equivalent to the =cut function
in Perl?
Thanks.
Jonathan
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the post
Judit -
Does this example provide any insight?
pop <- matrix(0,7,7)# Population matrix when t==0
pop[sample(1:length(pop),1)] <-1 #Population matrix when t==1
which(pop!=0,arr.ind=TRUE)
row col
[1,] 3 4
- Phil Spector
Dear all of you,
I would like if someone can help me to know how I can read a matrix looking for
positions in the matrix that are not null.
I have a matrix with most of the "cells" or "positions" null, I am trying to
know the invasion power of a plant (number 1 in the matrix when t==1).
I would
Is there a way to run TukeyHSD with Type III instead of Type I SS?
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-proje
Phoebe,
In addition to the barplot method below, you really can use plot() to
draw what it sounds like you are looking for IF you have a categorical
variable. To illustrate, try running the following code:
x=sample(c("Richard", "Minnie", "Albert", "Helen", "Joe", "Kingston"),
50, replace
Hi
On 1/12/2010 4:44 a.m., Peter Ehlers wrote:
On 2010-11-30 04:56, Ben Tupper wrote:
Hi,
I thought it might help if I posted the resulting images.
This is the pdf file where the map polygons are not clipped to the
plotting boundary.
http://dl.dropbox.com/u/8433654/test-map.pdf
And this is
thank you, I'll have a good look and come back to you if necessary
Dr. Iasonas Lamprianou
Assistant Professor (Educational Research and Evaluation)
Department of Education Sciences
European University-Cyprus
P.O. Box 22006
1516 Nicosia
Cyprus
Tel.: +357-22-713178
Fax: +357-22-590539
Honor
Hi Phoebe,
Try
x <- c(12, 33, 56, 67, 15, 66)
names(x) <- c('Richard','Minnie','Albert','Helen','Joe','Kingston')
barplot(x, las = 1, space = 0)
HTH,
Jorge
On Wed, Dec 1, 2010 at 2:51 PM, phoebe kong <> wrote:
> Hi,
>
> Can someone tell me how to draw a histogram for the following summary?
>
Antonio,
You need to compare the names of list(...) with the arguments you wish to
check. Here's one way to do so (note that I replaced c with d, because c
is a function):
f <- function(a, ...) {
argnames <- names(list(...))
# check whether b is an argument
if(!("b" %in% a
This would not draw a histogram.
On Wed, Dec 1, 2010 at 11:54 AM, Andrew Miles wrote:
> Try:
>
> plot (myCatVariable)
>
> Andrew Miles
> Department of Sociology
> Duke University
>
>
> On Dec 1, 2010, at 2:51 PM, phoebe kong wrote:
>
> Hi,
>>
>> Can someone tell me how to draw a histogram for th
Hi:
When the question is: 'Are there any functions in R for ?',
package sos is a good place to start:
library(sos) # install if necessary
findFn('MMeans')
Well actually, MMMeans and MMruns came up dry, but Michaelis-Menten is a
commonly occurring model in the biosciences:
findFn('Michaelis-
That did it. Thanks!
On Dec 1, 2010, at 12:38 PM, Paul Murrell wrote:
> Hi
>
> On 2/12/2010 4:55 a.m., Mark Ebbert wrote:
>> Dear R Gurus,
>>
>> I have a fairly simple problem, but I haven't been able to find the
>> answer on 'the google' or in the r-help archives.
>>
>> I am generating plots
Dear All,
I am using a Procrustes analysis to compare two NMDS ordinations for the
same set of sites. One ordination is based on fish data, the other is based
on invertebrate data. Ordinations were derived using metaMDS() from the
{vegan} library as follows:
fish.mds<-metaMDS(fish.data, distance=
Hi Ali,
Check
?split
HTH,
Jorge
On Wed, Dec 1, 2010 at 2:50 PM, Ali S <> wrote:
> I have a matrix with 3 years of data (2006, 2009, 2010). I am trying to
> split this matrix by year so that I have 3 separate matrices. My matrix
> looks like this:
> Q16.1 Year Gender Grade1 3 2006
Try:
plot (myCatVariable)
Andrew Miles
Department of Sociology
Duke University
On Dec 1, 2010, at 2:51 PM, phoebe kong wrote:
Hi,
Can someone tell me how to draw a histogram for the following summary?
Richard Minnie Albert Helen Joe Kingston
1233 56 6715
Hi,
Can someone tell me how to draw a histogram for the following summary?
Richard Minnie Albert Helen Joe Kingston
1233 56 6715 66
The summary tell that Richard has occurrence 12, Minnie has occurrence 33,
and so on. I would like to view this summ
I have a matrix with 3 years of data (2006, 2009, 2010). I am trying to split
this matrix by year so that I have 3 separate matrices. My matrix looks like
this:
Q16.1 Year Gender Grade1 3 2006 1 52 2 2006 0
53 3 2006 0 54 3 2006 0 5
Tingting Zhan jefferson.edu> writes:
>
> Hi all,
>
> I have the output of summary() of an lme object called "lme.exp1", for
> example
>
> #
> > summary(lme.exp1)
>
[snip]
> for the common variance parameter sigma. But if I need the covariance
>
I think both responses so far have missed the point, (assuming the O was a typo
for zero).
That is:
> seq(0:1)
[1] 1 2
when
> seq(0,1)
[1] 0 1
was intended.
Ray Brownrigg
On Wed, 01 Dec 2010, Ista Zahn wrote:
> So you are warning us that you must type zero instead of the letter O
> when we want
Hi:
On Wed, Dec 1, 2010 at 7:49 AM, CALEF ALEJANDRO RODRIGUEZ CUEVAS <
alejandro.rodriguez.cue...@gmail.com> wrote:
> Hello dudes.
>
> I'm developing VAR analysis based on suggestions made by Horváth in its
> paper Canonical Correlation Analysis and Wiener-Granger Causality Tests.
>
Poor Horvath
Hi
On 2/12/2010 4:55 a.m., Mark Ebbert wrote:
Dear R Gurus,
I have a fairly simple problem, but I haven't been able to find the
answer on 'the google' or in the r-help archives.
I am generating plots on both Windows and OS X where I need to
guarantee that the font used is Arial. In my plot com
Dear R-users,
I'm trying to work out a way to set default values for arguments in a function
which includes the optional argument '...'.
In practice, I have a 'plot' method for a function which specifies different
types of plots. Every different plot should have different default arguments
(fo
Thank you sir for your answer. I'll take a look at the original legend
function and modify it for my own purpose.
With regards,
Phil
--
View this message in context:
http://r.789695.n4.nabble.com/Question-regarding-legend-look-tp3067466p3067894.html
Sent from the R help mailing list archive at
David Winsemius wrote:
>
>
> ...BUT NOT SHOWN US CODE
> OR SAMPLE DATA.
>
Hi and thank you for your help.
For instance, here's my code:
legend("topleft", inset = .05, title="Water masses", pch = c(22,25,21), lty
= c(4,1,2), lwd = 1, c("North","Central","South"), horiz = F, pt.bg
=c("gray
On Dec 1, 2010, at 1:12 PM, Charles C. Berry wrote:
On Wed, 1 Dec 2010, cmccar...@bmcc.cuny.edu wrote:
Hi,
Suppose I have the following data
name score
Abel88
Baker 54
Charlie77
stored a table called myData.
I want to write a function that will create a table which i
Thank you all for your very fast replies.
I tried Henrique's method (see one of the above posts) , and it works
perfectly!
Thanks again!
--
View this message in context:
http://r.789695.n4.nabble.com/How-to-pass-selection-criteria-in-a-function-tp3067765p3067829.html
Sent from the R help mailin
Hi Arthur,
I was asking the same thing and came across the following (your need the "sna"
library).
http://students.washington.edu/mclarkso/documents/gplot%20Ver2.pdf
Take a look at the edge.lwd and vertex.cex examples of the function gplot. You
can use vectors for the different nodes.
Kind r
Try this:
newTable <- function(data, criteria) {
do.call(subset, list(data, substitute(criteria)))
}
On Wed, Dec 1, 2010 at 3:56 PM, wrote:
> Hi,
> Suppose I have the following data
>
> name score
> Abel88
> Baker 54
> Charlie77
>
> stored a table called myData.
>
>
>
On Wed, 1 Dec 2010, cmccar...@bmcc.cuny.edu wrote:
Hi,
Suppose I have the following data
name score
Abel??? 88
Baker? 54
Charlie??? 77
stored a? table called myData.
I want to write a function that will create a table which is a subset of myData
containing those have a score > 7
Hi,
Suppose I have the following data
name score
Abel 88
Baker 54
Charlie 77
stored a table called myData.
I want to write a function that will create a table which is a subset of myData
containing those have a score > 75.
I know I can do this with the following command:
On 01/12/2010 10:13 AM, Filoche wrote:
Hi everyone.
I have a quick question regarding the look of my legend in my plot. As you
can see in the next figure, I have 3 series.
http://r.789695.n4.nabble.com/file/n3067466/legend.png
However, I find rather difficult to differentiate the series 1 and
On Wed, 1 Dec 2010, Anna Berthinussen wrote:
Hi,
I receive the following warning message when I run a poisson GLM in R:
"glm.fit: fitted rates numerically 0 occurred"
The model summary is shown below. The variable 'Species' consists of
counts of different species ranging from 0 to 4. I suspec
Try this:
id <- 1:20
grade <- c(4,4,4,5,5,7,7,7,7,8,8,8,9,9,9,9,9,10,10,10)
sequence <- ave( id, grade, FUN=seq )
# if grade is not sorted
grade2 <- sample(grade)
sequence2 <- ave( id, grade2, FUN=seq )
cbind( grade2, sequence2 )
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermo
dear M.,
I do not know how to get the SE for the joinpoint (or breakpoint) from
your ljr fit. However you can find useful the segmented package which
works for any GLM (including the logistic one) and it returns
(approximate) StErr (and Conf Int) also for the joinpoint (breakpoint in
the segme
Hello
Anyone know how can I calculate the value of the beta parameter when I know
the number of cointegrating relationships between two variables.
I mean, I using the procedure: ca.jo I do the following:
summary (ca.jo (UR [, c (2.52)], type = "trace" ECDET = "trend", K = 2, spec
= "longrun"))
On Dec 1, 2010, at 11:51 AM, David Winsemius wrote:
On Dec 1, 2010, at 10:13 AM, Filoche wrote:
Hi everyone.
I have a quick question regarding the look of my legend in my plot.
As you
can see in the next figure, I have 3 series.
http://r.789695.n4.nabble.com/file/n3067466/legend.png
On Dec 1, 2010, at 10:13 AM, Filoche wrote:
Hi everyone.
I have a quick question regarding the look of my legend in my plot.
As you
can see in the next figure, I have 3 series.
http://r.789695.n4.nabble.com/file/n3067466/legend.png
However, I find rather difficult to differentiate the s
On Wed, 2010-12-01 at 21:23 +0800, yangwenjing wrote:
> Dear everyone,
>
> I am doing some work about species richness estimation. Nonparametric
> estimation (such as Chao1, Jacknife1) can be done just using function
> "specpool()" and "estimateR()" in package "vegan". The problem is that
> I can
Luana -
It's probably not the most efficient way, but here's
a solution that's not dependent on the grades being sorted:
grade <- c(4,4,4,5,5,7,7,7,7,8,8,8,9,9,9,9,9,10,10,10)
unlist(sapply(rle(grade)$lengths,function(x)seq(1,x)))
[1] 1 2 3 1 2 1 2 3 4 1 2 3 1 2 3 4 5 1 2 3
Hi Luana,
Try this:
ID <- 1:20
grade <- c(4, 4, 4, 5, 5, 7, 7, 7, 7, 8, 8, 8, 9, 9, 9, 9, 9, 10, 10, 10)
d <- data.frame(ID, grade)
d$Sequence <- do.call(c, sapply(rle(grade)$lengths, seq))
d
HTH,
Jorge
On Wed, Dec 1, 2010 at 11:08 AM, Luana Marotta <> wrote:
> Hello fellows,
>
> I would like
You can also explore the classical() function in the MiscPsycho package that
does item analysis.
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
> Behalf Of Iasonas Lamprianou
> Sent: Wednesday, December 01, 2010 4:52 AM
> To: Michael Be
Hi,
I receive the following warning message when I run a poisson GLM in R:
"glm.fit: fitted rates numerically 0 occurred"
The model summary is shown below. The variable 'Species' consists of
counts of different species ranging from 0 to 4. I suspect this may
have something to do with the warn
Try this:
> ord <- order(grade)
> ID <- Id[ord]
> grade <- grade[ord]
> sequence <- unlist(sapply(table(grade), FUN = function(x) 1:x),
use.names = F)
And as a general tip, it is much easier to work with related values like
ID and grade if they are in a data frame. Such as:
> dat <- data.frame
Here are two options that might work for you given the little bit you've said:
## If its the same parameter all 30 times
## say, for example, base = 4.5 to log
for(i in 1:30) {
print(log(1:10, base = 4.5))
}
## if they are different parameters, you could try
lapply(X = c(1.3, 3, 2.2, 4, 5), FUN
Hi everyone.
I have a quick question regarding the look of my legend in my plot. As you
can see in the next figure, I have 3 series.
http://r.789695.n4.nabble.com/file/n3067466/legend.png
However, I find rather difficult to differentiate the series 1 and 3
according to their line type (lty).
Dear R Gurus,
I have a fairly simple problem, but I haven't been able to find the answer on
'the google' or in the r-help archives.
I am generating plots on both Windows and OS X where I need to guarantee that
the font used is Arial. In my plot command I specify 'fontfamily="Arial"'. The
probl
On Wed, Dec 1, 2010 at 9:02 AM, Ben Bolker wrote:
>
>
> Peter Ehlers ucalgary.ca> writes:
>
>>
>> It might be a good idea not to use an outdated version of R.
>> I don't see your "problem" in R 2.12.0.
>>
>> Peter Ehlers
>>
>> On 2010-12-01 05:44, Jean.Coursol math.u-psud.fr wrote:
>> > Just for
Hi all,
I have the output of summary() of an lme object called "lme.exp1", for
example
#
> summary(lme.exp1)
Linear mixed-effects model fit by REML
Data: DATA
Log-restricted-likelihood: -430.8981
Fixed: fixed.exp1
Random effects:
Formu
Sorry this should have ben to the whole list:
Hadley
I think I've sorted it out in my head but for the record, and just to be sure...
I guess what I was expecting was that the width parameter would be independent
of binwidth. Thus a width parameter of 0.5 would always indicate an overlap of
hal
Hello fellows,
I would like to create a sequence for repeated numbers in a dataset. For
example:
ID <- c(1:20)
grade <- c(4,4,4,5,5,7,7,7,7,8,8,8,9,9,9,9,9,10,10,10)
Data:
ID Grade
1 4
2 4
3 4
4 5
5 5
6 7
7 7
8 7
9 7
(...)
I would like to create a variable "sequence":
Data:
In specifying a CFA model using the sem package, I got the following warning
message:
In sem.default(ram = ram, S = S, N = N, param.names = pars, var.names =
vars, :
Could not compute QR decomposition of Hessian.
Optimization probably did not converge.
This is the complete input (incl
On Dec 1, 2010, at 9:18 AM, Peter Ehlers wrote:
It might be a good idea not to use an outdated version of R.
I don't see your "problem" in R 2.12.0.
Furthermore isn't it a bit to be asking why the sign on a
number that is effectively zero happens to be negative? I would think
that FAQ 7
Hello dudes.
I'm developing VAR analysis based on suggestions made by Horváth in its
paper Canonical Correlation Analysis and Wiener-Granger Causality Tests.
That's the reason I'm looking for if there's any R package to develop Wiener
- Granger Causality Test.
Thanks a lot for your unvaluable he
Isn't it what do.call() does?
Ivan
Le 12/1/2010 16:39, Lamke a écrit :
Thank you so much Joshua. That's exactly what I am looking for.
What I wanted to do is to pass a parameter to a function and I have to run
the functions 30 times. Instead of typing them all out, I created a long
string of
Thank you so much Joshua. That's exactly what I am looking for.
What I wanted to do is to pass a parameter to a function and I have to run
the functions 30 times. Instead of typing them all out, I created a long
string of "f(a);f(b);f(c) ..." using paste() and use eval and parse to
evaluative t
I am having problems trying to get R to graph data input that is log-normal
on the horizontal (x) axis.
The data is log (base 10), and I am more interested in viewing the tails of
the distribution. The closest I can get with this is log on the vertical
(y) axis and linear on the horizontal axis.
Hi,
A new package Rd2roxygen has been released on CRAN:
http://cran.r-project.org/package=Rd2roxygen (source package
available; binaries to come)
This package can be useful for developers who used to document their
functions in the raw Rd files but want to switch to roxygen now using
inline roxyg
Peter Ehlers ucalgary.ca> writes:
>
> It might be a good idea not to use an outdated version of R.
> I don't see your "problem" in R 2.12.0.
>
> Peter Ehlers
>
> On 2010-12-01 05:44, Jean.Coursol math.u-psud.fr wrote:
> > Just for fun
> >
> > This small program gives negative Sum of Sq in a
thanks!
very useful
--
View this message in context:
http://r.789695.n4.nabble.com/Plot-a-matrix-recursively-tp3067283p3067443.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mail
Hi Dimtris and esteemed useRs,
I don't understand why i get this error message when attempting to use
merge() -
> temp <- merge(x, y[,17, drop=FALSE], by=rownames, sort=FALSE)
Error in as.vector(x, mode) :
cannot coerce type 'closure' to vector of type 'any'
It should work because:
> all(r
> However if you do:
> ggplot(data=dafr, aes(x = d1, fill=d2)) + geom_histogram(binwidth = 1,
> position = position_dodge(width=0.99))
>
> The position of first bin which goes from 0-2 appears to start at about 0.2
> (I accept that there is some "white space" to the left of this) while the
> pos
On Wed, 2010-12-01 at 07:32 -0500, Gabor Grothendieck wrote:
> On Sun, Nov 28, 2010 at 6:44 AM, Stephen Liu wrote:
>
> In the end it appears that there was no problem in accessing the gdata
> package or any of its components and that the poster's problem stemmed
> from issuing this series of com
On 12/01/2010 05:57 AM, nquer...@clinic.ub.es wrote:
> Hi everybody,
>
>
>
> I am a beginner in the steps of pre-processing and data analysis of
> non-targeted metabolomic profiling experiments.
> Anyone knows if there exists some tool for normalizing this type of data (raw
> data or XCMS
barbara.chaves wrote:
>
>
> I have a dataset with 4 subjects (see ID in example), and 4 treatment (see
> TRT in example) which are tested on 2 locations and in 3 blocs. By using
> Lattice dotplot, I made a graph that shows the raw data per location and
> per bloc. In that graph, I would like
Hi He Zhang,
>> Is the following right for extracting the scores?
>> ...
>> pca$loadings
>> pca$score
Yes.
But you should be aware that the function principal() in package psych is
standardizing your data internally, which you might not want. That is, the
analysis is being based on the correla
It might be a good idea not to use an outdated version of R.
I don't see your "problem" in R 2.12.0.
Peter Ehlers
On 2010-12-01 05:44, jean.cour...@math.u-psud.fr wrote:
Just for fun
This small program gives negative Sum of Sq in anova (with R versions
2.9.2 and 2.10.0) (linux 32 kernel 2.6.29
On 2010-12-01 04:02, barbara.cha...@bayer.com wrote:
Dear,
I have a dataset with 4 subjects (see ID in example), and 4 treatment (see
TRT in example) which are tested on 2 locations and in 3 blocs. By using
Lattice dotplot, I made a graph that shows the raw data per location and
per bloc. In tha
On 12/01/2010 02:43 PM, alcesgabbo wrote:
I plot the first column with the following function:
plot(m[,1],type="o", xaxt="n",ylim=c(min(m[,1:length(colnames(m))])-1,
max(m[,1:length(colnames(m))])+1))
for the other columns I use there functions:
lines(m[,2],type=\"o\")
lines(m[,3],type=\"o\")
Take a look at
?matplot
HTH,
Gerrit
On Wed, 1 Dec 2010, alcesgabbo wrote:
Hi,
I have the following matrix (named m):
key
sensor_date Laser_1 Laser_2
Laser_3
2010-09-30T15:00:12+02006
Hi everybody,
I am a beginner in the steps of pre-processing and data analysis of
non-targeted metabolomic profiling experiments.
Anyone knows if there exists some tool for normalizing this type of data (raw
data or XCMS matrix data)
in R repositories? Many thanks in advance.
Best r
Just for fun
This small program gives negative Sum of Sq in anova (with R versions
2.9.2 and 2.10.0) (linux 32 kernel 2.6.29.6-smp slackware 13.0) :
y = c(3.6, 5.0, 5.0, 4.6, 4.5, 4.3, 4.5, 5.1, 4.5, 4.3)
trans = as.factor(c("NT","NT","NT","NT","NT","T","T","T","T","T"))
lc = lm(y ~ trans)
l1
> David Winsemius
> on Mon, 29 Nov 2010 09:46:08 -0500 writes:
> On Nov 29, 2010, at 9:00 AM, pilchat wrote:
>> Hi guys,
>>
>> to make it easier, here is a simple case with the same issues. I use
>> the short function below to make the attached PS file.
>>
Dear everyone,
I am doing some work about species richness estimation. Nonparametric
estimation (such as Chao1, Jacknife1) can be done just using function
"specpool()" and "estimateR()" in package "vegan". The problem is that I can
not found any functions for parametric estimation (such as MMMe
Hi,
I have the following matrix (named m):
key
sensor_date Laser_1 Laser_2
Laser_3
2010-09-30T15:00:12+020063
1
2010-10-31T15:05:07+0100
On 2010-12-01 01:51, Seth Roberts wrote:
coplot() usually puts grid lines in the panels it makes. To see examples,
example(coplot).
How can I remove those grid lines?
>
> Seth Roberts
>
That's hard-coded in coplot. But you can easily make a
modified copy of coplot. Here are the steps:
1. make
Hi Sandy,
I don't have an answer to your question, but wanted to share another
possibility:
ggplot(data=dafr, aes(x = d1, fill=d2)) + geom_histogram(binwidth = 1,
position = position_identity(), alpha=.5)
Best,
Ista
On Wed, Dec 1, 2010 at 6:09 AM, Small Sandy (NHS Greater Glasgow &
Clyde) wrote
Peter Langfelder wrote:
>
> On Fri, Nov 26, 2010 at 6:55 AM, Derik Burgert wrote:
>> Dear list,
>>
>> running a hierachical cluster analysis I want to define a number of
>> objects that build a cluster already. In other words: I want to force
>> some of the cases to be in the same cluster from
Hi,
I am also doing PCA.
Is the following right for extracting the scores?
library(psych)
pca<-principal(data,nfactors=,rotate="varimax",scores=T)
pca$loadings
pca$score
Best regards,
He
On Tue, Nov 30, 2010 at 10:22 AM, Liviu Andronic wrote:
> Dear all
> I'm unable to find an example of extra
Hi Bill,
Please keep the R list copied.
I don't think you've given enough information to pin down your
problem. This works:
dat <- read.table(textConnection("dates returns
2000-1-4 -0.038344718
2000-1-50.00195
2000-1-60.000955702
2000-1-70.027090384
2000-1-100.011189966
2
On Sun, Nov 28, 2010 at 6:44 AM, Stephen Liu wrote:
> Hi folks,
>
> Win 7 64 bit
> R 32 bit
>
>> install.packages("gregmisc")
> Installing package(s) into
> ‘C:\Users\satimiswin764\Documents/R/win-library/2.12’
> (as ‘lib’ is unspecified)
> --- Please select a CRAN mirror for use in this session
On Wed, Dec 1, 2010 at 2:56 AM, Stephen Liu wrote:
> Hi folks,
>
> Win 7 64bit
> R 2.12.0 32bit
>
> Problem in reading Excel spreadsheets
>
There are quite a few alternatives for reading Excel spreadsheets.
They are listed here:
http://rwiki.sciviews.org/doku.php?id=tips:data-io:ms_windows
--
The following works nicely. Thank you.
plot(z, sf[[s]], type="l") #where the dataframe name in my case is "sf"
Sorry - I meant to say there were 24 rows in each column, so z=1:24
Gregory A. Graves, Lead Scientist
Everglades REstoration COoordination and VERification (RECOVER)
Restoration Scie
Dear,
I have a dataset with 4 subjects (see ID in example), and 4 treatment (see
TRT in example) which are tested on 2 locations and in 3 blocs. By using
Lattice dotplot, I made a graph that shows the raw data per location and
per bloc. In that graph, I would like to have a reference line per
On 01/12/2010 5:12 AM, Koray Kaya wrote:
Hello, I tried to use GEOQuery package of BioC. It does not download GSE. I investigated
problem, understood that the problem was about internal function "download".
Reccomendations about it mostly suggest switching any proxy off in R. I did, and nothing
On 30/11/2010 9:54 PM, randomcz wrote:
Hi guys,
How to pass an operator to a function. For example,
test<- function(a, ">", b)
{
return(a>b) #the operator is passed as an argument
}
Thanks,
It's much simpler than the other suggestions. Just pass the operator,
and treat it as a fun
[[alternative HTML version deleted]]
__
R-help@r-project.org 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,
Thanks Thierry
Using the position_dodge positioning option ought to work but there is
something wrong with the scaling when the binwidth is not = 1.
You can begin to see this with the example I sent:
If you do
ggplot(data=dafr, aes(x = d1, fill=d2)) + geom_histogram(binwidth = 1, position
= pos
Run
update.packages(checkBuilt=TRUE)
to get package binaries that were compield for your version of R.
Uwe Ligges
On 24.11.2010 04:24, john moran wrote:
Apologies for my previous effort in HTML which apparently was scrubbed
Dear R-users
I wonder if I could get advice on the above proble
Hi again,
There are also other packages to read and write xls files, including
gdata (read), dataframe2xls and writeXLS (write), which depend on Perl
or Python if I'm not mistaken.
I have no idea whether they work on Windows 64 bits.
Visit Crantastic http://crantastic.org/
Ivan
Le 12/1/2010
1 - 100 of 124 matches
Mail list logo