Hi again, this is a follow-up of my previous post :
I've noticed that
> attr(first.list[[1]]$size.ecdf,"call") gives ecdf(test.moms[, "m.pxs"]) and
> that there is no "test.moms" field in my records (a local variable named
> test.moms was used inside a function, but it was stored under another
On 24/10/2013 07:54, Olivier Eterradossi wrote:
Hi list,
In June 2013, using R.2.15.1 (i386) on Windows 7, I calculated a set of ecdf which I
stored in lists of records having a "size.ecdf" field, with following structure
:
Note that R 2.15.1 was not current even then.
[skip previous struc
We would like to announce the following stats course;
Course: Introduction to MCMC, Linear mixed effects models and GLMM with R
When: 20-24 January, 2014
Where: Bangor University, UK
Info: http://www.highstat.com/statscourse.htm
Flyer: http://www.highstat.com/Courses/Flyer2014_02Bangor.pdf
Kind
this works for me:
plot(1,main="\u2113")
cheers
Am 24.10.2013 01:39, schrieb Byron Dom:
>
>
> Original post: On 13/10/2013 18:53, Byron Dom wrote:
>
>>> Due to convention a script "l" - $$\ell$$ (LaTeX \ell) is used to
>>> represent a certain quantity in something I'm working on. I'm
>>> un
Thanks Jim, maybe now I can start learning.
Next I want to use my choice of n digits or width instead of round.
--
View this message in context:
http://r.789695.n4.nabble.com/labeling-abscissa-using-a-function-of-the-plotted-scale-tp4678075p4678946.html
Sent from the R help mailing list archiv
Dear, Patrick.
I'm using a workaround to work with RWeka and multicore package. The most
important point is not loading the RWeka package and use the namespace in a
encapsulated function. Take a look on this code:
# C4.5 classifier. Return the prediction for a test dataset.
cl.c45 = function(tra
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
> project.org] On Behalf Of wayne.w.jo...@shell.com
> Sent: 22 October 2013 18:52
> To: R-help@r-project.org
> Subject: [R] nls model definition help
>
> Hi fellow R users,
>
> I'm trying to fit a model u
Hi Arun,
That seemed to do the trick - thanks!!
Jonathan
On Wed, Oct 23, 2013 at 11:12 PM, arun wrote:
> HI,
>
> Better would be:
> res1 <- dcast(df,gene~case,value.var="issue",paste,collapse=",",fill="0")
>
> str(res1)
> #'data.frame':2 obs. of 4 variables:
> # $ gene : chr "gene1"
Dear all,
I was wondering whether (or better: how) I can use R to read recursively a
directory to get all the sub-folders and files located in the root folder and
put it into a tree like structure where the leaves are files and intermediate
nodes are the directories? The idea is that I'd like t
Hi all,
I am trying to plot a raster object (I can explain why but the point is that it
would be a raster objeçt)..
I have selected a small code to show you exactly the problem
require(raster)
test<-matrix(data=runif(1),nrow=100)
m<-raster(test)
plot(m,axes="FALSE")
axis(1,at=c(0,1),labels=
A wild guess -- take a look at the CRAN "phylohenetics" task view, as
that sounds like the sort of thing that might have tree generation and
manipulation functions.
... but you may do better with some non-R tool out there.
(Hopefully, you'll get a better response, though).
Cheers,
Bert
On Thu,
Thanks Jim, maybe now I can start learning.
Here is a run of my trying to learn:
> xvalue<-c(5.2,1.3,9.7,2.8,8.1,4.7,6.6,7.4)
> yvalue<-c(9,3,4,7,2,5,3,6)
> plot(xvalue,yvalue)
> axis(1,at=NULL,labels=1/xvalue,digits=5)
Error in axis(1, at = NULL, labels = 1/xvalue, digits = 5) :
'labels' is supp
One R package that might be of interest would be 'diagram':
http://cran.r-project.org/web/packages/diagram/
I would also agree with Bert here and would point you in the direction of
PSTricks, which can handle these sorts of complex figures. It would of course
require learning LaTeX, but that
perhaps as a somewhat tamer guess than Bert's, it's probably not
very hard to write a function that pieces together a dendrogram
(S3 class provided by the stats package), using the list.dirs and
list.files functions.
This could then be plotted via the plot function, and you could
even designate at
This is by no means an explanation, but as a possible approach to consider:
plot(0:1, 0:1, xaxt='n', asp=1)
axis(1,at=c(0,1),labels=c("a","b"))
plot(m,add=TRUE)
-Don
--
Don MacQueen
Lawrence Livermore National Laboratory
7000 East Ave., L-627
Livermore, CA 94550
925-423-1062
On 10/24/
I'm writing a function that needs an optional variable. If the variable is
given, then a subset must be made using that variable, but if the variable
is not given the subset must be done using all the values for that variable
(ignoring that variable).
Something like this:
*band <- function(file,
Dear all,
Running the following code produces an error in dummy.coef. I guess it might be
a bug..
> ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)
> trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)
> group <- gl(2, 10, 20, labels = c("Ctl","Trt"))
> weight <- c(ctl, trt)
Dear R-users,
I am trying to analyze data with ordinal independent variables and
binary response . Could some one suggest me an appropriate statistic
for this kind of analysis.
Thanks,
Regards,
Kishir
__
R-help@r-project.org mailing list
https://stat.e
Hi professors
I'm finding the best package in R for manage a big database about vessels. I
need to divided the travel of each boat into tracks. I have GPS points.
I would join these points in a line on the map. I make that with the function
"lines", but I can't divided by vessel or by track,
I
Berend Thanks. I am aware of the missing kb's in the download. However,
my problem is that it only happens with this package. I download other
packages manually and via command-line install.packages() and
Do not have any problems with completing full downloads of the entire file.
Perh
Hi All
I've gotten some awesome help getting a formular that finds the intersection of
two vectors. This works brilliantly, but I can't figure out how to make it run
over another factor. A simple example looks likes this:
df <- data.frame(
id = factor(rep(c("supply", "demand"), each =
On 24-10-2013, at 16:15, David Mora wrote:
> I'm writing a function that needs an optional variable. If the variable is
> given, then a subset must be made using that variable, but if the variable
> is not given the subset must be done using all the values for that variable
> (ignoring that vari
Kishor,
It sounds like you will need to use logistic regression. Although the following
URL might help you better understand logistic regression,
http://ww2.coastal.edu/kingw/statistics/R-tutorials/logistic.html it would
probably be most helpful if you would see help from a local statistician.
It is difficult to advise you when you don't provide reproducible data.
For example, you could include the first 6 lines of your data using
dput(head(mydata))
If you have GPS points stored as easting, northing, and if you have a
variable that identifies the vessel, and if your data are sorte
On Oct 23, 2013, at 7:53 PM, Long Vo wrote:
> Hi R users,
> Currently I want to fit a FIGARCH model to a dataset. The only package that
> allow for it that I could find is fGarch. However it seems that the FIGARCH
> model class fitting of this package has been moved to Oxmetrics. I tried to
> ins
Dear John,
Thanks for the information. It is quite helpful.
Regards,
Kishor
On Thu, Oct 24, 2013 at 9:24 PM, John Sorkin
wrote:
> Kishor,
> It sounds like you will need to use logistic regression. Although the
> following URL might help you better understand logistic regression,
> http://ww2.co
Using your csv file, I get no error running the following code.
SPIV2 <- read.csv("C:/Temp/SPI V2.csv")
heatmap(as.matrix(SPIV2), na.rm=TRUE)
Perhaps you modified SPIV2 in some way prior to running heatmap()?
Jean
On Tue, Oct 22, 2013 at 10:29 AM, David Romero wrote:
> Hi Jean
>
Sorry if this subject has been already dealt here.
Which are some common tests for nonparametric k-way ANOVA?
I have read about Kruskal-Wallis test as a kind of nonparametric one-way
ANOVA, but I have not found anything about a general-setting (I mean k-way)
nonparametric ANOVA.
Can you recommen
Have you tried to download from another CRAN mirror, e.g.
http://cran.r-project.org/web/packages/XML/
More mirrors at http://cran.r-project.org/mirrors.html
As already others said, it's very unlikely that this is not an issue
on your end.
/Henrik
On Thu, Oct 24, 2013 at 4:54 AM, Steven Dway
On Oct 24, 2013, at 11:38 AM, David Winsemius wrote:
>
> On Oct 23, 2013, at 7:53 PM, Long Vo wrote:
>
>> Hi R users,
>> Currently I want to fit a FIGARCH model to a dataset. The only package that
>> allow for it that I could find is fGarch. However it seems that the FIGARCH
>> model class fit
Rick,
This uses a different approach than what you propose, but it gets the job
done and perhaps you will find it helpful.
Jean
library(jpeg)
# download image from internet for use in example
# http://www.tootsie.com/wallpaper/wp_concord2_1344x1008.jpg
img <- readJPEG("c:/temp/wp_concord2_1344x1
On Oct 24, 2013, at 18:59 , Vicent Giner-Bosch wrote:
> Sorry if this subject has been already dealt here.
>
> Which are some common tests for nonparametric k-way ANOVA?
>
> I have read about Kruskal-Wallis test as a kind of nonparametric one-way
> ANOVA, but I have not found anything about a g
On 24/10/2013 18:25, Marc Schwartz wrote:
On Oct 24, 2013, at 11:38 AM, David Winsemius wrote:
On Oct 23, 2013, at 7:53 PM, Long Vo wrote:
Hi R users,
Currently I want to fit a FIGARCH model to a dataset. The only package that
allow for it that I could find is fGarch. However it seems that
Thanks. That works for me too.
I had worked that out based on Brian's response to my original post.
If you read thru my response to him in detail, you'll see:
> Here are a couple of examples using it that worked:
> > plot(1:10,xlab="\u2113")
> >
> plot(1:10,xlab="\u2113(\u2113 + 1)")
On Oct 24, 2013, at 1:11 PM, Prof Brian Ripley wrote:
> On 24/10/2013 18:25, Marc Schwartz wrote:
>>
>> On Oct 24, 2013, at 11:38 AM, David Winsemius wrote:
>>
>>>
>>> On Oct 23, 2013, at 7:53 PM, Long Vo wrote:
>>>
Hi R users,
Currently I want to fit a FIGARCH model to a dataset.
Obviously row 11 is incomplete.
It may contain nothing more than a space or an EOL.
Using the "fill=TRUE" argument to 'scan' or 'read.table' is your friend
here.
--
View this message in context:
http://r.789695.n4.nabble.com/Error-in-scan-tp4678976p4678982.html
Sent from the R help mailing li
Dear R users,
Can I add a text to a plot without x axis?
I saw a function that enables placing the text with mouse
text(location, "text to place", pos, ...)
locator(1)
but it did not work out.
--
Ahmed M. Attia
Research Assistant
Dept. Of Soil&Crop Sciences
Texas A&M University
ahmed.at...
Hi,
Try:
fun1 <- function(dat){
mat1 <- combn(colnames(dat1),2)
res <- sapply(seq_len(ncol(mat1)),function(i) {x1<- dat[,mat1[,i]];
wilcox.test(x1[,1],x1[,2])$p.value})
names(res) <- apply(mat1,2,paste,collapse="_")
res
}
set.seed(432)
dat1 <- as.data.frame(matrix(sample(18*10,18*10,replace=FALS
David,
You should cc r-help on all correspondence so others can follow the thread.
The error appears to be a result of the missing values yielded by the
dist() function on your matrix.
See https://stat.ethz.ch/pipermail/r-help/2009-July/204443.html.
Jean
On Thu, Oct 24, 2013 at 1:25 PM, David
Dear R users,
How can I customize x axis in a boxplot?
xaxt="n" to suppress x axis, but I would like to customize x axis for
different groups.
Thanks
--
Ahmed M. Attia
Research Assistant
Dept. Of Soil&Crop Sciences
Texas A&M University
ahmed.at...@ag.tamu.edu
Cell phone: 001-979-248-5215
___
Depends on what the original poster really wants with a "nonparametric
k-way ANOVA". If what they are really looking for is an approach that has
fewer distributional assumptions, then it is possible to perform
permutation tests for linear models parameterized for k-way anova that
would eliminated
You should provide an example of what you are trying to do. Try
this
> plot(0, 1, xlab="", xaxt="n")
> text(locator(1), "This is my text", xpd=TRUE)
The xpd= parameter controls the clipping region. The default
xpd=FALSE clips text to the plot region. You were not clear, but
if you were trying to
On Oct 24, 2013, at 7:00 AM, Hurr wrote:
> Thanks Jim, maybe now I can start learning.
> Here is a run of my trying to learn:
>> xvalue<-c(5.2,1.3,9.7,2.8,8.1,4.7,6.6,7.4)
>> yvalue<-c(9,3,4,7,2,5,3,6)
>> plot(xvalue,yvalue)
>> axis(1,at=NULL,labels=1/xvalue,digits=5)
> Error in axis(1, at = NULL
?axis
-
David L Carlson
Associate Professor of Anthropology
Texas A&M University
College Station, TX 77840-4352
-Original Message-
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On Behalf Of Ahmed Attia
Sent: Thursday, October
Hello,
There's a bug in your function, it should be 'dat', not 'dat1'. In the
line marked, below.
fun1 <- function(dat){
mat1 <- combn(colnames(dat),2) # Here, 'dat' not 'dat1'
res <- sapply(seq_len(ncol(mat1)),function(i) {x1<- dat[,mat1[,i]];
wilcox.test(x1[,1],x1[,2])$p.value})
Hi,
You may try:
?list.files()
nm1 <- list.files(pattern=".txt")
res <- lapply(nm1,function(x) {
ln1 <- readLines(x)
indx1 <- grep("DATE PROCESSED",ln1)
indx2 <- grep("[A-Z]",ln1)
There are several packages for track data on CRAN and they don't all
force the data to have "specific variables" beyond position,
date-time, track ID. With trip you can do this easily enough, but
depending on your needs adehabitatLT or move might be better. I wrote
trip so I know that you can do t
A good choice might be to store your data in objects of class
SpatialLinesDataFrame
Then you can assign any attributes you want to each line, such as your
vessel, type, destination, or status.
This will make it relatively easy to, for example, plot subsets of tracks
specified by different attr
Hi,
You may try
?rollmean from library(zoo)
?SMA from library(TTR)
dat1 <- structure(...
library(zoo)
res1 <- rollmean(dat1[-c(1:5),2],8)
library(TTR)
res2 <- SMA(dat1[-c(1:5),2],8)
all.equal(res1,res2[!is.na(res2)])
#[1] TRUE
A.K.
Hi,
How can I calculate running 8-hour averages of temperat
I agree that this seems to be a bug, but OTOH the help for dummy.coef()
*does* say:
> The method used has some limitations, and will give incomplete results
> for terms such as |poly(x, 2)|.
(Note that it says it will "give incomplete results" not "throw an error"!)
Looking into the code of d
On 10/25/2013 01:00 AM, Hurr wrote:
Thanks Jim, maybe now I can start learning.
Here is a run of my trying to learn:
xvalue<-c(5.2,1.3,9.7,2.8,8.1,4.7,6.6,7.4)
yvalue<-c(9,3,4,7,2,5,3,6)
plot(xvalue,yvalue)
axis(1,at=NULL,labels=1/xvalue,digits=5)
Error in axis(1, at = NULL, labels = 1/xvalue,
Thanks David, did I use incorrect english to call "labels = 1/xvalue" a
formula or function?
Seems to me an obvious first try would be to place the labels where
they would have been placed without the formula.
I haven't yet, but I will study from your suggestions.
--
View this message in conte
On Oct 24, 2013, at 1:51 AM, Alexandra Kuznetsova wrote:
> Dear all,
>
> Running the following code produces an error in dummy.coef. I guess it might
> be a bug..
>
>> ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)
>> trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)
>>
On Oct 24, 2013, at 3:56 PM, Hurr wrote:
> Thanks David, did I use incorrect english to call "e" a
> formula or function?
Incorrect R terminology. A function in R is a specific sort of object and a
formula is a different sort of object. R terminology for the use of "labels =
1/xvalue" inside
You could store your first approach in a function and lapply it to
your by_hour variable:
df <- data.frame(
hour = factor(rep(1:5,4)),
id = factor(rep(c("supply", "demand"), each = 10)),
price= c(5,7,9,11,13,15,17,19,21,23,
20,18,16,14,12,10,8,6,4,2 ),
quantity = c(3,5,7,13,19,31
Dear all,
Running the following code produces an error in dummy.coef. I guess it might be
a bug..
> ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)
> trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)
> group <- gl(2, 10, 20, labels = c("Ctl","Trt"))
> weight <- c(ctl, trt)
Hi Jean,
Problem result: the 1st column value must be complete for heatmap to run
correctly.
This solved the same problem I had with vegdist.
Thank you for your help.
Regards
David
De: Adams, Jean [mailto:jvad...@usgs.gov]
Enviado el: 24 de octubre de 2013 14:12
Para: David Ro
Sorry for the html code:
I'm writing a function that needs an optional variable. If the variable is
given, then a subset must be made using that variable, but if the variable
is not given the subset must be done using all the values for that variable
(ignoring that variable).
Something like this:
Hi,
I'm trying to install R on CentOS 6.4.
Following some instructions online, I've done this:
rpm -Uvh
http://download.fedoraproject.org/pub/epel/6/i386/epel-release-6-8.noarch.rpm
yum install R
But yum fails, with this (full output below):
Error: Package: R-core-2.10.0-2.el5.x86_64 (Rocks-6
I'm trying to place labels so as to avoid covering points (and other labels) in
an ordination plot. I've been trying to use orditkplot() for this purpose, but
get an error message, even when I try to replicate the example:
>require("vegan")
>data(varespec)
> ord <- cca(varespec)
> orditkplot(o
Hi,
Check out this function:-
pairwise.wilcox.test {package=stats}.
example(pairwise.wilcox.test)
On Fri, Oct 25, 2013 at 2:15 AM, Rui Barradas wrote:
> Hello,
>
> There's a bug in your function, it should be 'dat', not 'dat1'. In the line
> marked, below.
>
> fun1 <- function(dat){
> m
It looks much better than mine.
with p value adjustment:
p.adjust(fun1(dat1), method = "holm", n = 153)
#
dat1$id <- 1:10
library(reshape2)
dat2 <- melt(dat1,id.var="id")
with(dat2,pairwise.wilcox.test(value,variable))
with(dat2,pairwise.wilcox.test(value,variable,p.adj="none"))
A.K.
Hello,
Do you have headers installed (generally tcl-devel and tk-devel)?
Regards,
Pascal
On 25 October 2013 10:33, Michael Stauffer wrote:
> Hi,
>
> I'm trying to install R on CentOS 6.4.
>
> Following some instructions online, I've done this:
>
> rpm -Uvh
> http://download.fedoraproject.org/p
On 25/10/2013 02:33, Michael Stauffer wrote:
Hi,
I'm trying to install R on CentOS 6.4.
This is not the right list. But
- As the posting guide says, we only support current R here. R 2.10.0
is ancient, and other people seem to have found 3.0.1 RPMs for Centos 6.3.
- It seems your RPM is
Thank you for the reply!
I need to use it in the package for constructing the contrast matrices.
Renaming a variable like y2 <- y^2 will not be a solution for me...
Will try to look whether the function dummy.coef can be modified to work this
issue.
With best regards,
Alexandra
-Original
65 matches
Mail list logo
