[R] Binning a 2 column matrix by avarages of rows.

2011-12-17 Thread ali_protocol
Newbie here. Many apologies in advance for using the incorrect lingo. I'm new to statistics and VERY new to R. I have a nx2 matrix , I want to sort the values based on the average of 2 columns and put k lowest (or highest) values in bin1, second k high/low values in bin2, and so on (bins would

[R] Problem with reproducing log likelihood estimated with ghyp package

2011-12-17 Thread Steinar
I was playing around with the ghyp package and simulated series of t-distributed variables when suddenly i was not able to reproduce the log likelihood values reported by the package. When trying to reproduce the likelihood values, I summed the log(dt(x,v)) values and it worked with some simulated

[R] help me....

2011-12-17 Thread jurassic85
Life times of one brand of light bulb were observed. Life hours of nine light bulbs were 898, 720, 354, 405, 620, 54, 229, 306, 890. There were three more light bulbs that were still burning after 1000 hours. It is common to assume the life time of a light bulb follows an exponential distribution.

[R] Rejection metho..d...

2011-12-17 Thread Bryan
Generate random variables which flow f(x)=20*x*(1-x)^3 , 0x1 using rejection method with .g(x)=1, 0x1 how??? -- View this message in context: http://r.789695.n4.nabble.com/Rejection-metho-d-tp4208138p4208138.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] help me....

2011-12-17 Thread David Winsemius
On Dec 17, 2011, at 6:36 AM, jurassic85 wrote: Life times of one brand of light bulb were observed. Life hours of nine light bulbs were 898, 720, 354, 405, 620, 54, 229, 306, 890. There were three more light bulbs that were still burning after 1000 hours. It is common to assume the life

Re: [R] help me....

2011-12-17 Thread Stephen Sefick
This definatley looks like homework. On Sat 17 Dec 2011 05:36:03 AM CST, jurassic85 wrote: Life times of one brand of light bulb were observed. Life hours of nine light bulbs were 898, 720, 354, 405, 620, 54, 229, 306, 890. There were three more light bulbs that were still burning after 1000

Re: [R] Rejection metho..d...

2011-12-17 Thread David Winsemius
On Dec 17, 2011, at 7:58 AM, Bryan wrote: Generate random variables which flow f(x)=20*x*(1-x)^3 , 0x1 using rejection method with .g(x)=1, 0x1 how??? Just as with the posting adjacent to yours in my mail client, this looks to be homework. #---

Re: [R] Help for R parallelization

2011-12-17 Thread Dirk Eddelbuettel
Alvaro, On 17 December 2011 at 12:23, Alvaro Jose Abackerli wrote: | Dear Sir/Madam, | | I am a Brazilian researcher working in São Paulo – Brazil. I’ve being using R-software for many years with my students in some engineering applications, but up to now we've being using only in stand-alone

[R] time-varying parameters kalman filter estimation problem using FKF package

2011-12-17 Thread napps22
Dear R users, I am trying to carry out MLE of the time-varying CAPM using the FKF package. My approach so far has been to try and adapt the example given in the help file found using ?fkf which demonstrates the MLE of an ARMA(2,1) model. When I attempt to run my R code (given below) I get the

[R] auto.arima from the Forecast package

2011-12-17 Thread ludovic.theate
Hi, I've got a little problem using auto.arima. I run the following command auto.arima(drivers,ic=aic,d=1,D=1,max.order=10,max.p=5,max.q=5,max.P=5,max.Q=5,stepwise=FALSE,allowdrift=FALSE) and I get the following output : Series: drivers ARIMA(0,1,1)(5,1,1)[12]

Re: [R] Retrieving samples when ussing BRugs.

2011-12-17 Thread Uwe Ligges
On 14.12.2011 10:44, Thiago Guerrera wrote: Hi Sacha, I was referring to the BRugs package, where I usually proceed as follows: result = BRugsFit (...) param_samples = samplesSample(param)(1) But (1) takes a long time for big chains. So my question was specific to BRugs, if there an

Re: [R] Hi

2011-12-17 Thread Uwe Ligges
On 14.12.2011 23:11, Trying To learn again wrote: Hi all, I have solved my question Output- matrix(0, length(x), 1) for (i in 1:length(x)){ files- paste('KT', i, '.csv', sep = '') da1- read.csv(files) Output[i,1]-da1[1,2] } 2 comments: 1: use seq:_along(x) rather than 1:length(x)

Re: [R] R error

2011-12-17 Thread Uwe Ligges
On 15.12.2011 16:42, verse123 wrote: Hi guys, I am new to R and I am bascially trying to load a library that I installed and use external data that I have. When trying to use an R package called cummeRbund (http://compbio.mit.edu/cummeRbund/), I am doing: library(cummeRbund) Loading

Re: [R] slight documentation error in stats package arima

2011-12-17 Thread Uwe Ligges
On 15.12.2011 22:10, Jan Theodore Galkowski wrote: The documentation for the arima function in the package stats has a slight error. It references: Ripley, B. D. (2002) Time series in R 1.5.0. R News, 2/1, 2–7. [1]http://www.r-project.org/doc/Rnews/Rnews_2002-1.pdf This should be: Ripley,

Re: [R] Odp: Problem with package compilation

2011-12-17 Thread Uwe Ligges
On 16.12.2011 08:04, Petr PIKAL wrote: Hi Hi, I have a R package with some functions made all of then only with R code. I use the command R CMD build to build a package that I can install on linux, windows or mac, because all the code is only R code. But I have some problems with R

Re: [R] Extraction of rules from Support Vector Machines

2011-12-17 Thread Uwe Ligges
On 16.12.2011 08:44, Paul Smith wrote: Dear All, I am using the svm() function of package e1071 for creating Support Vector Machines prediction models. As far as I understand, there is no function in this package to extract rules of prediction. Is there some other package with such a

Re: [R] Binning a 2 column matrix by avarages of rows.

2011-12-17 Thread David Winsemius
On Dec 17, 2011, at 2:47 AM, ali_protocol wrote: Newbie here. Many apologies in advance for using the incorrect lingo. I'm new to statistics and VERY new to R. I have a nx2 matrix , I want to sort the values based on the average of 2 columns and put k lowest (or highest) values in bin1,

[R] odfWeave Windows binary

2011-12-17 Thread John C Nash
As I'm working with some folk who use Windows to prepare an article / documentation, I'd like to be able to know if we can use odfWeave. It seems there is no official binary at the moment for Windows. Does anyone have a working unofficial one, if possible both win32 and win64 flavours? My

[R] R package(s) for IRT scoring?

2011-12-17 Thread Prisciandaro, James J
Hello forum, I have developed an IRT model for a set of items on a large reference sample. I now want to use the IRT parameters that were estimated in the reference sample to score new individuals, one at a time. What R package(s) would allow me to do this? Also, several covariates were

[R] cp2tform equivalent

2011-12-17 Thread George Macaulay Shaw
I'm a new R user and am looking for an R equivalent to the Matlab 'cp2tform' function (see: http://www.mathworks.com/help/toolbox/images/ref/cp2tform.html). cp2tform accepts pairs of 2D control points and produces a least-squares optimized spatial transformation from one set of points to the

[R] Rattle install issues on windows 7 and R2.14.0

2011-12-17 Thread Muhammad Abuizzah
Hi all, Installed rattle and the required GTK+, but ran into problems. did some research on the net which recomened removing GTK+ and insalling it from within R from the package RGTK2 which I did but that didn't resolve the problem. I get two different errors one is noticed as soon as I

[R] odfWeave error

2011-12-17 Thread Frank Lawrence
I am new to using odfWeave but I have encountered a problem running both the example in the help file as well as another file.  I am not sure how to correct the error. First example: library(odfWeave) filein - 'c:\\my documents\\example01_in.odt' fileout - 'c:\\my documents\\example01.odt'

Re: [R] Zellig Error Message

2011-12-17 Thread Xu Wang
Hi Abraham, Zelig has its own mailing list. You should sign up for it there. The developers are usually great about responding. Xu On Fri, Dec 16, 2011 at 10:47 AM, Abraham Mathew abmathe...@gmail.comwrote: I'm trying to calculate predicted probabilities in R with Zelig and keep getting the

[R] Looking for package 'rgp'

2011-12-17 Thread Worik R
I am interested in the package 'rgp' But install.packages(rgp) Installing package(s) into ‘/home/worik/R/x86_64-pc-linux-gnu-library/2.12’ (as ‘lib’ is unspecified) Warning message: In getDependencies(pkgs, dependencies, available, lib) : package ‘rgp’ is not available It is in cran:

[R] Should data for the linear mixed model analysis meet the three assumptions of ANOVA?

2011-12-17 Thread Junli
Hi All, I am doing linear mixed model analysis for my multi-location experiment using R package lme4. I just wonder whether I should check my data first to see whether they meet the three assumptions of ANOVA, that is, independence, normality and homogeneity. I saw a lot of examples and the

Re: [R] Looking for package 'rgp'

2011-12-17 Thread David Winsemius
On Dec 17, 2011, at 9:30 PM, Worik R wrote: I am interested in the package 'rgp' But install.packages(rgp) Installing package(s) into ‘/home/worik/R/x86_64-pc-linux-gnu- library/2.12’ (as ‘lib’ is unspecified) Warning message: In getDependencies(pkgs, dependencies, available, lib) :