Re: [R] mgarch

2007-12-24 Thread YOUNG CHO
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

Re: [R] vector distribution

2007-12-24 Thread Ben Bolker
-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

Re: [R] Problems with EMACS-ESS on Ubuntu

2007-12-24 Thread Ben Bolker
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

Re: [R] "gam()" in "gam" package

2007-12-24 Thread Kunio takezawa
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] fields 4.1

2007-12-24 Thread Kunio takezawa
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

[R] Problems with EMACS-ESS on Ubuntu

2007-12-24 Thread Tom La Bone
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

Re: [R] "gam()" in "gam" package

2007-12-24 Thread Duncan Murdoch
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 >

Re: [R] "gam()" in "gam" package

2007-12-24 Thread Kunio takezawa
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

Re: [R] mgarch

2007-12-24 Thread Wensui Liu
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

Re: [R] ARIMA problem

2007-12-24 Thread Wensui Liu
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

[R] ARIMA problem

2007-12-24 Thread Sandeep Nikam
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

Re: [R] aggregation with two statistical functions - mean and variance

2007-12-24 Thread Gabor Grothendieck
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

[R] aggregation with two statistical functions - mean and variance

2007-12-24 Thread Agrarimmobilien
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

[R] Affy Package

2007-12-24 Thread Vallejo, Roger
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

Re: [R] R Compilation error on Ubuntu

2007-12-24 Thread Paul Gilbert
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

[R] mgarch

2007-12-24 Thread Young Cho
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

Re: [R] saving while loop values to one vector

2007-12-24 Thread Gabor Grothendieck
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 >

Re: [R] R Compilation error on Ubuntu

2007-12-24 Thread Takatsugu Kobayashi
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

Re: [R] saving while loop values to one vector

2007-12-24 Thread Syed Abid Hussaini
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

[R] geom_hline

2007-12-24 Thread Felipe Carrillo
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

Re: [R] How to remove some rows from a data.frame

2007-12-24 Thread affy snp
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(

Re: [R] How to remove some rows from a data.frame

2007-12-24 Thread Henrique Dallazuanna
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

Re: [R] How to remove some rows from a data.frame

2007-12-24 Thread affy snp
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, ")") + ) + } >

Re: [R] How to remove some rows from a data.frame

2007-12-24 Thread Henrique Dallazuanna
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

Re: [R] How to remove some rows from a data.frame

2007-12-24 Thread affy snp
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

Re: [R] curve fitting problem

2007-12-24 Thread Gabor Grothendieck
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(

[R] curve fitting problem

2007-12-24 Thread pieterprovoost
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

Re: [R] using solve.qp without a quadratic term

2007-12-24 Thread Martin Maechler
> "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

Re: [R] Help with read.zoo and transform

2007-12-24 Thread Gabor Grothendieck
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] Help with read.zoo and transform

2007-12-24 Thread Alex Park
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

Re: [R] Diagonal matrix with off diagonal elements

2007-12-24 Thread Martin Maechler
> "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

[R] [R-pkgs] anacor: yet another ca package

2007-12-24 Thread Jan de Leeuw
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

Re: [R] R appearance under linux

2007-12-24 Thread Jabba
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.

Re: [R] How to remove some rows from a data.frame

2007-12-24 Thread Moshe Olshansky
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