Dohso simple when you know how! That worked straightaway! Thank you
Martin!
--
View this message in context:
http://r.789695.n4.nabble.com/Setting-Derived-Class-Slots-tp4671683p4671737.html
Sent from the R help mailing list archive at Nabble.com.
___
Hi Simon.many thanks to you for your explanation. - I'll check out
setReplaceMethod. Point taken on nomenclature too...
--
View this message in context:
http://r.789695.n4.nabble.com/Setting-Derived-Class-Slots-tp4671683p4671738.html
Sent from the R help mailing list archive at Nabble.com.
Hi r-users,
I would like to use gev and my data (annual rainfall ) is as follows:
> head(dat,20) ABCDEFGHI J
1 45.1 41.5 58.5 50.1 46.0 49.1 37.7 49.1 59.8 54.0
2 50.3 39.8 49.4 56.4 49.4 48.8 42.1 49.8 49.4 58.3
3 41.7 39.3 44.6 39.1 35.7 41.5 40.8 4
Dear List,
This is Elaine.
I am using beta.part to calculate the beta diversity index.
The function of beta.part can generate three kinds of beta diversity indice
at one time.
However, it is Simpson index that I want.
The calculation result is composed of three rows and N/A columns (dist
class).
Hello,
You should say which package is gev() coming from. I believe it's from
evir, and what follows assumes that.
The problem is that 'dat' is a data.frame, a special type of list. It is
a list of vectors and gev is expecting just one vector. You could try
instead, to fit generalized extreme
Hello,
you may just first change the boxplot.stats directly:
boxplot.stats <- function (x, coef = NULL, do.conf = TRUE, do.out =
TRUE)
{
nna <- !is.na(x)
n <- sum(nna)
stats <- quantile(x, c(.05,.25,.5,.75,.95), na.rm = TRUE)
iqr <- diff(stats[c(2, 4)])
Hi R-Help,
When we tried to install R 2.8.1 in our AIX system, we faced some error,
please help.
OS level: AIX 7.1 64 bit
R 2.8.1
We had a guide about How to install R 2.8.1 in AIX from R source code. Here
is what I did:
1. Download the R source file, for example R-2.8.1.tar.gz , from
http://
Hi, I'm making a simple population pyramid using two back-to-back bar plots for
18 different age groups with totals for males and females.
Yesterday I achieved a fairly serviceable plot using the following.
library(ggplot2)
library(reshape2)
library(plyr)
#make sample data
df<-data.frame(ag=c(1:
Hello,
It's better if you Cc the list, the odds of getting more and better
answers are greater.
Have you tried with different block values?
Rui Barradas
Em 17-07-2013 13:33, Roslina Zakaria escreveu:
Hi Rui,
Yes, I am using gev from 'evir'
install.packages("evir")
library(evir)
I have tried
Solved it!! ;)
The problem was that the test-data contained factor-levels the
training-data didn't.
So when trying to fit this new factor-levels to a model which didn't have
this levels, the error occurred.
When excluding the factor-levels not used when fitting the model or when
taking care all l
- Forwarded Message -
From: Rui Barradas
To: Roslina Zakaria
Cc: 'r-help'
Sent: Wednesday, July 17, 2013 8:49 PM
Subject: Re: [R] error message in gev
Hello,
It's better if you Cc the list, the odds of getting more and better
answers are greater.
Have you tried with differen
Hello Arnaud,
Thank you for your pointer. However I need to more clarification.
I want to control the max. and min. values for the x-axis, as well as
number of vertical gridlines to be displayed. I tried the following:
MyData <- data.frame(Names1 = paste("XXX", 1:150), Values1 = 1:150 +
10, Val
Hi Rui,
Â
I have tried different block values and below are the error message I got:
 > out <- gev(dat, "A") Error in as.POSIXlt.default(attributes(data)$times) :
do not know how to convert 'attributes(data)$times' to class âPOSIXltâ
Â
out.list <- lapply(dat, gev, block = 10) Error in optim
Replied to on the R-sig-ecology list, with a request not to cross-post.
Briefly, the output is a list with named components, so can be
extracted as usual.
Sarah
On Wed, Jul 17, 2013 at 5:50 AM, Elaine Kuo wrote:
> Dear List,
>
> This is Elaine.
> I am using beta.part to calculate the beta diver
Hi
You can do a rotation and use gvisColumnChart instead gvisBarChart
plot(gvisColumnChart(MyData, xvar="Names1", yvar=c("Values1",
"Values2"),options=list(width=2500,height=1000)))
Michel
Le 17/07/2013 15:57, Christofer Bogaso a écrit :
Hello Arnaud,
Thank you for your pointer. However I nee
Hi,
YOu could try.
?split()
split(ats,ats$Project_NBR)
You also mentioned about two columns.
split(ats,list(ats$col1, ats$col2))
You should have provided an example dataset using ?dput() ( dput(head(data,10))
) for testing.
Also,
gsub("^-[^-]*-","","-005-190")
#[1] "190"
A.K.
Problem: I hav
Hello,
And what if you try it without a block value?
out.list <- lapply(dat, gev)
Rui Barradas
Em 17-07-2013 14:53, Roslina Zakaria escreveu:
Hi Rui,
I have tried different block values and below are the error message I got:
out <- gev(dat, "A")
Error in as.POSIXlt.default(attributes(data)
Thanks, I've now got the plot I was after. No idea what the stacking warning
is about, but it does not seem to compromise my output.
Gavin.
--
View this message in context:
http://r.789695.n4.nabble.com/Some-problems-with-back-to-back-bar-plots-in-ggplot2-tp4671744p4671765.html
Sent from the
> Hi, I'm making a simple population pyramid using two
> back-to-back bar plots for 18 different age groups with
> totals for males and females.
I can run both variants without error, though both throw the same warning
regarding stacking beein 'not well defined ..'
However, the axis scaling in
Help with outdated versions of R, and particularly with nonstandard patch
files, are not really on topic here. Posting in HTML format is not advised (see
posting guide) since it actually makes it harder to read R code, so text
coloring doesn't come through.
It appears to me that your patch fail
Hi Arnaud,
Thanks for your answer.
However I prefer to have Horizontal bar chart, because length of the
strings for x-axis are quite large, therefore it would be better to
put them vertically one-below-another instead horizontally.
Therefore I would really appreciate if someone points me how to
HI Helpers,
How could we use R and R packages licensed under GPL into commercial products?
Is it allowed to load a library and get the results from it and using the
results for commercial use? Thank you so much!
Regards,
Yan
[[alternative HTML version deleted]]
___
Many thanks, this works beautifully.
Louise
On Wed, Jul 17, 2013 at 7:42 AM, Paul Murrell wrote:
> Hi
>
> There are a couple of problems:
>
> 1.
> Your 'outline' is much bigger than it needs to be. For example, the
> following produces just Australia ...
>
> outline <- map("worldHires", region
Dear users,
I want to compute r-squared values from a glm regression using a gamma
distribution and an "identity" link-function, but find no such thing when
using the summary() or names() function. My next guess was to calculate it
by "hand", i.e.
r2 = (sum((estimate - xbar)^2) /sum((x-xbar)^2))
On Jul 17, 2013, at 10:56 AM, "Li, Yan" wrote:
> HI Helpers,
>
> How could we use R and R packages licensed under GPL into commercial
> products? Is it allowed to load a library and get the results from it and
> using the results for commercial use? Thank you so much!
>
> Regards,
> Yan
For
Hi Arun
I have two questions always about the question of symplify a dataframe
I would like
1) to transform the vector x1 into the vector y1
x1 <- c(1,1,1,-1000, 1,-1000, 1,1,1,1,1,1,-1000)
y1 <- c(1,1,1,1,2,2, 3,3,3,3,3,3,3)
2) to transform the vectors Deb
Hi,
No problem.
You could try:
printer = file("out.txt","w")
writeLines("This is line.",con=printer,sep=" ")
writeLines("The same line.",con=printer)
close(printer)
#or
cat(sprintf("This is line %d. ",1),file="out.txt",append=TRUE)
cat("The same line.",file="out.txt",append=TRUE)
A.K.
Thank
What is the function to do the truncation to a certain decimal digit of a
number. And the function approximation?
Thanks.
Francesco M
[[alternative HTML version deleted]]
__
HI,
One problem with using ?subst() would be it depends upon the number of digits,
characters etc.
For eg.
substring("-005-190",6)
#[1] "190"
substring("-0057-190",6)
#[1] "-190"
#whereas
gsub("^-[^-]*-","","-0057-190")
#[1] "190"
Probably, your dataset doesn't have that sort of problem.
d
Hello,
As for question (1), try the following.
y2 <- cumsum(c(TRUE, diff(x1) > 0))
identical(as.integer(y1), y2) # y1 is of class "numeric"
As for question (2) I'm not understanding it.
Hope this helps,
Rui Barradas
Em 17-07-2013 18:21, Arnaud Michel escreveu:
Hi Arun
I have two questio
Dear R Users---
I spent some time implementing the AS75 WLS regression algorithm in C
in order to learn the R interface to C. I did not find an easy, brief
C+R example of an algorithm that used persistent storage across calls.
they probably exist, but I could not find a nice brief template for
m
Hi R user,
I am wondering how we can open the data file which have .WTG extension?
Thanks for your help
cheers,
KG
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.
I am wondering what .WTG files have to do with R.
Care to elaborate? (Suggestion: Read and follow the posting guide)
Cheers,
Bert
On Wed, Jul 17, 2013 at 11:48 AM, Kristi Glover
wrote:
> Hi R user,
> I am wondering how we can open the data file which have .WTG extension?
> Thanks for your help
On Jul 17, 2013, at 1:48 PM, Kristi Glover wrote:
> Hi R user,
> I am wondering how we can open the data file which have .WTG extension?
> Thanks for your help
> cheers,
> KG
A quick search led me to http://wasptechnical.dk/Forum/viewtopic.php?id=433,
which suggests that the .WTG extension ma
Hi Bert,
This is the data related to climate. I can open it in Excel but the problem is
I do have 190 files .
thanks
> Date: Wed, 17 Jul 2013 12:00:58 -0700
> Subject: Re: [R] How to open .WTG file Extension in R?
> From: gunter.ber...@gene.com
> To: kristi.glo...@hotmail.com
> CC: r-help@r-p
That reduces the number of possibilities for the file structure to
ones that R can probably read, but we need to know more. Try opening
a sample file in a text editor to see if it is plain text. If so,
you can probably use read.table() or read.csv() to get the file into
R. If not, you will need one
Hello,
As for truncation, you can write a one line function:
dec_trunc <- function(x, digits = 0) trunc(x * 10^digits)/10^digits
As for approximation, maybe you're looking for ?signif.
Hope this helps,
Rui Barradas
Em 17-07-2013 18:41, Francesco Miranda escreveu:
What is the function to do
Hello,
So it's a XML file, like Marc said. See package XML.
install.packages('XML', dependencies = TRUE)
Hope this helps,
Rui Barradas
Em 17-07-2013 20:04, Kristi Glover escreveu:
Hi Bert,
This is the data related to climate. I can open it in Excel but the problem is
I do have 190 files .
Thank you for the question (1)
Sorry for the imprecision for the question (2) :
Suppose the date frame df
df1 <- data.frame(
Debut =c ( "24/01/1995", "01/05/1997" ,"31/12/1997", "02/02/1995"
,"28/02/1995"
,"01/03/1995", "13/03/1995", "01/01/1996", "31/01/1996") ,
Fin = c ( "30/04/1997", "30/12/
How can id do this efficiently in R ?
1 0 0
0 2 0
0 0 3
rotate right
0 0 1
0 2 0
3 0 0
rotate left
0 0 3
0 2 0
1 0 0
What I want to do is described here:
http://stackoverflow.com/questions/13049575/r-min-max-and-mean-of-off-diagonal-elements-in-a-matrix
but I want do to it for all off-anti-dia
Hi,
You could try:
df1[,1:2]<-lapply(df1[,1:2],as.character)
df2New<- data.frame(Deb=unique(with(df1,ave(Debut,INDX,FUN=function(x)
head(x,1,Fin=unique(with(df1,ave(Fin,INDX,FUN=function(x) tail(x,1)
identical(df2New,df2)
#[1] TRUE
A.K.
- Original Message -
From: Arnaud Michel
Hi,
I followed the code given in the link
http://www.forecastingfinancialrisk.com/3.html
EWMA<-matrix(nrow=T,ncol=5)
S<-cov(returns)
dim(S)
[1] 5 5
"EWMA[1,] = c(S)[c(1,4,2)] # extract the variances and covariance"
*here we will get a vector but we are equating with a matrix. hence
the foll
#or
library(plyr)
res<-ddply(df1,.(INDX),summarize,Debut=head(Debut,1),Fin=tail(Fin,1))
res$INDX<-factor(res$INDX,levels=unique(df1$INDX))
res[order(res$INDX),-1]
# Debut Fin
#3 24/01/1995 31/12/1997
#4 02/02/1995 12/03/1995
#1 13/03/1995 30/06/1995
#2 01/01/1996 31/01/1996
A.K.
---
hi,
Could anyone help me in solving the following error:
I have 5 stocks returns data (returns)
EWMA = matrix(nrow=T,ncol=5) # create a matrix to hold the
covariance matrix for each t
lambda = 0.94
S<-cov(returns) # initial (t=1)
covarian
# The following dataframes are the result of two analyses performed on
the same set of numeric data.
# The first analysis involved calculations that did not include zero values:
StatsUTAH = data.frame(MWtotaleesDue =
c(8.428571,2.496256,7,6.604472,1,17,3.593998,4.834573,12.02257),
On Jul 17, 2013, at 3:04 PM, Witold E Wolski wrote:
> How can id do this efficiently in R ?
>
> 1 0 0
> 0 2 0
> 0 0 3
>
>
> rotate right
> 0 0 1
> 0 2 0
> 3 0 0
> rotate left
> 0 0 3
> 0 2 0
> 1 0 0
>
> What I want to do is described here:
> http://stackoverflow.com/questions/13049575/r-min-m
Thanks Arun and Rui for your helps
Michel
Le 17/07/2013 22:20, arun a écrit :
#or
library(plyr)
res<-ddply(df1,.(INDX),summarize,Debut=head(Debut,1),Fin=tail(Fin,1))
res$INDX<-factor(res$INDX,levels=unique(df1$INDX))
res[order(res$INDX),-1]
# DebutFin
#3 24/01/1995 31/12/1997
#4 02/
This is most likely a silly question.
First I run the following:
require(car)
mod.ok <- lm(cbind(pre.1, pre.2, pre.3, pre.4, pre.5, post.1, post.2,
post.3, post.4, post.5,
fup.1, fup.2, fup.3, fup.4, fup.5) ~
treatment*gender, data=OBrienKaiser)
phase <- factor(rep(c("pr
I would like to "cut" a vector into groups of equal nr of elements.
looking for a function on the lines of cut but where I can specify
the size of the groups instead of the nr of groups.
--
Witold Eryk Wolski
__
R-help@r-project.org mailing list
http
In addition to what Petr said,
That message most likely means that you have some missing data where
you're not allowed to have missing data in your Datos.
Look and see if the mctp() function has an argument that controls what it
does with missing data.
-Don
--
Don MacQueen
Lawrence Livermore
Looking for robust smoothing methods available in R... a robust
version of lowess or smooth.spline.
Was searching for "median and smoothing" but the search engine isn't
calibrated .
--
Witold Eryk Wolski
__
R-help@r-project.org mailing list
https://stat
If you are happy with the character strings printed out then you can use
the capture.output function.
If you want the numbers (without needing to convert) then look at the print
method for the class of object that you are working with and see what the
code is. It is probably calling another funct
You could use the quantile function to choose the cut points and pass that
to cut.
Or you could sort the vector and just take the first n elements as the 1st
group, etc.
On Wed, Jul 17, 2013 at 3:43 PM, Witold E Wolski wrote:
> I would like to "cut" a vector into groups of equal nr of elements
HI,
Not sure whether this is what you wanted.
vec1<- 1:7
fun1<- function(x,nr) {((x-1)%/%nr)+1}
fun1(vec1,2)
#[1] 1 1 2 2 3 3 4
fun1(vec1,3)
#[1] 1 1 1 2 2 2 3
split(vec1,fun1(vec1,2))
A.K.
- Original Message -
From: Witold E Wolski
To: r-help@r-project.org
Cc:
Sent: Wednesday,
Hello,
Maybe ?stats::smooth
Hope this helps,
Rui Barradas
Em 17-07-2013 22:47, Witold E Wolski escreveu:
Looking for robust smoothing methods available in R... a robust
version of lowess or smooth.spline.
Was searching for "median and smoothing" but the search engine isn't
calibrated .
--
Wi
Dear List,
I'm running simulations using the glmnet package. I need to use an
'automated' method for model selection at each iteration of the simulation.
The cv.glmnet function in the same package is handy for that purpose.
However, in my simulation I have p >> N, and in some cases the selected
mo
Sorry, there was a mistake:
fun1 should be:
fun1<- function(x,nr) {((seq_along(x)-1)%/%nr)+1}
vec3<- c(4,5,7,9,8,5)
fun1(vec3,2)
#[1] 1 1 2 2 3 3
split(vec3,fun1(vec3,2))
A.K.
- Original Message -
From: arun
To: Witold E Wolski
Cc: R help
Sent: Wednesday, July 17, 2013 6:04 PM
Su
?qcut
On Jul 17, 2013 5:45 PM, "Witold E Wolski" wrote:
> I would like to "cut" a vector into groups of equal nr of elements.
> looking for a function on the lines of cut but where I can specify
> the size of the groups instead of the nr of groups.
>
>
>
>
> --
> Witold Eryk Wolski
>
> __
On Jul 17, 2013, at 4:43 PM, Witold E Wolski wrote:
> I would like to "cut" a vector into groups of equal nr of elements.
> looking for a function on the lines of cut but where I can specify
> the size of the groups instead of the nr of groups.
In addition to the other options, if the 'breaks'
Tena koe
Without reading your request in detail, I will suggest you look at ?merge. It
is often the answer when 'combine' is in the question.
Peter Alspach
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of bcrombie
Sent: Thursday,
Hello,
1) In the provided example, you have 2 stock returns and
> EWMA = matrix(nrow=T, ncol=3)
See the number of columns (2+1)
Please modify the number of columns according to the number of stock
returns.
2) As you have 5 stock returns, the following cannot work
> EWMA[1,] = c(S)[c(1,4,2)]
Y
Hello,
> require(sos)
> findFn("median and smoothing")
Regards,
Pascal
2013/7/18 Rui Barradas
> Hello,
>
> Maybe ?stats::smooth
>
> Hope this helps,
>
> Rui Barradas
>
> Em 17-07-2013 22:47, Witold E Wolski escreveu:
>
> Looking for robust smoothing methods available in R... a robust
>> ver
Thanks a lot for the pointer!
After I downloaded the source code and saw the innards of the print
function, I know what to do now.
Thanks again,
Gang
On Wed, Jul 17, 2013 at 6:00 PM, Greg Snow <538...@gmail.com> wrote:
> If you are happy with the character strings printed out then you can use t
Another option is
library(locfit)
?locfit
There is also a web page from the help which gives a link to Clive
Loaders book
Regards
Duncan
Duncan Mackay
Department of Agronomy and Soil Science
University of New England
Armidale NSW 2351
Email: home: mac...@northnet.com.au
At 08:23 18/07/2013
On Jul 17, 2013, at 5:26 PM, Axel Urbiz wrote:
Dear List,
I'm running simulations using the glmnet package. I need to use an
'automated' method for model selection at each iteration of the
simulation.
The cv.glmnet function in the same package is handy for that purpose.
However, in my simul
On 18/07/13 02:36, Chris89 wrote:
Dear users,
I want to compute r-squared values from a glm regression using a gamma
distribution and an "identity" link-function, but find no such thing when
using the summary() or names() function. My next guess was to calculate it
by "hand", i.e.
r2 = (sum((es
Hi Rui,
Â
Yes, it works beatifully this time.
Â
Thank you so much for your help.
From: Rui Barradas
Cc: R help forum
Sent: Wednesday, July 17, 2013 11:28 PM
Subject: Re: [R] error message in gev
Hello,
And what if you try it without a block value?
o
Hi,
Not sure if this is what you wanted:
#If columns are arranged in the same order in both data.frames.
lst1<-lapply(seq_len(ncol(StatsUTAH)),function(i)
{x1<-cbind(StatsUTAH[,i],sStatsUTAH[,i]);row.names(x1)<-row.names(StatsUTAH);colnames(x1)<-c("zeroNO","zeroYES");x1})
names(lst1)<- colnames(
68 matches
Mail list logo
