Hi
r-help-boun...@r-project.org napsal dne 04.06.2009 16:11:29:
> I want to make a log-log plot with a regression line, but I can't figure
> out what I'm doing wrong. What I'm trying is:
>
> plot(mass,area, log="xy", pch=as.numeric(food))
> abline(lm(mass~area))
abline(lm(area~mass))
>
> or
I'm all done now. The "max2" version below is what I went with in the
end for my proposed change to caret::nearZeroVar (which used the "sort"
method). Max Kuhn will make it available on CRAN soon. It speeds up
that routine by a factor 2-5 on my test cases and uses much less
memory. For what
Dear R
expert
I am a
student and I am currently conducting a research project on the Modeling Loss
Index Triggers to price Cat Bonds: Application of the risk of hurricanes in
USA.
I need to
solve with R (especially with EM algorithm) this specific problem below. CRAN
Package archive doesn't se
On Thu, 2009-06-04 at 15:21 +0100, Daniel Mail wrote:
> Hello,
>
> I need to know how to import ARIMA coefficients. I already determined
> the coefficients of the model with other software, but now i need to
> do the forecast in R.
>
>
> > fit <- arima(x, order=c(1, 0, 1), fixed=c(-.172295,
venkata kirankumar wrote:
Hi all,
I am new to R-project my problem is I tried to get subset from two different
tables its giving error
but if i m tring for geting results from one table its working
actually i have to take values from two tables with applying different
conditions on two tables
A mix of Petr and William's answers will sort you out. You had two problems:
the base 10/e problem and the axes swapping problem.
?dput
I didn't know about that one either - very useful!
andydol...@gmail.com
2009/6/5 Petr PIKAL
Hi
>
> r-help-boun...@r-project.org napsal dne 04.06.2009 16:1
Hi,
I am trying to install R 2.9.0 on AIX 5.3. During the 'make' run I get
the following error:
Error: Invalid DESCRIPTION file
Invalid Priority field.
Packages with priorities 'base' or 'recommended' or 'defunct-base' must
already be known to R.
See the information on DESCRIPTIO
Duncan Temple Lang a écrit :
Olivier Cailloux wrote:
Dear list,
I am trying to use the SSOAP package to access a very simple SOAP
service, included as a demonstrator in the axis2 Apache package:
"Version" (takes no input parameters and should return a string with
the version number). The w
A new package is now available via CRAN, called metafor.
The metafor package consists of a collection of functions for conducting
meta-analyses in R. Fixed- and random-effects models (with and without
moderators) can be fitted via the general linear (mixed-effects) model. For 2x2
table data, th
Irina Foss wrote:
> Dear,
>
> I am using "cwt "function from Rwave package to perform Morlet wavelet
> analysis.
> d1<-ts(d1deltat=1,start=c(1960,1))
> library(Rwave)
> D1-cwt(d1,noctave=8,nvoice=1)
>
> After I did Morlet wavelet analysis by using function "wavelet" from dplR
> package.
>
>
Hi,
While doing the ADF test in R using the following command I am getting the
error and the result..">
x.ct=ur.df(rev$REVENUE,start=1,end=length(rev$REVENUE),frequency=1)
Error in ur.df(rev$REVENUE, start = 1, end = length(rev$REVENUE), frequency
= 1) :
unused argument(s) (start = 1, end = 4, fr
Thanks guys, that's information if very much appreciated. Both
solution are better than mine which was to use capture.output after
using the cat function.
re: collapse argument in paste() -- I've always wondered what that
argument was for, i thought it was basically doing what the sep
argument doe
Hello,
I'm just looking for an easy way to find the positions of a complete
sequence in a bigger vector.
For example :
c("a","z","e") in c("a","z","e","r","t","a","z","a","z","e","c")
and the result should be
1 8
that is the positions of the beginning of the complete sequence.
I tried with %in%,
use gregexpr and paste
> aze <- paste(c("a", "z", "e"), collapse = "")
> sequence <- paste(c("a","z","e","r","t","a","z","a","z","e","c"), collapse =
> "")
> gregexpr(aze, sequence, fixed = TRUE)
[[1]]
[1] 1 8
attr(,"match.length")
[1] 3 3
HTH,
--sundar
On Fri, Jun 5, 2009 at 6:22 AM, Ptit_Ble
Dear R users,
I had somehow expected that read.table() would treat the column specified by
the row.names argument as of class character. That seems to be the only
sensible class allowed for a column containing row names. However, that does
not seem to be the case, as the following example shows:
Dear all,
is there any possibility to react on a process termination signal in R?
In a cluster (using RMPI and snow) if the main R process is aborted due
to runtime limitations I want to properly shutdown all spawned node
processes, e.g. via mpi.close.Rslaves.
Any hints are welcome,
Christi
Look into the R packages for missing data by Joe Schafer. These
missing data routines use EM. See here for starters. He also wrote a
book:
http://www.stat.psu.edu/~jls/misoftwa.html
Janet
__
R-help@r-project.org mailing list
https://stat.ethz.
Dear all,
I am using the RWeka package to append several arff files. Although it
works the resulting arff files always have "@relation R_data_frame",
and I have to change this manually to my desired relation name. Can
the package accomplish this for me instead?
Thank you,
Wil Koetsier
_
Hello,
Does someone know what are theta.mod coefficients in the output of mda()
function in the mda package?
Thanks for your help
Amélie Vaniscotte,
Department of Chrono-environment, UMR UFC/CNRS 6249 aff. INRA
University of Franche-Comté
amelie.vanisco...@univ-fcomte.fr
__
Hi every one
I perform a simple linear regression
lm(a b + c + d , data = data1)
How to say to R to perform and print the regression with restricting the
coefficient
of the variable c to be equal to 0.1. In the model print, I want to
show the p-values of all my coefficients.
Thank you in adv
On Jun 5, 2009, at 11:10 AM, Axel Leroix wrote:
Hi every one
I perform a simple linear regression
lm(a b + c + d , data = data1)
How to say to R to perform and print the regression with restricting
the coefficient
of the variable c to be equal to 0.1.
?lm
Examine material on specifyin
koyel chakrabartti wrote:
> While doing the ADF test in R using the following command I am getting the
> error and the result..">
> x.ct=ur.df(rev$REVENUE,start=1,end=length(rev$REVENUE),frequency=1)
> Error in ur.df(rev$REVENUE, start = 1, end = length(rev$REVENUE), frequency
> = 1) :
> unused a
Dear UseRs,
A simple class inheritance example:
setClass("test",representation(a="numeric"))
setMethod("initialize","test",
function(.Object,x,...){
print("Initialization!!!")
callNextMethod(.Object,a=x,...)
})
new("test",x=23)
[1] "Initializa
Hi Vitalie --
Vitalie S. wrote:
> Dear UseRs,
>
> A simple class inheritance example:
>
>> setClass("test",representation(a="numeric"))
>> setMethod("initialize","test",
> function(.Object,x,...){
> print("Initialization!!!")
> callNextMethod(.Object,a=x,...
Markus Loecher gmail.com> writes:
> [...]
> > x
> ID X1 X2
> 1 010007787048271871 1 4
> 2 1007109516820319 2 3
> 3 10094843652996959 3 2
> 4 010145176274075487 4 1
>
> The first column was not read in as a string, which mangled the IDs.
> I could use colClasses expli
Le vendredi 05 juin 2009 à 14:03 +0200, Viechtbauer Wolfgang (STAT) a
écrit :
> A new package is now available via CRAN, called metafor.
>
> The metafor package consists of a collection of functions for
> conducting meta-analyses in R. Fixed- and random-effects models (with
> and without moderator
Here is one approach (this could be wrapped into a function if you are applying
it repeatedly):
> one <- c("a","z","e")
> two <- c("a","z","e","r","t","a","z","a","z","e","c")
>
> which( sapply( 1:(length(two)-length(one)+1),
+ function(i) isTRUE( all.equal( one, two[ i + 0:(length(one)-1) ] )
I have a large dataset that contain duplicate records. How do I identify and
remove duplicate records?
Chris Anderson
707.315.8486
www.sassydeals4u.com
Free info for small business owners. Click here to find great products geared
for
Dear R-helpers,
I have been trying to figure out how to plot a graph with an axis label
consisting of a mixture of Latin, Greek and subscript characters.
Specifically, I need to write A[beta]{1-42}, where A is Latin script A,
[beta] is Greek lower case beta and {1-42} is subscript '1-42'.
I can us
Try this:
d <- data.frame(a = c(1, 1, 2, 3), b = c(10, 10, 9, 8))
unique(d)
On Fri, Jun 5, 2009 at 1:38 PM, Chris Anderson wrote:
> I have a large dataset that contain duplicate records. How do I identify
> and remove duplicate records?
>
>
> Chris Anderson
> 707.315.8486
> www.sassydeals4u.co
Chris,
How large is large? How may columns?
"Duplicate" across all columns of just some?
Henrique gave you simple R answer. Perhaps doing in SQL is more efficient?
eg
SELECT DISTINCT
FROM ;
HTH,
Jim Porzak
TGN.com
San Francisco, CA
www.linkedin.com/in/jimporzak
use R! Group SF
I am woking on Windows and using read.xls to read Excel.xls file.
But it needs to install perl on my machine first, and then give the
following until it can work:
read.xls(file,sheet=i,perl="E:\\Perl\\bin\\bin\\perl.exe");
I am wondering if there is other way to read Excel.xls file without needi
I think I would use Sundar's solution but here are a couple
more just in case:
> x <- c("a", "z", "e")
> y <- c("a", "z", "e", "r", "t", "a", "z", "a", "z", "e", "c")
>
> #1
> which(apply(embed(y, 3), 1, identical, rev(x)))
[1] 1 8
> #2
> library(zoo)
> which(rollapply(zoo(y), 3, identical, x, al
Check out:
http://wiki.r-project.org/rwiki/doku.php?id=tips:data-io:ms_windows
On Fri, Jun 5, 2009 at 1:37 PM, jlfmssm wrote:
> I am woking on Windows and using read.xls to read Excel.xls file.
> But it needs to install perl on my machine first, and then give the
> following until it can work:
>
Markus Loecher-4 wrote:
>
> Dear R users,
> I had somehow expected that read.table() would treat the column specified
> by
> the row.names argument as of class character. That seems to be the only
> sensible class allowed for a column containing row names. However, that
> does
> not seem to be
Chris Anderson wrote:
I have a large dataset that contain duplicate records. How do I identify and
remove duplicate records?
Here's one way:
> aq <- airquality[sample(NROW(airquality), replace=TRUE),]
> any(duplicated(aq))
[1] TRUE
> which(duplicated(aq))
[1] 2 15 34 44 45 47 49 5
Jonathan Williams wrote:
Dear R-helpers,
I have been trying to figure out how to plot a graph with an axis label
consisting of a mixture of Latin, Greek and subscript characters.
Specifically, I need to write A[beta]{1-42}, where A is Latin script A,
[beta] is Greek lower case beta and {1-42} is
Hi
I want to create a new generic method, but I end up with an error
(evaluation nested too deeply). see the transcript below.
The function beginYear.Fun() works, but not beginYear.
I have no idea why.
Any ideas welcome,
Rainer
> setClass("fun", representation(x = "numeric"))
[1] "fun"
> new(
Anyone know how to get p-values for the t-values from the coefficients
produced in vglm?
Attached is the code and output — see comment added to output to show
where I need p-values
+ print(paste("** Using VGAM function gamma2 **"))
+ modl2<-
vglm(MidPoint~Count,gam
Hi Robert
As Brian Ripley wrote in a message to the list earlier today,
you will want a more recent version of libcurl.
RCurl can be easily made to check the availability of all the different
options to handle old versions of libcurl, but this slows down the
installation for everyone. But more
Hi there,
I have 3 numeric variables: day (e.g., 05), month (e.g., 06), year (e.g.,
2009).
I would like to create a (string) variable of the following form:
month/day/year (e.g., 06/05/2009).
I would be grateful to anyone who could point me toward a solution.
Sincerely,
Marc
=
?paste
paste(month,day,year, sep="/")
On Jun 5, 2009, at 4:56 PM, Marc Belisle wrote:
Hi there,
I have 3 numeric variables: day (e.g., 05), month (e.g., 06), year
(e.g.,
2009).
I would like to create a (string) variable of the following form:
month/day/year (e.g., 06/05/2009).
I would b
Dear R users,
I was trying to find a library for the Chi-plot and Kendall plot, as described
in the paper "Everything you always wanted to know about Copula Modeling but
you were afraid to ask", but unsuccessfully. I was wondering if anyone have
already implemented these 2 functions in R?
Ma
Dear R users,
I was trying to find a library for the Chi-plot and Kendall plot,
as described in the paper "Everything you always wanted to know about
Copula Modeling but you were afraid to ask", but unsuccessfully. I
was wondering if anyone have already implemented these 2 functions in R
and c
Dear UseRs,
Does anyone know when exactly the validity is checked in S4? Documentation is silent:(.
Here is a small example:
setClass("test1",representation(a="numeric"))
setMethod("initialize","test1",
function(.Object,...){
a<-runif(1) ## here slot "a" is initial
Hi the list
Strictly speaking, this is not a "R" question, but I need the
information for the
creation of a package. My question is about vocabulary: What is the
difference between
"time series", "longitudinal data" and "trajectories"?
Sincerely
Christophe
__
Olivier Cailloux wrote:
Duncan Temple Lang a écrit :
Olivier Cailloux wrote:
Dear list,
Sorry about that, my e-mail has been sent too soon by mistake (also from
an incorrect exp. address). I was planning to add further details.
So I'm running Debian Lenny. sessionInfo() gives:
R version
On Jun 5, 2009, at 3:15 PM, Steven Matthew Anderson wrote:
Anyone know how to get p-values for the t-values from the
coefficients produced in vglm?
Attached is the code and output — see comment added to output to
show where I need p-values
+ print(paste("** Using VGAM functi
On Fri, Jun 5, 2009 at 4:09 AM, Allan Engelhardt wrote:
> I'm all done now. The "max2" version below is what I went with in the end
> for my proposed change to caret::nearZeroVar (which used the "sort" method).
> Max Kuhn will make it available on CRAN soon. It speeds up that routine by
> a fac
On Jun 5, 2009, at 3:56 PM, Marc Belisle wrote:
Hi there,
I have 3 numeric variables: day (e.g., 05), month (e.g., 06), year
(e.g.,
2009).
I would like to create a (string) variable of the following form:
month/day/year (e.g., 06/05/2009).
I would be grateful to anyone who could point me
Martin Morgan wrote:
> Hi Vitalie --
>
> Vitalie S. wrote:
>
>> Dear UseRs,
>>
>> A simple class inheritance example:
>>
>>
>>> setClass("test",representation(a="numeric"))
>>> setMethod("initialize","test",
>>>
>> function(.Object,x,...){
>> print("Initializa
Thanks Prof. Ripley,
I installed the current version curl-7.19.2,
still got the same problem.
I am installing RCurl_0.97-3 and my system
is x86_64/CentOS 5.
Looking forward to more suggestions/ideas.
--- On Fri, 6/5/09, Prof Brian Ripley wrote:
From: Prof Brian Ripley
Subject: Re: [R]
Marc Belisle wrote:
> Hi there,
>
> I have 3 numeric variables: day (e.g., 05), month (e.g., 06), year (e.g.,
> 2009).
>
> I would like to create a (string) variable of the following form:
> month/day/year (e.g., 06/05/2009).
>
> I would be grateful to anyone who could point me toward a solution.
>
My apologies if this is not the correct venue for this type of question.
I am looking for a software consultant to assist in our project to
integrate R with an existing Java data collection application. So far,
I have not been able to find a directory of software consultants who
have R experi
Rainer M Krug wrote:
Hi
I want to create a new generic method, but I end up with an error
(evaluation nested too deeply). see the transcript below.
The function beginYear.Fun() works, but not beginYear.
I have no idea why.
Any ideas welcome,
Rainer
setClass("fun", representation(x = "num
At least for chi-plot:
http://lmgtfy.com/?q=genest+%20copula+r-project+%22chi-plot%22
Genest and Favre cite R in that article. Have you contacted them?
On Jun 5, 2009, at 5:29 PM, Zhao Nan wrote:
Dear R users,
I was trying to find a library for the Chi-plot and Kendall plot,
as described in
Rainer M Krug writes:
> Hi
>
> I want to create a new generic method, but I end up with an error
> (evaluation nested too deeply). see the transcript below.
> The function beginYear.Fun() works, but not beginYear.
>
> I have no idea why.
>
> Any ideas welcome,
>
> Rainer
>
>
>> setClass("fun", re
How do you fit a Weibull distribution in R?
__
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,
Hi,
Is anyone aware of a rank-based, non-paired test such as the Krustal-Wallis,
that can accommodate weights?
Alternatively, would it make sense to simulate a dataset by duplicating
observations in proportion to their weight, and then using the Krustal-Wallis
test?
thanks!
Dylan
___
"Vitalie S." writes:
> Dear UseRs,
>
> Does anyone know when exactly the validity is checked in S4?
> Documentation is silent:(.
It gets called as part of initialize,ANY-method, i.e., at the bottom
of your 'callNextMethod'. It is NOT called when there are no arguments
in addition to .Object, e.
At 04:02 PM 6/5/2009, Christophe Genolini wrote:
Hi the list
Strictly speaking, this is not a "R" question, but I need the
information for the
creation of a package. My question is about vocabulary: What is the
difference between
"time series", "longitudinal data" and "trajectories"?
Sincere
Wacek Kusnierczyk writes:
> Martin Morgan wrote:
>> Hi Vitalie --
>>
>> Vitalie S. wrote:
>>
>>> Dear UseRs,
>>>
>>> A simple class inheritance example:
>>>
>>>
setClass("test",representation(a="numeric"))
setMethod("initialize","test",
>>> function(.Object
On Fri, 5 Jun 2009, Dylan Beaudette wrote:
Is anyone aware of a rank-based, non-paired test such as the Krustal-Wallis,
that can accommodate weights?
You don't say what sort of weights, but basically, no.
Whether you have precision weights or sampling weights, the test will no
longer be distr
Is there a circular shift function in R similar to the shift function in
Octave? From my reading of the manual there does not appear to be one, but
if I am mistaken, can someone point out the relevant section?
--
View this message in context:
http://www.nabble.com/Circular-shift-function--tp2389
http://lmgtfy.com/?q=fit+weibull+distribution+R
Hit #8:
http://cran.r-project.org/doc/contrib/Ricci-distributions-en.pdf
On Fri, Jun 5, 2009 at 7:21 PM, wrote:
>
> How do you fit a Weibull distribution in R?
>
> __
> R-help@r-project.org mailing lis
On Friday 05 June 2009, Thomas Lumley wrote:
> On Fri, 5 Jun 2009, Dylan Beaudette wrote:
> > Is anyone aware of a rank-based, non-paired test such as the
> > Krustal-Wallis, that can accommodate weights?
>
> You don't say what sort of weights, but basically, no.
>
> Whether you have precision weig
Martin Morgan wrote:
[...]
>
>> it sounds approximately ok that defining a subclass fails if its
>> initalizer does not provide a value required by the initializer of the
>> superclass. pardon me my ignorance, but here's an obvious question:
>> what if i *do* want to extend a class that has an
?str
It's a pretty fundamental function for understanding R objects.
> str(summary(vglm.D93))
Formal class 'summary.vglm' [package "VGAM"] with 43 slots
..@ coef3 : num [1:5, 1:3] 3.04 -4.54e-01 -2.93e-01
1.00e-15 9.33e-16 ...
.. ..- attr(*, "dimnames")=List of 2
.. .. ..$ : chr
Wacek Kusnierczyk wrote:
> Martin Morgan wrote:
>
> [...]
>
>>> from a contiguous sequence of integers, starting at 1 (this doesn't seem
>>> to be a particularly extravagant idea either):
>>>
>>> setClass('foo', representation(id='integer'))
>>> setMethod('initialize', 'foo', local({
>>>
An interview with Chief Marketing Officer of SAS Institute, Jim Davis.
Here is an extract-
http://www.decisionstats.com/2009/06/05/interview-jim-davis-sas-institute/
*Ajay -It is rare to find a major software company that has zero involvement
with open source movement (or as I call it with peer-
Dear ryan,
See the third example in
require(survival)
?survreg
HTH,
Jorge
On Fri, Jun 5, 2009 at 6:21 PM, wrote:
>
> How do you fit a Weibull distribution in R?
>
> __
> R-help@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-
Dear R experts..
I would highly appreciate your help in troubleshooting errors as shown
below..
When I type R at shell prompt, I get the following error:
The version of R in 64bit Linux is: R version 2.9.0 (2009-04-17)
***
On Fri, Jun 05, 2009 at 08:27:45PM -0700, Santosh wrote:
>
> When I run the following:
> nm /usr/lib64/R/library/stats/libs/stats.so | grep _gfortran_copy_string
>
You probably need to install libraries related to gfortran.
__
R-help@r-project.org mail
On Sat, 6 Jun 2009, Zeljko Vrba wrote:
On Fri, Jun 05, 2009 at 08:27:45PM -0700, Santosh wrote:
When I run the following:
nm /usr/lib64/R/library/stats/libs/stats.so | grep _gfortran_copy_string
You probably need to install libraries related to gfortran.
Yes, specifically libgfortran. Loo
74 matches
Mail list logo
