Please don't lie! You falsely claimed:
In R version 2.7.1, Brian Ripley adopted Wichmann's 1984 correction in
RNG document.
The period in the help has been unchanged since at least Jan 2000 (since
Random.Rd has not had any changes to those lines since then).
Your coincidence calculations m
R-devel is the list for enhancement requests.
On Tue, 19 Aug 2008, Richard M. Heiberger wrote:
Currently (R-2.7.1 on Windows)
bringToTop(stay=TRUE)
when run before a device is opened, gives the error message
bringToTop(stay=TRUE)
Error in bringToTop(stay = TRUE) :
can only bring windows de
On Mon, 18 Aug 2008, Zornitsa Luleva wrote:
Hello,
I am trying to call a FORTRAN subroutine within R and something really
strange happens:
Did you read the posting guide? It suggests the R-devel list for
questions about 'C, C++ etc', and Fortran is part of the 'etc'.
I have a dll-library
Currently (R-2.7.1 on Windows)
bringToTop(stay=TRUE)
when run before a device is opened, gives the error message
> bringToTop(stay=TRUE)
Error in bringToTop(stay = TRUE) :
can only bring windows devices to the front
>
In the same circumstances, a call to par() opens a device.
My request is
On 19/08/2008, at 12:55 PM, ascentnet wrote:
I have a simple X, Y data frame that I am trying to run regression
analysis
on. The linear regression looks great, but when I use lm(formula =
y ~
poly(x, degree = 5)) I get the same coeffecients. So for example
if I use
degree =3 my formula
You can only have one set of column names and one set of row names per
matrix. This is also true for data.frame, but for data.frame you can
at least use one of the columns to store a 2nd set of "row names".
Typically people store extra annotation like what you are asking for
externally and then a
Please forgive the cross posting. I sent this to Medstats and have not received
any response.
I am looking for a good written explanation of components of variance (a.k.a
variance components), in particular as they are used in genetic analyses. Does
anyone have a suggestion? I am interested in
I have a simple X, Y data frame that I am trying to run regression analysis
on. The linear regression looks great, but when I use lm(formula = y ~
poly(x, degree = 5)) I get the same coeffecients. So for example if I use
degree =3 my formula would look like y = 4.2 x^3 + 3.2x^2 + 2.1x + 1.0 and
Vose Consulting is pleased to announce the following course:
Quantitative Risk Analysis with R
Fort Collins, Colorado, USA
October 27-30, 2008
This 4-day course will cover the core principles of quantitative risk
analysis and the most important risk modeling principles, methods and
techniques. Th
A, I'm doing this because of the project I'm working on. The gene
names in clinical genomic usually have two specifications for one
entry. And I was just wondering whether there is a simple way of
solving that problem.
2008/8/19 jim holtman <[EMAIL PROTECTED]>:
> The row names are easy since t
Michael
set.seed(1)
boxplot(list(A=rnorm(50, 2), B=runif(50)*2))
text(2, 3, 'Text here')
Presumably this is not what you mean. Can you be a little more
explicit?
Peter Alspach
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of Michael Grant
> Sent
Hi Alberto,
Please disregard my previous note - I probably had a black-out!!!
--- On Tue, 19/8/08, Alberto Monteiro <[EMAIL PROTECTED]> wrote:
> From: Alberto Monteiro <[EMAIL PROTECTED]>
> Subject: [R] A doubt about "lm" and the meaning of its summary
> To: r-help@r-project.org
> Received: Tue
Hi Alberto,
In your second case the linear model
y = a*x + b + error
does not hold.
--- On Tue, 19/8/08, Alberto Monteiro <[EMAIL PROTECTED]> wrote:
> From: Alberto Monteiro <[EMAIL PROTECTED]>
> Subject: [R] A doubt about "lm" and the meaning of its summary
> To: r-help@r-project.org
> Recei
Hans,
I was not aware of the Orange library. Since import via PMML
is not possible I guess I could simply build the model in Python
using this library. I'll give this a shot and compare it with
the approach via Rpy.
Thanks for the help!
Bob
On Sun, 2008-08-17 at 23:31 -0700, Hans W. Borche
Is there anyone who can help me with the following? It's hard to describe in a
few words, but hope you'll understand.
I attempt to model refugee mortality over time using Monte Carlo. Before doing
repetitive simulations, my function should go through only one survey/row at
the time, but instea
Is there a easy way to insert text inside a boxplot as can be done with
plot using text(...)?
Thanks,
MCG
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
Hello,
I am trying to setup a menu to open files and directories. I have the
following code:
opendir<-function() {
dirname<<-tclvalue(tkchooseDirectory())
}
openfile<-function() {
filename<<-tclvalue(tkgetOpenFile())
}
require(tcltk)
t1<-tktoplevel()
topMenu<-tkmenu(t1)
tkconfigure(
Hello,
I have compared the potentials of R and Stata about GLMM, analysing the dataset
'ohio' in the package 'faraway' (the same dataset is analysed with GEE in the
book 'Extending the linear model with R' by Julian Faraway).
Basically, I've tried the 2 commands 'glmmPQL' and 'glmer' of R and the
I am trying to use the type="b" in the plot statement so that there
are indeed points where there is a real data point and then the linear
interpolation is between "points" that are seperated by NA. The
symbol shows that only where the symbol is is where there is real
data, and I planned on relayi
well i see i've been really foolish, sorry & thank you!
Juan
On Mon, Aug 18, 2008 at 9:09 PM, Henrik Bengtsson <[EMAIL PROTECTED]>wrote:
> ...and see saveObject() and loadObject() in R.utils, which might do
> what you expected, e.g.
>
> > x <- stats::runif(20)
> > y <- list(a=1, b=TRUE, c="oops"
...and see saveObject() and loadObject() in R.utils, which might do
what you expected, e.g.
> x <- stats::runif(20)
> y <- list(a=1, b=TRUE, c="oops")
> saveObject(list(foo=x,bar=y), "xy.Rdata")
> xy <- loadObject("xy.Rdata")
> str(xy)
List of 2
$ foo: num [1:20] 0.1154 0.6534 0.0876 0.5977 0.684
Read the help for load(); note the ``Value'' section.
cheers,
Rolf Turner
On 19/08/2008, at 11:51 AM, Juan Manuel Barreneche wrote:
Hello list! i write because i'm having a strange problem with the
"save"
function, here an example:
x <- stats::runif(20)
y <- list(
Hello list! i write because i'm having a strange problem with the "save"
function, here an example:
> x <- stats::runif(20)
> y <- list(a = 1, b = TRUE, c = "oops")
> save(x, y, file = "xy.Rdata")
> xy <- load("xy.Rdata")
> xy
[1] "x" "y"
why does this happens, when:
> x
[1] 0.96586833 0.6034050
>> Doesn't a line plot inherently display a set of linear interpolations?
>
> Yes. And your point is?
>
> Compare:
>
> x <- 1:10
> y <- rep(1:2,5)
> y[5] <- NA
> y0 <- approx(x,y,xout=1:10)
> plot(y0,type="l")
>
> with
>
> foo(x,y)
>
> where
>
> foo <- function(x,y,...) {
> plot(x,y,type="n",...)
...or use the more generic substitute() to replace parts of an expression, e.g.
i <- 3;
xlab1 <- substitute(g[idx], list=list(idx=i));
xlab2 <- bquote(g[.(i)]);
stopifnot(identical(xlab1, xlab2));
/Henrik
On Mon, Aug 18, 2008 at 3:18 PM, Marc Schwartz
<[EMAIL PROTECTED]> wrote:
> on 08/
Hello,
The simple example below plots three species by year and site. However,
I would like the points to occur in groups horizontally, rather than
on top of each other.
I thought jitter might work, but haven't been successful in getting it
to do what I want.
Any ideas? Thanks in advanc
Hi,
see reassignInPackage() of R.utils - that's where the code snipped
comes from. I wrote that so that I could patch/optimize certain
function (before they where patch in the public release). I haven't
used it to patch S4 methods, so I don't know it will work for them.
Here is an example of how
On Sun, Aug 17, 2008 at 2:36 AM, Dieter Menne
<[EMAIL PROTECTED]> wrote:
> Michael Braun MIT.EDU> writes:
>
>>
>> Dieter:
>>
>> Thank you for your response. As you requested, I created a self-
>> running example, pasted below. It may be a little wordier than I
>> would like, but it runs.
>
> ..
Thank you for the list of references. Do you know of any "free"
references available online? I'll have to find my library card :-)
My motivation is to try to correct for a "time on board" bias we see in
our surveys. Not surprisingly, riders who are only on board a short
time don't attempt/fini
Hi,
I am having a problem with a fixed-effects regression using the "plm"
package. A regression of this form runs just fine:
plm(y ~ x*z, data=datasubset, model="within")
but when i try to add a lag in there, and run it on the same dataset,
like this:
plm(y ~ x*z + lag(x,-1)*z, data=datasu
On 19/08/2008, at 10:15 AM, hadley wickham wrote:
On Mon, Aug 18, 2008 at 3:54 PM, Rolf Turner
<[EMAIL PROTECTED]> wrote:
On 19/08/2008, at 1:04 AM, stephen sefick wrote:
The real data are counts of aquatic insects at distinct locations
on a
river continuum on the way down a river, so th
on 08/18/2008 05:00 PM Nanye Long wrote:
> Hi all,
>
> I want to do plot() in a loop to make 10 graphs, so I have some code like
>
> for (i in 1:10) {
>plot(... ... , xlab = expression(g[i]) )
> }
>
> I expect g_1, g_2, and so on appear on x labels, but it simply prints
> g_i for each graph.
On Mon, Aug 18, 2008 at 3:54 PM, Rolf Turner <[EMAIL PROTECTED]> wrote:
>
> On 19/08/2008, at 1:04 AM, stephen sefick wrote:
>
>> The real data are counts of aquatic insects at distinct locations on a
>> river continuum on the way down a river, so the when a point point on
>> the graph is missing
Nanye Long said the following on 8/18/2008 3:00 PM:
Hi all,
I want to do plot() in a loop to make 10 graphs, so I have some code like
for (i in 1:10) {
plot(... ... , xlab = expression(g[i]) )
}
I expect g_1, g_2, and so on appear on x labels, but it simply prints
g_i for each graph. Does
Hi all,
I want to do plot() in a loop to make 10 graphs, so I have some code like
for (i in 1:10) {
plot(... ... , xlab = expression(g[i]) )
}
I expect g_1, g_2, and so on appear on x labels, but it simply prints
g_i for each graph. Does anybody know how to get around this problem?
Thanks.
N
Hi,
since you are new to R, you might not be aware of some nice feature at
the homepage of R: the task view
http://cran.r-project.org/web/views/
For various topics, somebody kindly provided an overview of packages and
functions for the specific topic.
Also for cluster analysis:
http://cran.r-
Using Mathematica, I've performed a hierarchical clustering
of 3,107 U. S. counties based on 1995-2000 intercounty
migration, using an algorithm of my own devising
(see http://arxiv.org/abs/0807.1550). I can also generate
the associated 3,107 x 3,107 matrix of ultrametric distances.
Unfortunately
On 19/08/2008, at 1:04 AM, stephen sefick wrote:
The real data are counts of aquatic insects at distinct locations on a
river continuum on the way down a river, so the when a point point on
the graph is missing (lost sampling equipment for that month) the
points downstream (toward the right) a
See dmt function in csampling package.
On Mon, Aug 18, 2008 at 3:15 PM, Andreia da Silva Meyer
<[EMAIL PROTECTED]> wrote:
> Hi,
>
> I am looking for an efficient way to compute the values of the density
> function of a multivariate T distribution - something like "dmvnorm", but
> for T distr. Does
Hi,
I am looking for an efficient way to compute the values of the density
function of a multivariate T distribution - something like "dmvnorm", but
for T distr. Does this exist somewhere?
Many thanks,
Andréia.
__
R-help@r-project.org mailing list
htt
Hi,
I have forgotten something: in case you don't want to produce postscript
or pdf files, have also a look at the nice article by Paul Murrell in R
News, Vol 4. No. 2:
http://cran.r-project.org/doc/Rnews/Rnews_2004-2.pdf
Best,
Roland
Roland Rau wrote:
Hi,
milton ruser wrote:
we choose a
Hi,
milton ruser wrote:
we choose an family of fonts. I tryed par(family="times")
without success.
what about this:
pdf("plot1Times.pdf", family="Times")
plot(1,1)
dev.off()
pdf("plot1NotTimes.pdf")
plot(1,1)
dev.off()
Maybe you check also the help for
?postscriptFonts
I hope this helps yo
On Mon, Aug 18, 2008 at 9:27 AM, Ben Bolker <[EMAIL PROTECTED]> wrote:
> Fränzi Korner oikostat.ch> writes:
>> Have you got an answer? I am now working with R 2.7.1.. This R-version does
>> not give the scale parameter in the summary of an lmer-object.
>> Therefore, I would like to calculate the
On Mon, Aug 18, 2008 at 11:20 AM, Aaron Mackey <[EMAIL PROTECTED]> wrote:
> I have a fairly large model:
>> length(Y)
> [1] 3051
>> dim(covariates)
> [1] 3051 211
> All of these 211 covariates need to be nested hierarchically within a
> grouping "class", of which there are 8. I have an accessor
Hi,
Maybe you can specify the length of vector:
seq(from=0, by=4, l = 4)
seq(from=0, by=4, l = 5)
On Mon, Aug 18, 2008 at 3:38 PM, Timothy W. Hilton <[EMAIL PROTECTED]> wrote:
> Hello,
>
> Using seq, I would like to specify a minumum end value, rather than a maximum
> end value. For example,
Dear all,
I know that it is a know issue, but I would like to change
the type of font on my plot, and I am not sure the
rigth way. I would like to use Times New Roman font,
but according to the "par()" help, some device allow
we choose an family of fonts. I tryed par(family="times")
without succes
On Mon, 2008-08-18 at 14:31 +0100, Denise Xifara wrote:
> Thank you very much Stephen, but how will aggregate deal with months that
> fall outside annual quarters? eg, one extra month at the end of the dataset?
[Without your data I'm kind of guessing at the exact format and problem,
but the exampl
Hello,
Using seq, I would like to specify a minumum end value, rather than a maximum
end value. For example, rather than
> seq(from=0, to=10, by=4)
[1] 0 4 8
I would like to obtain
[1] 0 4 8 12
I can do that with
> by.value = 4
> seq(from=0,by=by.value,to=ceiling(10 / by.value)*by.value)
[1
I have a conceptual problem (sort of). Maybe there's a simple solution,
maybe not.
First, let me explain the test case that goes ok.
Let x be a (fixed) n-dimensional vector. I simulate a lot of linear models,
y = m x + c + error, then I do a lot of regressions. As expected, the
estimated value
Your reading, in increasing order of difficulty/mathematical details,
might be Lohr's "Sampling"
(http://www.citeulike.org/user/ctacmo/article/1068825), Korn &
Graubard's "Health Surveys"
(http://www.citeulike.org/user/ctacmo/article/553280), and Sarndal et.
al. Survey Math Bible
(http://www.citeul
Sorry a correction to my last posting. I had accumulate switched
between prod and cumprod and I had also forgotten to included time for
conversion from list back to matrix for cumprod. Now as Chuck stated
the results for Reduce are about the same or worse than a loop.
regards--jeff
Thanks for the tips on inline, jit and Reduce. The latter was exactly
what I wanted although the loop
is still the fastest for the simple product (accumulate=TRUE for
reduce). With regards to Moshe's comment,
I was just surprised by the timing difference. I tend to use apply
without giving it
I have a fairly large model:
> length(Y)
[1] 3051
> dim(covariates)
[1] 3051 211
All of these 211 covariates need to be nested hierarchically within a
grouping "class", of which there are 8. I have an accessory vector, "
cov2class" that specifies the mapping between covariates and the 8 classes
I'm trying to learn how to calibrate/postStratify/rake survey data in
preparation for a large survey effort we're about to embark upon. As a
working example, I have results from a small survey of ~650 respondents,
~90 response fields each. I'm trying to learn how to (properly?) apply
the aforemen
Robin,
You can use Lag (upper case "L") from package Hmisc.
Nestor Arguea
On Aug 18, 2008, at 8:40 AM, Williams, Robin wrote:
Dear all,
I am having difficulties using the seemingly-simple function lag.
I have a dataframe with several weather variables (maxitemp,
windspeed, rainfall etc), and
Hi,
this is clearly a Bioconductor-specific issue. Please post your
question to the Bioconductor mailing list instead.
/Henrik
On Mon, Aug 18, 2008 at 12:20 AM, Tae-Hoon Chung <[EMAIL PROTECTED]> wrote:
> Hi, All;
>
> I simply followed the exonmap vignette but came up with an unexpected
> erro
Henrik Bengtsson stat.berkeley.edu> writes:
>
> The few times I want to patch a function like this, I use:
>
> unlockBinding(name, env);
> assignInNamespace(name, value, ns=pkgName, envir=env);
> assign(name, value, envir=env);
> lockBinding(name, env);
>
> /Henrik
>
OK, I admit th
Hi Lorenzo,
I may (?) have left something out. It isn't clear what "~" is supposed to
mean; perhaps it is just a spacer, or perhaps you meant the following:
plot(1,1)
legend("topright", expression(paste(R[g] %~~% k^{1/d[f]^{small}},~5<=k,
{}<=15)))
HTH, Mark.
Mark Difford wrote:
>
> Hi Loren
Whoops, the final var estimator var(f(Y)) should have N^4 in the
denominator not N^2
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of Doran, Harold
> Sent: Monday, August 18, 2008 10:47 AM
> To: Stas Kolenikov
> Cc: r-help@r-project.org
> Subject: Re
Hallo,
i want to fit a GARCH model with a extern regressor (without arma
components), so i found the following function in package fGarch. I tryed
out a lot of things but usually I get this Error.
> garchFit(formula=y~x, formula.var=~garch(1,1),data=w)
Error in .garchFit(formula.mean, formula.var
Applejus <[EMAIL PROTECTED]> wrote:
> I am trying to convert the following command from SPLUS to R:
>
> graphsheet(pages = TRUE)
>
> Does anyone have an idea what is the equivalent in R?
I would use
windows(record = TRUE)
which won't give a tabbed graphics device, but will allow you to
page
It also turns out that in educational testing, it is rare to consider
the sampling design and to estimate design-consistent standard errors. I
appreciate your thoughts on this, Stas. As a result, I was able to bring
to my mind more transparency into what R's survey package as well as SAS
proc surve
Hi Lorenzo,
>> ...but I would like to write that 5<=k<=15.
This is one way to do what you want
plot(1,1)
legend("topright", expression(paste(R[g]~k^{1/d[f]^{small}}~5<=k, {}<=15)))
HTH, Mark.
Lorenzo Isella wrote:
>
> Dear All,
> I am sure that what I am asking can be solved by less than a
You want to apply lag to time series objects and its likely that
you don't have that.
The dyn (and also dynlm) packages have facilities for peforming lagged
regressions on various time series objects.
If you don't want to use time series objects you can do the lag yourself so
for vector x this wi
On 8/16/08, Doran, Harold <[EMAIL PROTECTED]> wrote:
> In terms of the "design" (which is a term used loosely) the schools were
> not randomly selected. They volunteered to participate in a pilot study.
Oh, that's a next level of disaster, then! You may have to work with
treatment effect models, o
Zornitsa Luleva wrote:
> Hello,
>
> I am trying to call a FORTRAN subroutine within R and something really
> strange happens:
>
> I have a dll-library, that I load with dyn.load('mpbvv.dll'). I have checked
> the [Ordinal/Name Pointer] Table for the function within the library that I
> want to call
Fränzi Korner oikostat.ch> writes:
> Have you got an answer? I am now working with R 2.7.1.. This R-version does
> not give the scale parameter in the summary of an lmer-object.
>
> Therefore, I would like to calculate the scale parameter by hand, but I did
> not find out how.
>
One quick an
Hi Fränzi,
mod.sum <- summary(model)
mod.sum @ sigma
Hank
On Aug 18, 2008, at 9:34 AM, Fränzi Korner wrote:
Hi Oliver
you once asked this question at the R help list.
Hi all,
I try to fit a glmm model with binomial distribution and I would to
verify that the scale parameter is close to
Dear All,
I am sure that what I am asking can be solved by less than a
one-liner, but I am having a hard time to get this right.
Inside a legend environment, I do not have any problem if I issue the command:
expression(""*R[g]*"~"* k^{1/d[f]^{small}}*", "*k<= {} *"15")
but I would like to write th
Hello,
I am trying to call a FORTRAN subroutine within R and something really
strange happens:
I have a dll-library, that I load with dyn.load('mpbvv.dll'). I have checked
the [Ordinal/Name Pointer] Table for the function within the library that I
want to call - it is there (objdump - p mpbvv.dll
Dear all,
I am having difficulties using the seemingly-simple function lag.
I have a dataframe with several weather variables (maxitemp,
windspeed, rainfall etc), and the response variable (admissions). The
dataset is fairly large (1530 observations). I simply want to model the
response again
Hi Oliver
you once asked this question at the R help list.
Hi all,
I try to fit a glmm model with binomial distribution and I would to
verify that the scale parameter is close to 1...
the lmer function gives the following result :
Estimated scale (compare to 1 ) 0.766783
Bu
Thank you very much Stephen, but how will aggregate deal with months that
fall outside annual quarters? eg, one extra month at the end of the dataset?
2008/8/18 stephen sefick <[EMAIL PROTECTED]>
> ?aggregate
> may do what you want
>
> On Mon, Aug 18, 2008 at 8:19 AM, Denise Xifara
> <[EMAIL PROT
Ajay ohri wrote:
I have a data file in a KML format which is Google Earth's format for
geographic data. I need to import it into a csv file . How can I do that.
KML format is just like XML format .example below
AFAIK you will need to write your own function (using the XML package)
to extract
-- begin included message ---
Dear useRs,
I'm trying to estimate the proportion of individuals with a without a certain
recurring ailment at several times points. The data are of the survival type,
with "start"-"stop" dates and whether the individual had the ailment in that
interval.
Some cases
Is there anyone who can help me with the following? It's hard to describe in a
few words, but hope you'll understand.
I attempt to model refugee mortality over time using Monte Carlo. Before doing
repetitive simulations, my function should go through only one survey/row at
the time, but inste
?aggregate
may do what you want
On Mon, Aug 18, 2008 at 8:19 AM, Denise Xifara
<[EMAIL PROTECTED]> wrote:
> Dear R users,
>
> I have a dataframe where column is has countries, column 2 is dates
> (monthly) for each countrly, the next 10 columns are my factors where I have
> measurements for each c
Hi,
and thanks for the advices. I actually tried using 'lme', but haven't
figured out the appropriate syntax to get the random error correct with
repeated measures.
let's say i have 3 independent variables, with 'aov', i would write
something like: aov(dep_var~(indep_var1*indep_var2*indep_var3) +
The real data are counts of aquatic insects at distinct locations on a
river continuum on the way down a river, so the when a point point on
the graph is missing (lost sampling equipment for that month) the
points downstream (toward the right) are valid to show a trend.
Aren't they?
Stephen Sefic
Hi,
Bob Muenchen's book is now out on release. Its terrific for a step by step
tutorial for doing stuff in R, even if all you know is SAS and SPSS. Check
this out here http://rforsasandspssusers.com/
and at
http://www.amazon.com/SAS-SPSS-Users-Statistics-Computing/dp/0387094172/ref=pd_bbs_sr_1?i
Hi,
I have a data file in a KML format which is Google Earth's format for
geographic data. I need to import it into a csv file . How can I do that.
KML format is just like XML format .example below
Which R module with deal with an OR problem (like transportation problem
using geo coded data as in
On Mon, 18 Aug 2008, [EMAIL PROTECTED] wrote:
Ok - there is an easy work-around by just to passing the skeleton to
the objective function. But would it be possible to clarify this point
in the relist() help file please (specificaly the example I would
think), implies that attributes are copied
Trends in ge are handled by the "trend"argument.
In particular, for the so called "external trend" we use varioables which
can be at the geodata object or another object.
For instance a model fitting call could be somethoing like:
ML <- likfit..., trend= ~covar1+ covar2, ...)
where covar1 and cov
Dear R users,
I have a dataframe where column is has countries, column 2 is dates
(monthly) for each countrly, the next 10 columns are my factors where I have
measurements for each country and for each date. I have attached a sample
of the data in csv format with the data for 3 countries.
I wou
Many thanks Achim, much appreciated.
A simple, naive implementation of granger testing (see the R code below)
which does not explicitly accont for collinearities appears to give the
same results as the implementation in MSBVAR, so presumably, they do
something similar.
> granger <-function(d,
Hi All,
Another nls related problem (for background, I'm migrating a complicated
modelling package from S-plus to R).
Below I've reduced this to the minimum necessary to demonstrate my problem
(I think); the real situation is more complicated.
Two similar selfStart functions, ssA and ssB.
The
Try this:
sapply(c("csv$", "pdf$"), function(x)length(dir(path=Dir, pattern = x)))
On Mon, Aug 18, 2008 at 4:59 AM, Kurapati, Ravichandra (Ravichandra)
<[EMAIL PROTECTED]> wrote:
>> Let's say,
>
>>
>
>> Dir<-"/var/www/html/celnet/users/cxadmin/FlowDuration_DataVolume/fld_0
>
>> 1_
>
>> 08_18_08"
Hi Rhelpers,
Thanks a lot, Duncan, for your help. Yes, your solution was very much what I
was looking for. I have learnt a lot of tricks from the code snippet that you
provided. I would appreciate it if you could also show me how to create a
legend for the plot. The text for the legend can be :
Hello R-list-members,
I'm trying to model ARMA(0,2) & GARCH(1,1) process using the code below, but
according to my textbook, the estimated parameters are wrong. The MA-parameters
should be negative. (I've got the same problem using garchFit()). Can anyone
tell me what I'm doing wrong? And how ca
Hi all,
I need to test symmetry of distributions. How can I do it using R?
thank you!
giov
--
View this message in context:
http://www.nabble.com/testing-symmetry-of-distibution-tp19028200p19028200.html
Sent from the R help mailing list archive at Nabble.com.
Hi,
Im using geoR and I'm trying to do some predictions, based on an external
trend.
I'm having some problems specifying my model.control, specifically how do I
define my model, and also the source of the covariate data at the prediction
locations?
I am assuming that the covariate data at the p
Hi ...
Sorry, an "e" was erroneously "elided" from Ripley...
Mark Difford wrote:
>
> Hi Nikolaos,
>
>>> My question again is: Why can't I reproduce the results? When I try a
>>> simple anova without any random factors:
>
> Lack of a "right" result probably has to do with the type of analys
Hi Nikolaos,
>> My question again is: Why can't I reproduce the results? When I try a
>> simple anova without any random factors:
Lack of a "right" result probably has to do with the type of analysis of
variance that is being done. The default in R is to use so-called Type I
tests, for good rea
> Let's say,
>
> Dir<-"/var/www/html/celnet/users/cxadmin/FlowDuration_DataVolume/fld_0
> 1_
> 08_18_08"
>
> There will be pdf and/or csv files.
> I want to know how many *.pdf and/or csv files are there in that
directory.
>
> How can I get using R commands.
>
Thanks
K.Ravicha
Hello,
I would like my Sweave document to display messages generated with the
message() command. Sweave ignores this output and I did not find any Sweave
option, nor any property that I could change through the options() command
for changing this behavior.
Here is a chunk example.
<>=
test <- fun
On Mon, 18 Aug 2008, Kurapati, Ravichandra (Ravichandra) wrote:
> Let's say,
>
> Dir<-"/var/www/html/celnet/users/cxadmin/FlowDuration_DataVolume/fld_01_
> 08_18_08"
>
> There will be pdf and/or csv files.
> I want to know how many *.pdf files are there in that directory.
>
> How can I get using R
length(Sys.glob(file.path(Dir, "*.pdf)))
or use the 'pattern' argument of list.files()
or use system(paste("ls", file.path(Dir, "*.pdf), "| wc -l"), intern=TRUE)
On Mon, 18 Aug 2008, Kurapati, Ravichandra (Ravichandra) wrote:
Hi
Lets say,
Dir<-"/var/www/html/celnet/users/cxadmin/FlowDur
2008/8/18 Kurapati, Ravichandra (Ravichandra)
<[EMAIL PROTECTED]>:
> Dir<-"/var/www/html/celnet/users/cxadmin/FlowDuration_DataVolume/fld_01_
> 08_18_08"
>
> There will be pdf and/or csv files.
>
>I want to know how many *.pdf files are there in that
> directory.
> How can
Hi, All;
I simply followed the exonmap vignette but came up with an unexpected
error.
Here are what I did:
> library(exonmap)
Loading required package: affy
Loading required package: Biobase
Loading required package: tools
Welcome to Bioconductor
Vignettes contain introductory material. To v
Hei there,
I am trying to build a new R package and after running
R CMD check and
R CMD build
I have now my pkg.tar.gz file.
The problem is that if i try to install it using
install.packages("pkg.tar.gz",repos=NULL)
i get the error
Warning in install.packages("pkg/INLA_1.0.tar.gz", repos = NULL) :
1 - 100 of 101 matches
Mail list logo
