On Sat, 12 Feb 2011, Joshua Wiley wrote:
AFAIK CRAN does not explicitly require use of a GPL license (though
obviously the license needs to allow at least free access to the
packages and source).
GPL and the other licences are about distribution. CRAN requires a
licence to allow it to distri
you can just probably add
+ labs(x="Time Elapsed,y="Predicted Probability")
Felipe D. Carrillo
Supervisory Fishery Biologist
Department of the Interior
US Fish & Wildlife Service
California, USA
http://www.fws.gov/redbluff/rbdd_jsmp.aspx
- Original Message
> From: Shige Song
> To:
Dear Colleagues,
In the following simple ggplot2 code:
m <- ggplot(d.fig, aes(time, prob))
m + stat_summary(fun.data = "median_hilow", conf.int = .95, geom =
"smooth") + facet_wrap(~ Cohort, nrow=1) + coord_cartesian(ylim = c(0,
.03))
--
I did a fourier transform on a function in time domain to get the following
functions in frequency domain (in latex):
$Y_1[\omega] = \frac{1}{1-\phi_1 e^{-jw}}$
$Y_2[\omega] = \frac{1}{1-(\phi_1 + \phi_2)e^{-jw} +\phi_1\phi_2e^{-2jw}}$
How do I find the spectrum of this function for given $\phi
Does it make sense for you to combine the 2 data sets and do a 2-way anova with
treatment vs. control as one factor and experiment number as the other factor?
Then you could test the interaction and treatment number factor to see if they
make a difference.
--
Gregory (Greg) L. Snow Ph.D.
Stat
On Sat, Feb 12, 2011 at 4:08 PM, Paka yag wrote:
>
> Hello
>
> I have model that is: lm(y~ lag(x, -1) + lag(z, -1)
> so basically a time series regression with exogen variables
> And I want to make rolling out of sample forecasts, meaning that:
>
> I first use a subsample (e.g. 1990 -1995) for e
Hello
I have model that is: lm(y~ lag(x, -1) + lag(z, -1)
so basically a time series regression with exogen variables
And I want to make rolling out of sample forecasts, meaning that:
I first use a subsample (e.g. 1990 -1995) for estimating, then I perform a one
step ahead forecast, then I a
Thanks for the hint. Helped a lot
Gabor Grothendieck wrote:
On Sat, Feb 12, 2011 at 7:48 AM, Jose-Marcio Martins da Cruz
Functions which work with ts typically assume that a full cycle is
represented by 1 unit so if a full cycle is a week then a week must be
one unit and a day must be 1/7th
That worked great! Thanks!
--
View this message in context:
http://r.789695.n4.nabble.com/Data-manipulation-tp3302717p3303001.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailm
AFAIK CRAN does not explicitly require use of a GPL license (though
obviously the license needs to allow at least free access to the
packages and source).
On Sat, Feb 12, 2011 at 11:57 AM, Tal Galili wrote:
> Hi all,
> I'm not sure who to ask this, so I'm posting this here.
>
> I just ran:
> req
Hi,
I am new to R and I want to display longitude and latitude for a projected
map.
map.axes won't do it since it only works for unprojected map.
Can anyone help me with this?
Thanks,
Sybil
--
View this message in context:
http://r.789695.n4.nabble.com/use-map-axes-for-a-projected-map-tp3302918
thanks that worked out
tina :-)
On Fri, 11 Feb 2011 12:21:19 -0500, Ista Zahn
wrote:
Hi,
one approach is to modify
getAnywhere(print.Anova.mlm)
to return the information you want.
Best,
Ista
On Fri, Feb 11, 2011 at 7:16 AM, Bettina Kulle Andreassen
wrote:
hi,
i am not able to extract
Ummm...
The simple answer to your question:
i want the
estimates of
coefficients of each level along with the intercept term(general
effect).can anyone please help me?
is no, because it's matrhematically impossible. You need to consult a
local statistician or read up on linear mo
--- On Sun, 13/2/11, Shing Hing Man wrote:
> From: Shing Hing Man
> Subject: Re: [R] Linked List in R
> To: "Martin Morgan"
> Date: Sunday, 13 February, 2011, 4:49
> Hi Martin,
> Thanks for the pointer ! It seems quite involved to
> create a
> a linked list in the oo way. I need to read up
On Sat, Feb 12, 2011 at 11:57 AM, Tal Galili wrote:
> Hi all,
> I'm not sure who to ask this, so I'm posting this here.
>
> I just ran:
> require(KernSmooth)
> And got (I bolded the text):
> Loading required package: KernSmooth
> KernSmooth 2.23 loaded
> *Copyright M. P. Wand 1997-2009*
> Warning
On 02/12/2011 11:08 AM, Shing Hing Man wrote:
> Hi,
> I am trying to create a linked list in R by defining a class Node which has
> a instance variable Node.
>
> setClass("Node", representation(rate="numeric", nextNode="ANY"))
>
> The above works. But the following gives me a warning message.
Hi all,
I'm not sure who to ask this, so I'm posting this here.
I just ran:
require(KernSmooth)
And got (I bolded the text):
Loading required package: KernSmooth
KernSmooth 2.23 loaded
*Copyright M. P. Wand 1997-2009*
Warning message:
package 'KernSmooth' was built under R version 2.12.1
What do
Why would you do this, when lists are fundamental types in R?
---
Jeff Newmiller The . . Go Live...
DCN: Basics: ##.#. ##.#. Live Go...
Live: OO#.. Dead: OO#.. Playing
Research Engineer (Solar/Batteries O.O#. #.O#. with
Hi,
I am trying to create a linked list in R by defining a class Node which has a
instance variable Node.
setClass("Node", representation(rate="numeric", nextNode="ANY"))
The above works. But the following gives me a warning message.
setClass("Node", representation(rate="numeric", nextNode="
Will this do it for you:
> x <- read.table(textConnection(" A B C
+ 1 1 a 1999
+ 2 1 b 1999
+ 3 1 c 1999
+ 4 2 c 2001
+ 5 2 d 2001
+ 6 3 a 2004
+ 7 3 b 2004"), header = TRUE)
> closeAllConnections()
> # add a tenure column
> x$tenure <- ave(x$C, x$B, FUN = function(yr) yr - min(yr)
mathijsdevaan [Sat, Feb 12, 2011 at 03:00:18PM CET]:
>
> Hi,
>
> I have a dataset with info on individuals (B) that have been involved in
> projects (A) during multiple years (C). The dataset contains three columns:
> A, B, C. Example:
>A B C
> 1 1 a 1999
> 2 1 b 1999
> 3 1 c 1999
>
On Feb 12, 2011, at 11:16 AM, Ista Zahn wrote:
Well, I may have been too hasty there. Actually the question was
slightly improved. I think you are looking for ?contr.sum
Another issue will arise besides what the switch of the contrasts
-- The intercept is not a "general effect" but rather a
See ?contrasts and ?contr.sum
By default, R uses contr.treatment, you want contr.sum
Christophe
On Sat, Feb 12, 2011 at 6:57 AM, ATANU wrote:
>
> i am trying to fit a linear model with both continuous covariates and
> factors. When fitted with the intercept
> term the first level of the factor
Well, I may have been too hasty there. Actually the question was
slightly improved. I think you are looking for ?contr.sum
Best,
Ista
On Sat, Feb 12, 2011 at 3:55 PM, Ista Zahn wrote:
> I already responded to your identical post yesterday, asking you to
> actually ask a question. "can you please
I already responded to your identical post yesterday, asking you to
actually ask a question. "can you please help me?" is not enough.
Best,
Ista
On Sat, Feb 12, 2011 at 5:57 AM, ATANU wrote:
>
> i am trying to fit a linear model with both continuous covariates and
> factors. When fitted with th
Another option is to use StatTransfer which is a commercial data
exchange product from Circle Systems:
http://www.stattransfer.com
That's the basis of the import/export in S-PLUS:
http://www.stattransfer.com/stattransfer/developers.html
I wouldn't be surprised if it's al
Hi!
is there a way in R to check whether the outcome of two different
experiments is statistically distinguishable or indistinguishable? More
preciously, I used the wilcoxon test to determine the differences between
controls and treated subjects for two different experiments. Now I would
like to
On Sat, Feb 12, 2011 at 7:48 AM, Jose-Marcio Martins da Cruz
wrote:
>
> This question is surely trivial, sorry. I'm afraid I'm misunterpreting the
> information I got with the documentation, and I'm a little bit confused. I'm
> just an engineer with some little skills in statistics.
>
> Well, I ha
Hi,
I have a dataset with info on individuals (B) that have been involved in
projects (A) during multiple years (C). The dataset contains three columns:
A, B, C. Example:
A B C
1 1 a 1999
2 1 b 1999
3 1 c 1999
4 2 c 2001
5 2 d 2001
6 3 a 2004
7 3 b 2004
I am interested in the
On 11.02.2011 16:14, xin shi wrote:
Dear:
I am recnetly trying to install some libraries. However, I found this issue for
both my laptop and desktop even I uninstall and install it again.
I even can not update the R now.
I wonder if you have the similar issue.
Do you have a) internet co
Hi,
instead of sub, use substr in R, also look for 'which' and 'factor' in the
manual.
E.g.:
?substr
?which
> ?factor
I hope this could help you to rewrite the SPSS syntax in R.
Good luck!
Best,
Gergely
On Sat, Feb 12, 2011 at 2:08 PM, beky wrote:
>
> There is a code from SPSS Syntax
>
>
You might want to take a look at Bob Muenchen's book "R for SAS and SPSS
Users"
There is an 80 page preview at .. http://rforsasandspssusers.com/
Additionally, there is lots of documentation for getting started with R on
the CRAN website http://cran.r-project.org/
HTH
Pete
--
Vi
This question is surely trivial, sorry. I'm afraid I'm misunterpreting
the information I got with the documentation, and I'm a little bit
confused. I'm just an engineer with some little skills in statistics.
Well, I have a time series - 600 days long - with some weekly
periodicity inside. So
There is a code from SPSS Syntax
do if sub(ATVK,2,2)="01".
comp strata=1.
else if sub(ATVK,2,2)>="05" and sub(ATVK,2,2)<="27".
comp strata=3.
else if sub(ATVK,4,2)>"20" or sub(ATVK,6,2)>"20".
comp strata=4.
else if sub(ATVK,4,2)>"00".
comp strata=2.
end if.
value labels strata 1 "R 2 "Li" 3 "M"
I'm looking for Markov switching models, written in R. Does anyone know about
programs?
Thanks,
Nanda
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mail
Hello Antony,
Thank you very much for your detailed answer!
All of your points are valid and interesting to reflect upon.
Regarding your note number (6), it's a very good point - I didn't think of
it.
It might be argued that because interactive data visualization can be faster
then analytical prog
Dear Sabine
In negbin(aod), the deviance is calculated by:
# full model
logL.max <- sum(dpois(x = y, lambda = y, log = TRUE))
# fitted model
logL <- -res$value
dev <- -2 * (logL - logL.max)
(the log-Lik contain all the constants)
As Ben Bolker said, whatever the formula used for deviance, diff
Dear Gabriel,
it is possible. You can define a new class or just an object of an
already existing class.
Did you look at:
library(distr)
vignette("newDistributions")
as well as
?AbscontDistribution
?DiscreteDistribution
Please let me know if you have further questions!
Best
Matthias
On 11.
i am trying to fit a linear model with both continuous covariates and
factors. When fitted with the intercept
term the first level of the factor is treated by R as intercept and the
estimate of the effects of remaining levels(say i th level) are given as
true estimate of i th level - estimate o
I want to create natural cubic spline on my predictors and fit them to a linear
model. The DF and Degree of all the splines are the same. I noticed that I
can setup the lm() formula as y~ns(x1, df=5) + ns(x2, df=5).
My problem is that I have a huge number of predictors (hundreds). ns however
Hello
T
I want to order some calculation "result", there will be lots of
"result" that need to calculate and order
PS: the "result" is just a intermediate varible and ordering them is
the very aim
# problem:
# For fixed NT and CT, and some pair (c,n). order the pair by co
Yes, this point i can understand your suggestion and I should read "HOW TO ASK
GOOD QUESTIONS".
I'm a just new mailing list, also a R user for researching on graduate school.
Any person can make mistake, also it can be effect to other people, however a
good person is teaching good way.
Anyway, a
On Fri, 2011-02-11 at 20:51 -0500, Axel Urbiz wrote:
> Dear users,
>
> I'll appreciate your help with this (hopefully) simple problem.
>
> I have a model object which was fitted to inputs X1, X2, X3. Now, I'd like
> to use this object to make predictions on a new data set where only X1 and
> X2 a
43 matches
Mail list logo
