Hi
Following on David's rate argument
try (with modifications of pch and grid)
rate <- 1/4
shape = 8
rate = c(rep(1/4,100),rep(1/3,100))
x = rgamma(200,shape,rate)
groups = gl(2,100,200,labels=LETTERS[1:2])
dat = data.frame(x=x, gp=groups)
qqmath(~ x, data = dat,
groups = gp,
pc
I have a complicated, multi-part question. My apologies if I do not make
myself clear. I am also a fairly novice R user, so forgive me if this seems
rudimentary. I want to calculate a index of colocation for whale dive data
and prey distribution data. This entails:
Calculating a frequency distribu
Hello,
I've recently upgraded to R 2.11.1 on Mac OS X 10.8.4. Now when I
start R up I get the following message:
Error in identical(call[[1L]], quote(doTryCatch)) :
7 arguments passed to .Internal(identical) which requires 5
Error in normalizePath(dirname(pkgpath), "/", TRUE) :
3 argumen
"Vallejo, Roger" asked:
I would like to know if we can estimate Rg between two binary traits
(disease status: alive vs. dead) with the R package.
My data: we have 100 full-sib (FS) families,
and two random samples (each with n= 200 FS fish) from each FS family
were evaluated for disease resis
Dear pieter,
On Fri, 21 Jun 2013 10:12:26 +0200
pieter wrote:
> Hey,
> I don't really get the output of the performed Levene's test.
> So in short.
> two data sets, one with 75 elements, one with 283.
> And I want to check whether they have the same variance.
>
> I did a LeveneTest
> And this i
Hi,
May be this helps:
###Added more lines of fake data
fn_hp<- read.table(text="
id date wl_m wet cuml_day
585 fn 2012-03-03 0.1527048 1 1
586 fn 2012-03-04 0.2121408 1 2
587 fn 2012-03-05 0.1877568 1 3
On 22/06/13 01:48, Stefano Sofia wrote:
Dear R users,
I have a seasonal time series of period 4 (my_ts).
I would like to apply to my_ts a Seasonal ARMA(2,0) process (only the seasonal
part), with Seasonal AR coefficients respectively 0.54 and 0.5.
I tried to use the arima command from the stats
In a similiar manner with additions
dat <- data.frame(x=x,y=y,z=z)
dotplot(y~z+x,data = dat, pch = 16, type= c("p","g"))
HTH
Duncan
Duncan Mackay
Department of Agronomy and Soil Science
University of New England
Armidale NSW 2351
Email: home: mac...@northnet.com.au
At 07:00 22/06/2013, yo
On Jun 21, 2013, at 1:21 PM, Saptarshi Guha wrote:
> Hello,
>
> Consider the following transcript
>
>> x=NULL
>> x$date=10
>> x
> $date
> [1] 10
>
> or
>
>> x=NULL
>> x[10]=10
>> x
> [1] NA NA NA NA NA NA NA NA NA 10
>
>
> Wouldn't one expect that x remain NULL despite the further additions
Maybe something like this?
dataset <- data.frame(x=c(y, y), y=c(x, z), g=rep(1:2, each=58))
dotplot(x~y, groups=g, data=dataset)
-
David L Carlson
Associate Professor of Anthropology
Texas A&M University
College Station, TX 77840-4352
-Original Message---
Well, feel free to make up any semantics that you like for a language
you write, but R already has its own and tells you what to expect:
>From ?NULL:
NULL can be indexed (see Extract) in just about any syntactically
legal way: whether is makes sense or not, the result is always NULL.
Objects with
dotplot is a lattice function. "add" is an argument to some base
graphics. Never the twain shall meet. There is no "add" argument to
dotplot -- did you read the lattice Help for dotplot (found under
xyplot)??
I don't understand what you mean by "overlay" dotplots, but lattice
does this via conditi
Hello,
Consider the following transcript
> x=NULL
> x$date=10
> x
$date
[1] 10
or
> x=NULL
> x[10]=10
> x
[1] NA NA NA NA NA NA NA NA NA 10
Wouldn't one expect that x remain NULL despite the further additions (i.e
the x$date, x[10]) etc? Is coercing appropriate here?
Cheers
Saptarshi
Dear All,
wonder if you would provide your insights on the following: the code:
library(lattice)
y <-c(1:58)
x <-runif(58,5,10)
z <-runif(58,8,12)
dataset <-data.frame(y,x)
dotplot(y ~ x, data = dataset)
dataset <-data.frame(y,z)
dotplot(y~z,data = dataset,col="red")
I would like to overlay the
On Jun 21, 2013, at 3:39 AM, Stefan Sobernig wrote:
> Dear list,
>
> I have already posted to sci.stat.consult, however, I was hoping that there
> might be some experience in the R community on the following question:
>
>> I am wondering whether anyone has experience with variants of (parametr
Ah, I see.
You could try something like this ... create a list of your data at the
outset, then just use the for() loop to create the plots.
Jean
# fake data with two Seamounts ... one with a missing depth
SChla <- structure(list(Seamount = c("Atlantis", "Atlantis", "Atlantis",
"Atlantis", "A
Hi, Christian.
It sounds like to the maintainer feels that the package needs improvements,
i.e., it may not function as expected. However, you may contact him
directly using maintainer("randomForest").
Best,
Benjamin
On Fri, Jun 21, 2013 at 9:36 AM, Christian Kampichler <
christian.kampich...@
Hi
What about delivering further information about data, code, R and OS version.
Most probably you are tryind to index some object in cycle above its dimension.
Regards
Petr
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
> project.org] On Behalf Of
HI,
May be this also works:
plot(x,col=as.numeric(factor(x)))
A.K.
- Original Message -
From: Rui Barradas
To: Suparna Mitra
Cc: r-help@r-project.org
Sent: Friday, June 21, 2013 11:32 AM
Subject: Re: [R] How to define desired numbers to a vector based on the
present numbers
Hello,
Does it work if you ignore NAs when calculating the max?
ylim=c(0,max(sm.mean$Total_Ch + 0.1, na.rm=TRUE))
Jean
On Fri, Jun 21, 2013 at 7:30 AM, Tinus Sonnekus wrote:
> Hi all,
>
> So this code works for three of the six seamounts. The problem I am
> encountering is the fact that some of the
Read 'An Introduction to R', found on the R web site. It is short and worth
reading in its entirety. Section 2.5 on missing values answers your question:
2.5 Missing values
In some cases the components of a vector may not be completely known. When an
element or value is "not available" or
You really really need to do your homework.
?NA
-- Bert
On Fri, Jun 21, 2013 at 8:59 AM, Yasin Gocgun wrote:
> Hi,*
>
> *I would like to set the entries of an array that appear as "NA" to a
> certain integer, but did not find a way to do so. The related code is given
> below:
>
> A <- array(-10
<>
what kind of error? how to solve it?
[[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.
There exists an option corr.bias for making regression analyses with the
R package randomForest. The manual warns “Experimental. Use at your own
risk.” What does corr.bias precisely do and why can it be dangerous to
use it? As far as I know the bias correction is based on a linear
regression be
You should tell us where hdr.den() comes from as it is not one of
the default packages and there are thousands of packages. Some
searching suggests you are using package hdrcde.
Looking at the source code for hdr.den() it appears that the package
author used the default palette in selecting colors
Hi,*
*I would like to set the entries of an array that appear as "NA" to a
certain integer, but did not find a way to do so. The related code is given
below:
A <- array(-1000 , dim=c(100) );
for(i in 1: 100 ){
A[i] =B[i,57] - B[i,5];
}
A [ which( A == "NA") ] <- 900;
The "NA"s still appe
Rui,
package is hdrcde, and thanks for the input,
Andras
--- On Fri, 6/21/13, Rui Barradas wrote:
> From: Rui Barradas
> Subject: Re: [R] hdr.den plot colors help
> To: "Andras Farkas"
> Cc: r-help@r-project.org
> Date: Friday, June 21, 2013, 11:51 AM
> Hello,
>
> It's not your first post,
Hello,
It's not your first post, you should know to say from which package is
the function you're using.
As for the question the answer is no, those values are hard coded in the
source. What you can do is to download the source and change file hdr.den.R
Hope this helps,
Rui Barradas
Em 21
Yes I have already changed that as this was only a small vector, my real
vector is quite big..
thanks a ton.
Bets wishes,
Mitra
On 21 June 2013 23:43, Clint Bowman wrote:
> I suspect the OP may want
>
> rep(1:length(unique(x)), rle(x[order(x)])$lengths)[**order(order(x))]
>
> to allow for vari
I suspect the OP may want
rep(1:length(unique(x)), rle(x[order(x)])$lengths)[order(order(x))]
to allow for variable numbers of unique values.
Clint
Clint BowmanINTERNET: cl...@ecy.wa.gov
Air Quality Modeler INTERNET: cl...@math.utah.edu
Department of
Thanks a lot Rui,
It helps :)
I can modify this according to my need.. I didn't know about rle function.
Thanks again,
Mitra
On 21 June 2013 23:32, Rui Barradas wrote:
> Hello,
>
> I'm not sure I understand. You want to attribute a color number 1:7 to
> each element of your vector? Maybe the fo
Your "rate" variable is not a single value, but a vector of 200!.
I'm surprised you get anything at all. You have to pick one rate
value to use to create the values for the x-axis. For example run
> rate <- 1/4
And then run the qqmath() function.
Or, change the rate variable in the ggamma() func
Hello,
I'm not sure I understand. You want to attribute a color number 1:7 to
each element of your vector? Maybe the following will do.
x <- scan(text = "
43 43 43 43 0 39 13 39 50 39 39 23 23 32 0 13 32 23 32 23 0 13 13 0
")
cols <- rep(1:7, rle(x[order(x)])$lengths)[order(order(x))]
pl
Hello R experts,
I want to define desired numbers to a vector based on the present
numbers. Can anybody please help me?
Obviously I found worst ways to do it, but I believe there must be any
better way.
I have vector as
> X
[1] 43 43 43 43 0 39 13 39 50 39 39 23 23 32 0 13 32 23 32 23 0 13
Dear All,
I have the following code:
a <-rnorm(5000,10,2)
hdr.den(a,prob = c(25, 50, 75,95))
is there a way to change the colors of the horizontal bars representing the
given highest density regions from the current white, green, red and black to
colors of the shades of grey going from lighter
Hello,
I am trying to superpose on a single panel two QQ plots with the lattice
qqmath function.
Here is a reproducible example of the problem I am facing:
# Generate data
shape = 8
rate = c(rep(1/4,100),rep(1/3,100))
x = rgamma(200,shape,rate)
groups = gl(2,100,200,labels=LETTERS[1:2])
# Pl
I think you have to assume there could be other coding errors as
well, such as misspellings or abbreviating Street as St. You
probably will need to use sub() to correct A2 and possibly A1 before
trying to merge. To figure where the problems are, you might try
something like this. The last command l
Hi,
Try:
dat1<- read.csv("mydata.csv",sep="\t",row.names=1) #sep could be different in
your case.
dat1
# A B C D E
#Jhon 1 3 6 9 8
#Kelly 2 3 4 6 9
#Audrey 3 5 6 9 7
A.K.
Hi,
I am calling my csv data with 37 columns and 16 rows into R. However I really
need to keep coloumns name
Dear R users,
I have a seasonal time series of period 4 (my_ts).
I would like to apply to my_ts a Seasonal ARMA(2,0) process (only the seasonal
part), with Seasonal AR coefficients respectively 0.54 and 0.5.
I tried to use the arima command from the stats package, but this code is not
correct:
a
Hi all,
So this code works for three of the six seamounts. The problem I am
encountering is the fact that some of the seamounts doesn't have all the
depth categories (i.e missing the deep category) and thus gives me this
error "Error in plot.window(...) : need finite 'ylim' values."
How can I fix
Hello,
Just to add to the last post, you can try to get that p-value using what
that output gives you:
df1 <- 1
df2 <- 356
F_stat <- 10.948
pf(F_stat, df1, df2, lower.tail = FALSE) # same value
Hope this helps,
Rui Barradas
Em 21-06-2013 12:11, Rui Barradas escreveu:
Hello,
The Pr(>F) i
Hello,
The Pr(>F) is a regular p-value, with the usual meaning.
Rui Barradas
Em 21-06-2013 11:05, pieter escreveu:
k.
Thanks.
Sorry i'm not the biggest statician, but I still find the output write
out not the clearest of all. Anyway, that is probably a personal
problem. :-)
Just one question.
Dear list,
I have already posted to sci.stat.consult, however, I was hoping that
there might be some experience in the R community on the following question:
I am wondering whether anyone has experience with variants of (parametric) Gini
coefficients (or alternatives) which allow for reflect
> I don't really get the output of the performed Levene's test.
One answer would be to read John Fox's rather good book, for which car is the
corresponding R package.
Another is to look up leven'es test; a reasonably authoritative web reference
for the calculation is
http://www.itl.nist.gov
Once you have this matrix:
mx <- cbind(x1, x2, x3, x4)
> > x1 x2 x3 x4
> > [1,] 22 14 24 10
> > [2,] 27 12 1 2
> > [3,] 23 13 6 2
> > [4,] 27 6 21 2
> > [5,] 14 29 12 19
> > [6,] 5 14 11 29
> > [7,] 10 10 27 25
> > [8,] 16 29 28 10
> > [9,] 22 22 19 7
> > [10,] 30 20 5 22
> > [11,] 2
Hello everybody
I have this problem: I need to match an addresses database F1 with the
information contained in a toponymic database F2.
The format of F1 is given by three columns and 800 rows, with the
columns being:
A1. Street/Road/Avenue
A2. Name
A3. Number
Consider for instance Avenue J. Ke
I received the same error message. With me, the problem was that I had the
cells in my CSV file in Excel formatted as percentages. R read them as
factors and not as numeric values. After changing the cell format from
percentage to general in Excel, the problem was solved.
--
View this message in
hi
I'm trying to understand (a little) the code behind summary.gam, and have
the Biometrika article referred to in the Help. But am stuck early on.
In the code, starting at line 167:
if (est.disp)
rdf <- residual.df
else rdf <- -1
res <- testStat(p, Xt, V, df[i], type = p.type,
res.df = rdf)
Hey,
I don't really get the output of the performed Levene's test.
So in short.
two data sets, one with 75 elements, one with 283.
And I want to check whether they have the same variance.
I did a LeveneTest
And this is the output:
Levene's Test for Homogeneity of Variance (center = median)
On Jun 20, 2013, at 16:21 , Annette Timm wrote:
> Hi All
> I am having some trouble making a test for linearity of my passing bablok
> regression.As far as I am concerned i need to make cusum test for linearity
> to get a p-value. I can make a cusum plot, but I can not figure out how to
> test
Hi Greg
On 21/06/13 08:02, Greg Dropkin wrote:
> hi
>
> I'm trying to understand (a little) the code behind summary.gam, and have
> the Biometrika article referred to in the Help. But am stuck early on.
>
> In the code, starting at line 167:
>
> if (est.disp)
>rdf <- residual.df
> else rdf <-
Hello
My colleagues asked me to install "R" with module "shiny".
R version 3.0.1 compiled fine on Solaris 10 and is running well.
I tried to install "shiny". With the dependencies "Rcpp" should be installed
before. But the compile step did fail. See below.
The initial error message is "Err
52 matches
Mail list logo
