Re: [R] suggestions argument in rbga function in genalg package

2012-04-08 Thread Frank Schwidom
On Thu, Sep 22, 2011 at 06:46:23PM +, Joseph Boyer wrote: > Would someone be so kind as to provide example code where they use the > suggestions argument in the rgba function > In genalg? I can't get it to work. > > The following code works just fine: > > GenFit <-rbga(Lower, Upper, evalFun

Re: [R] Problems with CRAN package Ryacas

2012-04-08 Thread Jeff Newmiller
check your firewall? --- Jeff NewmillerThe . . Go Live... DCN:Basics: ##.#. ##.#. Live Go... Live: OO#.. Dead: OO#.. Playing Res

Re: [R] Summary Statistics Help

2012-04-08 Thread bobo
I've got this solved via Talks Stat mod.1<-lm(Patents~FHouse, data=datpat) summary(mod.1) anova(mod.1) xtable(mod.1) -- View this message in context: http://r.789695.n4.nabble.com/Summary-Statistics-Help-tp4541923p4542103.html Sent from the R help mailing list archive at Nabble.com. _

[R] Creating a dummy variable

2012-04-08 Thread arunkumar1111
HI I want to create a dummy variable for a dataset which has various dates Eg I've one year data. The user can choose any date range startdate : 2012-01-01 enddate: 2012-02-01 startdate: 2012-06-01 enddate=2012-07-01 the columns of dataset= X ,Y ,date Thanks and re - Thanks in Adv

Re: [R] Error in integrate(int.fn, lower = 0, upper = Inf) : evaluation of function gave a result of wrong length

2012-04-08 Thread Guaramy
Sorry but i didn't understand. Thanks for your answer. -- View this message in context: http://r.789695.n4.nabble.com/Error-in-integrate-int-fn-lower-0-upper-Inf-evaluation-of-function-gave-a-result-of-wrong-length-tp4541250p4542101.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] return one vector that is the sum of a list of vectors

2012-04-08 Thread David Winsemius
On Apr 9, 2012, at 12:01 AM, Frank Tamborello wrote: Dear R Help, I am attempting to write a function that takes a list of variable groups and a vector of numbers (e.g., "jtsv" would indicate one group of variables, "jtsv1", "jtsv2", etc, each of which name a variable in a data frame),

Re: [R] Avoid loop with the integrate function

2012-04-08 Thread Berend Hasselman
On 09-04-2012, at 01:26, Navin Goyal wrote: > Hi, > I am not sure I follow the scalar part of the suggestion. Could you please > elaborate it a bit more? Each row has a different value of score for each > subject at each time point. That last remark has nothing to do with it. The point is tha

Re: [R] return one vector that is the sum of a list of vectors

2012-04-08 Thread Bert Gunter
?get ?rowSums ?rowsum a "group" is not a defined data structure in R, btw. You have a much better chance of getting precise answers if you ask precise questions. See the Posting Guide. -- Bert On Sun, Apr 8, 2012 at 9:01 PM, Frank Tamborello wrote: > Dear R Help, > > I am attempting to write a

[R] return one vector that is the sum of a list of vectors

2012-04-08 Thread Frank Tamborello
Dear R Help, I am attempting to write a function that takes a list of variable groups and a vector of numbers (e.g., "jtsv" would indicate one group of variables, "jtsv1", "jtsv2", etc, each of which name a variable in a data frame), and returns a list of variable group sums. For example, give

[R] rwmetrop

2012-04-08 Thread Yuanyuan Tang
Hi, all: Can anybody check where is wrong with my code? I tried a lot of times, but did not find an error. The parameters' estimator is not accurate. It's a simple model about a multiple regression, with five covariates. rwmetrop is supposed to give a much more accurate estimand. Thanks a lot. rm

Re: [R] Error in integrate(int.fn, lower = 0, upper = Inf) : evaluation of function gave a result of wrong length

2012-04-08 Thread David Winsemius
On Apr 8, 2012, at 1:28 PM, Guaramy wrote: Hi, i am writing a function to plot a pdf functions of a Generalized normal laplace distribution. The code is this { y = x-rho*mu cf.fn = function(s){ cplex = complex(1,0,1) temp1 = alpha*beta*e

[R] Error using PGMM function in the PLM package

2012-04-08 Thread Taylor White
Good day fellow R users: I have routinely received the following message when attempting to estimate a GMM model for a somewhat square panel (N = 20, T = 9-27, Obs = 338) using the pgmm function in the plm package: Error in function (..., deparse.level = 1) : number of rows of matrices must ma

[R] Summary Statistics Help

2012-04-08 Thread bobo
Hi, I would really appreciate all the help I can get. Unfortunately, I am really new to statistics! I hope you guys don't mind this. I am trying to find significance levels, beta, R, R squared, adjusted R squared, standard error and t test. FILE http://r.789695.n4.nabble.com/file/n4541923/datpat

[R] How To Setup hunspell in R

2012-04-08 Thread yixiang
Hello, I am doing some text mining and text analysis in R. However, I have NO idea how to install hunspell into R and use the aspell() function correctly, I keep running into this error: Error in aspell(f, "Rd") : No suitable spell check program found. Come someone give me a detailed step by

Re: [R] Changing grid defaults

2012-04-08 Thread Brett Presnell
First I'd like to thank Ilaik and Paul for their responses. I seem to be getting what I need by using LaTeX's tikz package and including tikz=true in the options for the figure chunks using vcd, e.g., <>= ... R commands ... @ I have actually used this before in other documents, but not for this

Re: [R] Drawing a line in xyplot

2012-04-08 Thread wcheckle
i appreciate all the interest in my question, and thank you Elai for both of your suggestions, which work very well. Elai, your last code was particularly simple and helpful to generate the figure i was looking for. x11(height=8,width=11) par(lend=2) xdat = data.frame(mortality =c(5, 8,7,5,8,10,11

Re: [R] Drawing a line in xyplot

2012-04-08 Thread David Winsemius
On Apr 8, 2012, at 7:24 PM, John Maindonald wrote: > This can be simplified by using the layering abilities that Felix > Andrews > made available in latticeExtra. These are too little known. These > pretty > much make it unnecessary to resort to trellis.focus(), at least in > such > cases

Re: [R] Avoid loop with the integrate function

2012-04-08 Thread Navin Goyal
Hi, I am not sure I follow the scalar part of the suggestion. Could you please elaborate it a bit more? Each row has a different value of score for each subject at each time point. Also I simplified the example but there are other terms that change over each row for each subject. Does this mean t

Re: [R] Drawing a line in xyplot

2012-04-08 Thread John Maindonald
This can be simplified by using the layering abilities that Felix Andrews made available in latticeExtra. These are too little known. These pretty much make it unnecessary to resort to trellis.focus(), at least in such cases as this. These layering abilities are too little known: library(latti

Re: [R] quantmod getOptionChain Not Work

2012-04-08 Thread Sparks, John James
Michael, I have not had time to look at this for a while but still wanted to say thanks for looking into it and sending this solution. By the way, Jeff mentioned that the version of quantmod on the SVN (0.3.18) works for this. I tried to figure out how to download that version, but found the doc

[R] CRAN (and crantastic) updates this week

2012-04-08 Thread Crantastic
CRAN (and crantastic) updates this week New packages * bigdata (0.1) Maintainer: Han Liu Author(s): Han Liu, Tuo Zhao License: GPL-2 http://crantastic.org/packages/bigdata The big data package is a collection of scalable methods for large-scale data analysis. * boss (1.

Re: [R] assigment operator question

2012-04-08 Thread Thomas Lumley
On Sun, Apr 8, 2012 at 11:57 AM, Duncan Murdoch wrote: > On 12-04-07 4:51 PM, Henrik Bengtsson wrote: >> >> On Sat, Apr 7, 2012 at 1:30 PM, Mark Heckmann >>  wrote: >>> >>> Hello, >>> >>> using the<<- assignment operator I do not understand why the following >>> does not work. >>> >>> l<<- list()

[R] multithreaded low level functions

2012-04-08 Thread Benedict Holland
Hi all, I wanted to make sure that I am getting the most from my 8 core computer running R. I have a distance function like this: distance <- function(pointX, pointY, vecX, vecY) { #convert everything to radians pointX = pointX * pi/180 pointY = pointY * pi/180 vecX = vecX * pi/180 vecY

[R] Error in integrate(int.fn, lower = 0, upper = Inf) : evaluation of function gave a result of wrong length

2012-04-08 Thread Guaramy
Hi, i am writing a function to plot a pdf functions of a Generalized normal laplace distribution. The code is this { y = x-rho*mu cf.fn = function(s){ cplex = complex(1,0,1) temp1 = alpha*beta*exp(-sigma*s^2/2) temp2 = (alpha-

Re: [R] How to produce serveral plots with pairs of vectors

2012-04-08 Thread slyrs66
Thank you, David. What I actually would like to do in an more "elegant" way is the following boxplot (meinspss$attr_diff_gesamt ~ meinspss$R1_02) boxplot (meinspss$finanz_diff_gesamt ~ meinspss$R2_02) boxplot (meinspss$leist_diff_gesamt ~ meinspss$R3_02) boxplot (meinspss$soz_diff_gesamt ~ meinsp

[R] Difference between spec.pgram & spec.ar

2012-04-08 Thread Bazman76
Hi there, Can someone explain what the difference between spec.pgram and spec.ar is? I understand that they attempt to do the same thing one using an AR estimation of the underlying series to estimate teh sensity the other using the FFT. However when applied to teh same data set they seem to be g

Re: [R] axis labels not showing

2012-04-08 Thread Ivan Allaman
It's very simple my boy!! Do you already to play with "mar"? So..Try to change the values this object. For example, par(..., mar=c(*5*,2,2,2)). Bye! -- View this message in context: http://r.789695.n4.nabble.com/axis-labels-not-showing-tp4541268p4541354.html Sent from the R help mailing list

Re: [R] Error in xy.coords(x, NULL, log = log) : (list) object cannot be coerced to type 'double'

2012-04-08 Thread Bazman76
that's it! Thank you. -- View this message in context: http://r.789695.n4.nabble.com/Error-in-xy-coords-x-NULL-log-log-list-object-cannot-be-coerced-to-type-double-tp4524046p4541340.html Sent from the R help mailing list archive at Nabble.com. __ R-he

Re: [R] filling the matrix row by row in the order from lower to larger elements

2012-04-08 Thread Dimitri Liakhovitski
Agreed, Rui provided a very elegant vectorized solution - and I am very thankful to him. Unfortunately (for myself), I am not as proficient in vectorization - otherwise, I would not have asked the question. Why I did not follow on the original hint to use apply? For this reason (quote): > However -

Re: [R] Problems with CRAN package Ryacas

2012-04-08 Thread Gabor Grothendieck
On Sun, Apr 8, 2012 at 2:56 PM, Kjetil Halvorsen wrote: > First: >> sessionInfo() > R version 2.15.0 (2012-03-30) > Platform: x86_64-pc-linux-gnu (64-bit) > > locale: >  [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C >  [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8 >  [5] LC_MONETARY=en_US.UT

Re: [R] filling the matrix row by row in the order from lower to larger elements

2012-04-08 Thread ilai
On Sun, Apr 8, 2012 at 8:26 AM, Dimitri Liakhovitski wrote: > Sorry, I didn't have time to check the speed, indeed. > However - isn't apply the same as a loop, just hidden? > D. > Yes ?apply is a loop but not the same as ?for, see "Intro to R". As Bert Gunter pointed out the issue here was not sp

[R] Problems with CRAN package Ryacas

2012-04-08 Thread Kjetil Halvorsen
First: > sessionInfo() R version 2.15.0 (2012-03-30) Platform: x86_64-pc-linux-gnu (64-bit) locale: [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C [3] LC_TIME=en_US.UTF-8LC_COLLATE=en_US.UTF-8 [5] LC_MONETARY=en_US.UTF-8LC_MESSAGES=en_US.UTF-8 [7] LC_PAPER=C LC_NAME=C

Re: [R] axis labels not showing

2012-04-08 Thread Vikram Chhatre
Indeed. Changing margins brought back the axis labels. Thank you. What I would really like to do is stack these plots atop each other so they share the common x-axis. Is there a way to do it with ggplot2? Vikram On Sun, Apr 8, 2012 at 1:33 PM, ilai wrote: > On Sun, Apr 8, 2012 at 11:40 AM,

Re: [R] axis labels not showing

2012-04-08 Thread ilai
On Sun, Apr 8, 2012 at 11:40 AM, Vikram Chhatre wrote: > Hello - > > I want to generate stacked plots with par(mfrow)) function.  However, > my axis labels aren't showing. > Your mar (2) are too narrow. You could increase back to the default or use the lines option in mtext to write labels closer

[R] independent set function from graph package

2012-04-08 Thread Diogo Alagador
Dear all, I have been using the graph package until last week, especially the functions independence.number and largest.independent.vertex.sets. This weekend I formatted my hard disk and when trying to install packages I noticied that graph package is no longer on R repository although can

[R] axis labels not showing

2012-04-08 Thread Vikram Chhatre
Hello - I want to generate stacked plots with par(mfrow)) function. However, my axis labels aren't showing. My script is here: http://pastebin.com/yXXeMQgb The plot is here: http://www.crypticlineage.net/rdisc/strplot.pdf Thank you for your time. Vikram __

Re: [R] xyplot() does not plot legends with "relation=free" scales

2012-04-08 Thread ilai
On Sun, Apr 8, 2012 at 3:42 AM, Kaveh Vakili wrote: > Hi all, > > I have this problem with lattice that xyplot() won't draw some of my axis > labels if the type (i.e. the relation argument) of scales is set as free. For > example, in the plot below, I would want it to also show: > > 1. the label

Re: [R] How do Sweave users collaborate with Word users?

2012-04-08 Thread Alexander Shenkin
This sounds neat, and I wish you the best of luck with it. While I'm sure it'll be great for folks who are already curious enough to dive into Sweave or knitr, my guess is that coaxing collaborators into using markdown with SVN or GIT instead of passing around a Word file with Track Changes is goi

Re: [R] Make package out of own function

2012-04-08 Thread Uwe Ligges
Export the default method (and otehr methods) from your NAMESPACE? See Writing R Extensions about the S3method directive. Note that print.f is the print method for a class f, while f is a generic function, both fs are unrelated. Uwe Ligges On 02.04.2012 12:02, Johannes Radinger wrote: He

Re: [R] Error in xy.coords(x, NULL, log = log) : (list) object cannot be coerced to type 'double'

2012-04-08 Thread Uwe Ligges
On 01.04.2012 20:15, Bazman76 wrote: Hi there, When I run the code below I get the error Error in xy.coords(x, NULL, log = log) :(list) object cannot be coerced to type 'double' Any tips how I can resolve this? The mra output is a list, hence plot.ts cannot handle it. You probably want to

Re: [R] How to produce serveral plots with pairs of vectors

2012-04-08 Thread David Winsemius
On Apr 8, 2012, at 8:51 AM, slyrs66 wrote: Dear R-Users, I have a newbie-question about producing several plots with five variable pairs within in one code (loop). Problem: I define two objects, erveryone has five vectors (columns /variables): dimensionen <- data.frame(meinspss$attr_diff_g

Re: [R] Parallel writes in R

2012-04-08 Thread Uwe Ligges
On 02.04.2012 21:35, Jonathan Greenberg wrote: R-helpers: I'm curious what support R has for parallel writes to a binary file? If I want to use snow to have each node write different "rows" of a flat binary file (possibly out of sequence), are there any tricks/issues I should be aware of? Th

Re: [R] Drawing a line in xyplot

2012-04-08 Thread ilai
On Sat, Apr 7, 2012 at 8:29 PM, wcheckle wrote: > Thank you David, the bwplot option does what I need: > However, I am interested in also learning how to do it in xyplot as well.  I > wasn’t able to follow the last two set of instructions That was me. Sorry for any confusion. wcheckle, these we

Re: [R] filling the matrix row by row in the order from lower to larger elements

2012-04-08 Thread Bert Gunter
Preoccupation with speed of execution is typically (certainly not always) misplaced. Provided you have used sensible basic vectorization, loops in whatever form work adequately. First get working code. Then, if necessary, parallelization, byte compilation, or complex vectorization strategies can be

Re: [R] How to produce serveral plots with pairs of vectors

2012-04-08 Thread Bert Gunter
Please do your homework first: Read at least the beginning of An Introduction to R or other R tutorial -- there are many. I say this because in your post below you do not seem to realize that R is case sensitive, and you do not seem to know the difference between a data frame and a column thereof a

Re: [R] newbie question: strategy

2012-04-08 Thread sys ott
thanks for the heads up, I am checking all my providers EULAs to see which one I can share online... I imagine that as long as it is delayed 24 hours, there might be no issues... we will see... either way, worst case, I lost $9 for a registration... best case... I can share most ETF's with the comm

Re: [R] filling the matrix row by row in the order from lower to larger elements

2012-04-08 Thread Dimitri Liakhovitski
Thank you very much, Rui. This definitely produces the result needed. Again, I have not checked the speed yet. input<-as.matrix(data.frame(a=c(5,1,3,7),b=c(2,6,4,8))) f <- function(x){ nr <- nrow(x) result <- matrix(0, nrow=nr, ncol=ncol(x)) colnames(result) <- colnames(x)

Re: [R] filling the matrix row by row in the order from lower to larger elements

2012-04-08 Thread Dimitri Liakhovitski
Sorry, I didn't have time to check the speed, indeed. However - isn't apply the same as a loop, just hidden? D. On Fri, Apr 6, 2012 at 6:59 PM, ilai wrote: > On Fri, Apr 6, 2012 at 4:02 PM, Dimitri Liakhovitski > wrote: >  This works great: > > Really ? surprising given it is the EXACT same for-

Re: [R] .Rhistory in R.app

2012-04-08 Thread David Winsemius
On Apr 8, 2012, at 12:56 AM, Berend Hasselman wrote: On 08-04-2012, at 03:31, nkschmidt wrote: Hi Maria, you need to name your history file somefilename.Rapp.history, example savehistory(file = "history.Rapp.history") /nkschmidt You replied to a message almost two years old, which to m

[R] How to produce serveral plots with pairs of vectors

2012-04-08 Thread slyrs66
Dear R-Users, I have a newbie-question about producing several plots with five variable pairs within in one code (loop). Problem: I define two objects, erveryone has five vectors (columns /variables): dimensionen <- data.frame(meinspss$attr_diff_gesamt, meinspss$finanz_diff_gesamt, meinspss$leist

[R] Coloring comments differently

2012-04-08 Thread Renger van Nieuwkoop
Hi I use outline mode in ESS and was wondering if it is possible to change the color of the comments depending on the level. E.g. # Green Red ### Blue Cheers Renger Modelworks Gewerbestrasse 15 3600 Thun - Switzerland +41 79 818 53 73 i...@modelworks.c

Re: [R] xyplot() does not plot legends with "relation=free" scales

2012-04-08 Thread kv
Sorry, the B data.frame() didn't copy/paste well. Here is a working version (hopefully) B<-structure(list(yval = c(0.7, 0.61, 0.65, 0.63, 6.08, 0.64, 5.68,6.77, 1.48, 7.71, 0.82, 1.15, 0.54, 1.01, 0.59, 4.84, 0.69, 0.71, 8.7, 0.48, 0.69, 4.81, 1.42, 1.19, 0.84, 4.89, 0.85, 0.67, 7.07, 0.66, 7.9

[R] xyplot() does not plot legends with "relation=free" scales

2012-04-08 Thread Kaveh Vakili
Hi all, I have this problem with lattice that xyplot() won't draw some of my axis labels if the type (i.e. the relation argument) of scales is set as free. For example, in the plot below, I would want it to also show: 1. the labels E1,...E6 below the 10th panel (i.e. 3rd row, 2 col)just as

[R] Need help interpreting output from rcorrp.cens with Cox regression

2012-04-08 Thread Peter Jepsen
Dear R-listers, I am an MD and clinical epidemiologist developing a measure of comorbidity severity for patients with liver disease. Having developed my comorbidity score as the linear predictor from a Cox regression model I want to compare the discriminative ability of my comorbidity measure w

Re: [R] Avoid loop with the integrate function

2012-04-08 Thread Berend Hasselman
On 08-04-2012, at 08:28, Navin Goyal wrote: > Dear R users, > I am running a loop with the integrate function. I have pasted the code > below. I am integrating a function from time=0 to the time value in every > row. > I have to perform this integration over thousands of rows with different > par

Re: [R] newbie question: strategy

2012-04-08 Thread R. Michael Weylandt
Just a heads up: I'm pretty sure Josh and the other QS developers have access to intra-day data, but they can't provide an example using it in the package because the EULAs of most data providers won't allow direct redistribution of their data. The examples included all go online and download data