First of all, thanks so much! In fact I had looked at the site
before. But, that site leads me to the following site:
http://www.vsthost.com/vstDocs/mgarchBEKK/
but it was not clear to me where I can find the package to download &
compile. Is there some special website I need go to in order t
-Halcyon- gmail.com> writes:
>
>
> Hi everyone,
>
> say i have a population (stable) with different amounts of animals in every
> ageclass (80 of age 1, 60 of age 2, etc) in a vector.
> Can anybody tell me how i can add gender (male or female) to all ageclasses?
> I want a 1:1 ratio of males
Tom La Bone gforcecable.com> writes:
>
>
> I installed EMACS-ESS on Windows XP-Pro and it worked with R perfectly right
> out of the box. I was impressed with how easy it is to use (I normally use
> Tinn-R). I then switched over to Ubuntu 7.10 and installed EMACS-ESS.
> Everything worked the sa
R-users
E-mail: r-help@r-project.org
Murdoch>You may have a binary install of the package; you need the source.
Look
Murdoch>in http://cran.r-project.org/src/contrib/gam_0.98.tar.gz, files
Murdoch>gam/src/backfit.f and gam/src/backlo.f
Thank you very much.
I obtained the source codes of backfi
R-users
E-mail: r-help@r-project.org
My R is R2.61 + fields 4.1 + spam0.31-1.
I tried Tps() in fields like
function ()
{
library(fields)
xx <- matrix(c(1,2,3,2,4,5,10,30,20,40,50,90), ncol=2)
print(xx)
yy <- c(5,1,2,3,4,5)
out1 <- Tps(xx,yy, lambda=1000)
print(out1$fitted.values)
print(out1
I installed EMACS-ESS on Windows XP-Pro and it worked with R perfectly right
out of the box. I was impressed with how easy it is to use (I normally use
Tinn-R). I then switched over to Ubuntu 7.10 and installed EMACS-ESS.
Everything worked the same as in Windows (which is my main reason for using
On 24/12/2007 6:59 PM, Kunio takezawa wrote:
> R-users
> E-mail: r-help@r-project.org
>
>>> I found the answer myself.
>>> '.Fortran("baklo",' in lo.wam() and .Fortran("bakfit",in
>>> s.wam() may carry out backfitting. But I cannot
>>> create an R code which gives the same results as those of
>
R-users
E-mail: r-help@r-project.org
>> I found the answer myself.
>> '.Fortran("baklo",' in lo.wam() and .Fortran("bakfit",in
>> s.wam() may carry out backfitting. But I cannot
>> create an R code which gives the same results as those of
>> "bakfit". If someone knows the detail of "bakfit" alg
Hi, Young,
pls see if attached is what you are looking for.
http://www.google.com/search?q=m-garch&domains=r-project.org&sitesearch=r-project.org&btnG=Google+Search
On 12/24/07, Young Cho <[EMAIL PROTECTED]> wrote:
> Is there a package or function for multivariate garch model in R? I am
> having a
Hi, Sandeep,
what is your specification of p, d, q? it is not surprising to have
all forecasted with same value.
On 12/24/07, Sandeep Nikam <[EMAIL PROTECTED]> wrote:
> Hi,
>
> This is regarding the ARIMA model.
>
> I am having time series data of stock of 2000 values. Using the ARIMA model
> in
Hi,
This is regarding the ARIMA model.
I am having time series data of stock of 2000 values. Using the ARIMA model
in R, I want the forcasted values for next 36 time points.
However when I run this model in R, I am getting same value for all 36 time
points.
I have tried to fit the data with ARI
See summaryBy in the doBy package.
On Dec 24, 2007 4:27 PM, Agrarimmobilien
<[EMAIL PROTECTED]> wrote:
> Hello,
>
> using the syntax
>
> aggregate(daten[,c(3,4)], list(A,B), mean)
>
> I'm getting the following data.frame:
>
>
>A BCD
> 1 351 6.16000
Hello,
using the syntax
aggregate(daten[,c(3,4)], list(A,B), mean)
I'm getting the following data.frame:
A BCD
1 351 6.16000 5
2 47131.24333 20
3 54126.81773 2
4 3 212.990
Dear R Users,
In the expresso function, which combination of these methods for data
pre-processing (when using affymetrix oligo arrays) is the best:
bgcorrect.metod = rma rma2 mas
normalize.method = qspline quantiles loess
pmcorrect.method = pmonly subtractmm mas
summary.method = liwong avgdif
Taka
You might want to try R-SIG-Debian for this kind of question. Here are
instructions from Dirk in response to my similar question a few days ago.
HTH,
Paul
_
On 12 December 2007 at 18:17, Paul Gilbert wrote:
| I'm trying to build R from source on Ubuntu Gutsy Gibbon. I've do
Is there a package or function for multivariate garch model in R? I am
having a hard time in locating one. Thanks for help in advance.
-Young
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman
On Dec 24, 2007 1:08 PM, Syed Abid Hussaini <[EMAIL PROTECTED]> wrote:
> Thanks Gabor!
> I never realized there could be a "diff" function. all i was trying to do
> was to calculate the
> iterated differences in a vector. I feel so embarassed, but so happy for
> using R.
> by the way, i used
>
Hi
I bought a new laptop HP dv9500 just a week ago and installed a Ubuntu
gutsy ribbon on this laptop. I wanted to install Fedora but there are
more threads on Ubuntu, so I decided to install Ubuntu. After hours of
struggle in configuring x server/graphic card stuff, I installed R for
Gutsy ribbon
Thanks Gabor!
I never realized there could be a "diff" function. all i was trying to do
was to calculate the
iterated differences in a vector. I feel so embarassed, but so happy for using
R.
by the way, i used
"wave2 <- sort (wave2, decreasing=F)" to sort the values non-decreasingly
because
Hi all:
I am trying to draw a horizontal line along the zero
"Y axis" value but since zero isn't showing therefore
the line is not drawn. If I set my intercept to 15
then it'll work but I dont want to hardcode it because
the Y axis parameters could be different for another
variable. I would like t
Hi Henrique,
I read in the matrix as x
and replicated your code:
f <- function(x)
{
cbind.data.frame(chr=unique(x$chr),
Start=min(x$pos),
End=max(x$pos),
Rows=nrow(x),
Pattern=paste("(", x$s1, x$s2, ")")
)
}
do.call(
My object df is this:
df <- structure(list(BAC = structure(c(13L, 3L, 8L, 14L, 12L, 4L, 2L,
5L, 7L, 9L, 11L, 10L, 6L, 1L), .Label = c("CTD-2003M22", "RP11-155C15",
"RP11-198H14", "RP11-210E16", "RP11-210F8", "RP11-218N6", "RP11-263L17",
"RP11-267M21", "RP11-340F16", "RP11-474G23", "RP11-68A1", "R
Dear Henrique,
Thanks a lot for the help! I got this:
> f <- function(x)
+ {
+ cbind.data.frame(chr=unique(x$chr),
+ Start=min(x$pos),
+ End=max(x$pos),
+ Rows=nrow(x),
+ Pattern=paste("(", x$s1, x$s2, ")")
+ )
+ }
>
Try this:
f <- function(x)
{
cbind.data.frame(chr=unique(x$chr),
Start=min(x$pos),
End=max(x$pos),
Rows=nrow(x),
Pattern=paste("(", x$s1, x$s2, ")")
)
}
do.call("rbind", lapply(lapply(split(df, paste(df$s1, df$s2
Thanks Moshe! I apologize for not being so clear about the
second part. Again, below is how the data looks like. The
pattern for columns s1 and s2 will be:
(-1 -1) (-1 0) (-1 1) (0 -1) (0 0) (0 1) (1 -1) (1 0) (1 1)
104131 57 631 30566833 15 107
Th
Your model is singular. Varying m and log(l) have the same
effect: log(ir) = log(k) + m * log(l) * ox
Also with plinear you don't specify the linear coefficients but
rather an X matrix whose coefficients represent them:
If we use this model instead:
ir = k * exp(m * ox)
Then:
> mod0 <- lm(log(
I'm trying to fit a function y=k*l^(m*x) to some data points, with reasonable
starting value estimates (I think). I keep getting "singular matrix 'a' in
solve".
This is the code:
ox <- c(-600,-300,-200,1,100,200)
ir <- c(1,2.5,4,9,14,20)
model <- nls(ir ~ k*l^(m*ox),start=list(k=10,l=3,m=0.004
> "ML" == Mark Leeds <[EMAIL PROTECTED]>
> on Fri, 21 Dec 2007 19:38:19 -0600 (CST) writes:
ML> I was playing around with a simple example using solve.qp ( function is
in the quadprog package ) and the code is below. ( I'm not even sure there if
there is a reasonable solution bec
zoo objects are intended to represent time series and are
based on vectors and matrices, like ts objects, not data frames.
See ?zoo
Create your new columns as numeric variables:
> library(zoo)
> library(chron) # need chron 2.3-16 for month.day.year
> z <- zoo(matrix(1:24, 6, 4), Sys.Date() + seq(
R
I get a daily feed of data over the internet that I keep in various .csv
files.
I have built a function that reads that data into R for me:
getMarketData<-function(market)
{
library(zoo)
pathname<- "C:/DATA/"
files<-c("AN_REV.csv","AX_REV.csv","BN_REV.csv")
ma
> "PD" == Peter Dalgaard <[EMAIL PROTECTED]>
> on Sat, 22 Dec 2007 11:29:30 +0100 writes:
PD> Chris Stubben wrote:
>> Also try the odiag function in the demogR package
>>
>> odiag( 1:5, -1)
>> [,1] [,2] [,3] [,4] [,5] [,6]
>> [1,]000000
anacor-0.9.0 is on CRAN (by De Leeuw and Mair)
anacor does correspondence analysis and canonical correspondence
analysis. It can make row plots, column plots,
joint plots, Benzécri plots, regression plots,
and transformation plots. Where appropriate, plots can
be in 3d using either rgl or scatterp
Tyler Smith ha scritto:
> On 2007-12-21, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote:
>> Dear R users,
>>
>> I have just moved to R2.6.1 under Opensuse linux 10.3. I used to
>> work with R under XPpro. Is it "normal" to have a visual aspect of R
>> under linux different ?
>
> Yes, that's normal.
To answer your firs question try
M[-which( M$s1 == 0 & M$s2 == 0),]
For the second question, you must start with the more
precise definition of the grouping criterion.
--- affy snp <[EMAIL PROTECTED]> wrote:
> Hello list,
>
> I have a data frame M like:
>
> BAC chrpos
34 matches
Mail list logo