[R] Weighted Average application on Summary Dataset

Hi, I have 2 huge datasets - May and Jun - a miniscule sample of one is given below. I am trying to do 2 things with these datasets. I need to verify if the weighted average of variable A for a Reason in Jun is same/different from the same for May. To do this I am first computing the weighted ave

Re: [R] HOW to install RSQLite database

Yes i am asking how to install RSQLite packagein windows.Please help on this regards On 6/11/10, david.jessop [via R] < ml-node+2251498-1601505055-288...@n4.nabble.com> wrote: > > Are you asking how to install the RSQLite package or how to create a > SQLite database? The two are somewhat distinc

Re: [R] How to output text to sink from data frame line by line without column names

OK I see how to remove the line numbers[1] etc by using cat instead of print, but cannot work out how to remove the column names from the data frame output On 13/06/2010 4:21 PM, Nevil Amos wrote: I want to output a text file assembeld from various soruces within r ( actually as a genepop file

[R] How to output text to sink without initial line number [1], and data frame line by line without column names

I want to output a text file assembeld from various soruces within r ( actually as a genepop file) the output should be formatted line 1 "text comment" line 2:n selected column names from data frame

Re: [R] help with R

li li-13 wrote: > > Hi all, >I want to solve the following equation for x with rho <- 0.5 > >pnorm(-x)*pnorm((rho*dnorm(x)/pnorm(x)-x)/sqrt(1-rho^2))==0.05 > >Is there a function in R to do this? > Or if you wish to try different values of rho f <- function(x, rho) { pnorm(

Re: [R] Scope and sapply

I was careless. Here is a better example of what I am trying to. With the '<<-' you offered. ?<<- That was exactly what I needed, thankyou. cheers Worik N <- 10 ## x simulate a return series x <- runif(N)-.5 ## Build an array of cumulative returns of a portfolio starting with $1 as it chan

[R] using latticeExtra plotting confidence intervals

I am wanting to plot a 95% confidence band using segplot, yet I am wanting to have groups. For example if I have males and females, and then I have them in different races, I want the racial groups in different panels. I have this minor code, completely made up but gets at what I am wanting, 4 rand

Re: [R] Sim function

Hi: The book has an accompanying R package called arm. Within that package is a function named sim, which could be what you're looking for... HTH, Dennis On Sat, Jun 12, 2010 at 6:02 PM, Noah Silverman wrote: > I'm reading Gellman's book "Data Analysis Using Regression and > Multilevel-Hierarch

[R] Break in the y-axis

Hello all, I have been having trouble getting a break in my y-axis. All of my data points are up around 100-200, but the graph has to start at zero, so i would like to remove all the white space using a break symbol. I have been able to get the break and labels to be correct, however, I can't

Re: [R] mob (party package) question

Achim - Thank you so much for your suggestions. Tudor -- View this message in context: http://r.789695.n4.nabble.com/mob-party-package-question-tp2252500p2253141.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org m

Re: [R] Rgui crashed on Windows XP Home

On 2010-6-12 5:58, Duncan Murdoch wrote: Jinsong Zhao wrote: Hi there, I just installed R 2.11.1 on my PC, which runs a Windows XP Home. The installation is successful, however, when I double click on the R icon, I get the following error message: R for Windows GUI front-end has encountered a

Re: [R] Scope and sapply

On Sat, Jun 12, 2010 at 9:22 PM, Worik R wrote: > I am puzzled by the scope rules that apply with sapply. > > If I want to modify a vector with sapply I tried... > > N <- 10 > vec <- vector(mode="numeric", length=N) > test <- function(i){ >  vec[i] <- i > } > sapply(1:N, test) > vec > > but it not

Re: [R] mob (party package) question

Achim: Thanks a lot for your suggestions. Tudor > Date: Sun, 13 Jun 2010 00:26:08 +0200 > From: achim.zeil...@uibk.ac.at > To: tudor_bo...@hotmail.com > CC: r-help@r-project.org > Subject: Re: [R] mob (party package) question > > On Fri, 11 Jun 2010, tudor wrote: > > > > > Dear useRs

Re: [R] Scope and sapply

Hi, On Sat, Jun 12, 2010 at 9:22 PM, Worik R wrote: > I am puzzled by the scope rules that apply with sapply. > > If I want to modify a vector with sapply I tried... > > N <- 10 > vec <- vector(mode="numeric", length=N) > test <- function(i){ >  vec[i] <- i > } > sapply(1:N, test) > vec You didn

[R] Scope and sapply

I am puzzled by the scope rules that apply with sapply. If I want to modify a vector with sapply I tried... N <- 10 vec <- vector(mode="numeric", length=N) test <- function(i){ vec[i] <- i } sapply(1:N, test) vec but it not work. How can this be done? Worik [[alternative HTML versio

Re: [R] help with R

Dear Hannah, If I understand you correctly, you want a solution when rho is 0.5; if so, > f <- function(x){ + pnorm(-x)*pnorm((0.5*dnorm(x)/pnorm(x)-x)/sqrt(1-0.5^2)) - 0.05 + } > uniroot(f, c(-3, 3)) $root [1] 0.8031289 $f.root [1] -1.565857e-06 $iter [1] 11 $estim.prec [1] 6.103516e-05

[R] Sim function

I'm reading Gellman's book "Data Analysis Using Regression and Multilevel-Hierarchical Models" In Chapter 7 (and later), he makes frequent referent to a function names "sim". I can't find the function anywhere, not in my standard R install, or in any of the packages. Doe anyone have a refe

[R] help with R

Hi all, I want to solve the following equation for x with rho <- 0.5 pnorm(-x)*pnorm((rho*dnorm(x)/pnorm(x)-x)/sqrt(1-rho^2))==0.05 Is there a function in R to do this? Thank you very much! Hannah [[alternative HTML version deleted

Re: [R] Fast way to compute largest eigenvector

You can use the power method for computing the dominant eigenvector. A more sophisticated approach (for large matrices) is the Lancsoz algorithm for Hermitian matrices, which is based on the power method. The `arpack' function in the "igraph" package uses the more general Arnoldi iteration, wh

Re: [R] Can one get a list of recommended packages?

On 06/12/10 05:27 PM, Douglas Bates wrote: On Sat, Jun 12, 2010 at 8:37 AM, Dr. David Kirkby wrote: R 2.10.1 is used in the Sage maths project. Several recommended packages (Matrix, class, mgcv, nnet, rpart, spatial, and survival) are failing to build on Solaris 10 (SPARC). Have you checked

Re: [R] mob (party package) question

On Fri, 11 Jun 2010, tudor wrote: Dear useRs: I try to use mob from the party package (thanks Achim and Co.!) to model based recursive partition a data set. The model is a logistic regression specified with model=glinearModel and family=binomial(). Running mob results in a few warnings of th

Re: [R] Handling character string

^[[:space:]]* indicates that you want to match as many (zero or more) spaces at the beginning of the string. The '^' anchors the match at the beginning. [[:space:]]*$ says to match as many (zero or more that is what "*" asks for) at the end of the string; "$" anchors the search at the end of the

[R] Configuration of host address for database access

Hello everybody, Whatever I use as the host address of my PostgreSQL database server, I always get a connection if I'm accessing the database from R on the same machine and none if I'm trying to do so from a remote client. PostgreSQL is running on a windows xp machine as the RDBMS, JDBC, R, DBI

Re: [R] Clustering algorithms don't find obvious clusters

Henrik, Given your initial matrix, that should tell you which authors are similar/dissimilar to which other authors in terms of which authors they cite. In this case authors 1 and 3 are most similar because they both cite authors 2 and 4. Authors 2 and 3 are most different because they

[R] calling a function with new inputs every 1 minute

I have inputs to a function which are changing all the time - I pull these values from the internet. I then apply a function to the values. What I'd like to do is automate the process so it runs every one minute and adds the output of the function as a new element of a vector. Pseudo code: at

Re: [R] Clustering algorithms don't find obvious clusters

Dave, I used daisy with the default settings (daisy(M) where M is the matrix). Henrik On 11 June 2010 21:57, Dave Roberts wrote: > Henrik, > >The clustering algorithms you refer to (and almost all others) expect > the matrix to be symmetric. They do not seek a graph-theoretic solution, >

Re: [R] Handling character string

Thanks Jim for this reply. This is the way what I was looking for. However would you please explain me the meaning of ^[[:space:]]*" or '[[:space:]]+$'? When should I use "^" or "*" or "+$"? Thanks for your time. --- On Sat, 6/12/10, jim holtman wrote: From: jim holtman Subject: Re: [R] Handl

[R] who know how to program Bartlett lewis model using R

I am a graduate student working on disaggregation rainfall using Bartlett Lewis model but i have a problem on how can I make a program to estimate the model parameters and using it for the simulation. who can help of those people who are famous in R.   thanks for all in advance

Re: [R] Compiling R with multi-threaded BLAS math libraries - why not actually ?

On Sat, Jun 12, 2010 at 3:16 PM, Tal Galili wrote: > Hello Douglas, > > Thank you for the BLAST!=BLAS correction (I imagine my slip was due to some > working I have done recently with an RNA analysis software called BLAST). > > Also, thank you for the very interesting posting here and in your repl

Re: [R] Handling character string

This is probably what you want: > sub("^[[:space:]]*", "", " Now is the time") [1] "Now is the time" > You need to anchor it at the beginning with '^' On Sat, Jun 12, 2010 at 10:29 AM, Megh Dal wrote: > Thanks Erik for you reply. You have pointed correctly I want to remove the > "space" at

Re: [R] generating ordered, random decimal fractions

I must be missing something; can't you just use, e.g., one of: sort(round(runif(10), 3)) sort(sample(1000, 10, TRUE)/1000) ? -Peter Ehlers On 2010-06-12 7:49, kurt_h...@nps.gov wrote: Greetings How do I do this in R? Checking the Cran site produces a bewildering array of packages that

Re: [R] Transforming simulation data which is spread acrossmanyfiles into a barplot

Thanks Gabor - I was able to use that for my purposes. On 11 June 2010 16:27, Bert Gunter wrote: > So two time series? Fair enough. But less is more. Plot them as separates > series of points connected by lines, different colors for the two different > series. Or as two trellises plots. You may

[R] Fast way to compute largest eigenvector

Hello all, I was wondering if there is a function in R that only computes the eigenvector corresponding to the largest/smallest eigenvalue of an arbitrary real matrix. Thanks Minh -- Living on Earth may be expensive, but it includes an annual free trip around the Sun. ___

Re: [R] Can one get a list of recommended packages?

On 06/12/10 03:31 PM, Tal Galili wrote: Hello David, I am not sure I understood your question. Sorry, perhaps I should have rephrased it better. Are you asking what are the packages that the R release comes with? Sort of. When R is configured, there is an option --with-recommended-packa

Re: [R] Problems with BRugs

With BRugsFit() of the current verison of BRugs, the not given inits are not generated automatically (which may change with the next version). Hence you need to do it manually as in: library("BRugs") modelCheck("test.bug") modelData(bugsData(Data)) modelCompile(numChains=1) modelInits(bugsInit

Re: [R] Can one get a list of recommended packages?

On Sat, Jun 12, 2010 at 8:37 AM, Dr. David Kirkby wrote: > R 2.10.1 is used in the Sage maths project. Several recommended packages > (Matrix, class, mgcv, nnet, rpart, spatial, and survival) are failing to > build on Solaris 10 (SPARC). Have you checked the dependencies for those packages? Some

Re: [R] Logic with regexps

Thanks, Gabor, for the initiation to perly regexps! I've only been used to extended ones till now. A pity, perhaps, that "perl=TRUE" is not an option for the likes of browseEnv(), help.search(), list.files() and ls() (which take extended regexps), but one can always assign the output and then grep(

Re: [R] Handling character string

Thanks Erik for you reply. You have pointed correctly I want to remove the "space" at the 1st place (if any). In the mean time I have looked into the function sub() and there seems to be one example that mimics my problem : > str <- '   Now is the time      '> sub('[[:space:]]+$', '', str)[1] "  

Re: [R] Displaying "homogeneous groups" in aov post-hoc results ?

Thank you very much Hadley, exactly what I was looking for. Tal Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talgalili.com (Hebrew) | www.biostatistics.co.il (Hebrew) | www.r-statistics.com

[R] nonparametric density and probability methods

Hello, I tried to post this earlier, but it seems that it did not appear on the list. If you've rec'd 2 m I'm trying to calculate non-parametric probabilities using the np package and having some difficulties. OS is Windows, R version 2.11.1 Here is what I've done so far. library(np) veg <- da

Re: [R] Displaying "homogeneous groups" in aov post-hoc results ?

Try multcompView Hadley On Sat, Jun 12, 2010 at 8:42 AM, Tal Galili wrote: > Hello dear R-help mailing list, > > A friend of mine teaches a regression and experimental design course and > asked me the following question. > > She is trying to find a way to display the "homogeneous groups" (after

Re: [R] Can one get a list of recommended packages?

Hello David, I am not sure I understood your question. Are you asking what are the packages that the R release comes with? Or are you asking what recommended packages one should have when installing R? (There is a good list to start with here

Re: [R] Compiling R with multi-threaded BLAS math libraries - why not actually ?

Hello Douglas, Thank you for the BLAST!=BLAS correction (I imagine my slip was due to some working I have done recently with an RNA analysis software called BLAST). Also, thank you for the very interesting posting here and in your reply to David's post. My current conclusion from this thread are

Re: [R] Logic with regexps

On Sat, Jun 12, 2010 at 5:38 AM, Ted Harding wrote: > Greetings, > The following question has come up in an off-list discussion. > Is it possible to construct a regular expression 'rex' out of > two given regular expressions 'rex1' and 'rex2', such that a > character string X matches 'rex' if and

Re: [R] meta analysis with repeated measure-designs?

Dear Gerrit, the most appropriate approach for data of this type would be a proper multivariate meta-analytic model (along the lines of Kalaian & Raudenbush, 1996). Since you do not know the correlations of the reaction time measurements across conditions for the within-subject designs, a simpl

[R] generating ordered, random decimal fractions

Greetings How do I do this in R? Checking the Cran site produces a bewildering array of packages that I can't seem to find to load. Surely the main program has this function? Cheers Kurt *** Kurt Lewis Helf, Ph.D. Invertebrate Ecologist

[R] Displaying "homogeneous groups" in aov post-hoc results ?

Hello dear R-help mailing list, A friend of mine teaches a regression and experimental design course and asked me the following question. She is trying to find a way to display the "homogeneous groups" (after performing tukey test on an aov object). here's an example for what she means by "homog

Re: [R] Compiling R with multi-threaded BLAS math libraries - why not actually ?

On Sat, Jun 12, 2010 at 6:18 AM, Tal Galili wrote: > Hello Gabor, Matt, Dirk. > > Thank you all for clarifying the situation. > > So if I understand correctly then: > 1) Changing the BLAST would require specific BLAST per computer > configuration (OS/chipset). It's BLAS (Basic Linear Algebra Subr

[R] Can one get a list of recommended packages?

R 2.10.1 is used in the Sage maths project. Several recommended packages (Matrix, class, mgcv, nnet, rpart, spatial, and survival) are failing to build on Solaris 10 (SPARC). We would like to be able to get a list of the recommended packages for R 2.10.1, but ideally via a call to R, so it is

Re: [R] extended Kalman filter for survival data

Thanks Christos. BayesX was strongly recommended also by Fahrmeir and other people of that field. So it is clearly where I need to look at! Christophe iPhone.fan Le 12 juin 2010 à 04:46, Christos Argyropoulos a écrit : > > If you mean this paper by Fahrmeir: > http://biomet.oxfordjourn

Re: [R] Logic with regexps

Thanks, Brian. I had indeed overlooked grepl() (too busy delving into the syntax summary)! That is certainly a useful shortcut of the construction I had used. Your "Not in general" implies that using grep() twice (in this example; more times in more complex combinations) is inevitable -- which of

Re: [R] Logic with regexps

I think you have missed grepl(), e.g. X[grepl(rex1, X) & !grepl(rex2, X)] grepl is a fairly recent addition (2.9.0) that is used extensively in R's own text-processing operations (e.g. help files, utilities such as 'R CMD check'). On Sat, 12 Jun 2010, ted.hard...@manchester.ac.uk wrote: Gr

[R] help with npcdens function package np

Hello, I'm trying to calculate non-parametric probabilities using the np package and having some difficulties. OS is Windows, R version 2.11.1 Here is what I've done so far. library(np) veg <- data.frame(factor(Physiogomy), meanAnnualDepthAve, TP) attach(veg) : for clarification dim(veg) retur

Re: [R] points marking

Hi, Well Thanks for letting me know that pch is of no use with segments petr. I am using lend as it suits to me more as gregory suggested , but I am not getting imite??? think I try to fix it with some other method also, as I have to deal more with the symbols in this case, But I want to the know

Re: [R] Compiling R with multi-threaded BLAS math libraries - why not actually ?

Hello Gabor, Matt, Dirk. Thank you all for clarifying the situation. So if I understand correctly then: 1) Changing the BLAST would require specific BLAST per computer configuration (OS/chipset). 2) The advantage would be available only when doing _lots_ of linear algebra So I am left wonderin

[R] meta analysis with repeated measure-designs?

Dear all, I am trying to run a meta analysis of psycholinguistic reaction-time experiments with the meta package. The problem is that most of the studies have a within-subject designs and use repeated measures ANOVAs to analyze their data. So at present it seems that there are three non-optimal

[R] Logic with regexps

Greetings, The following question has come up in an off-list discussion. Is it possible to construct a regular expression 'rex' out of two given regular expressions 'rex1' and 'rex2', such that a character string X matches 'rex' if and only if X matches 'rex1' AND X does not match 'rex2'? The desi

[R] mob (party package) question

Dear useRs: I try to use mob from the party package (thanks Achim and Co.!) to model based recursive partition a data set. The model is a logistic regression specified with model=glinearModel and family=binomial(). Running mob results in a few warnings of the type: In glm.fit ... algorithm did

Re: [R] extended Kalman filter for survival data

If you mean this paper by Fahrmeir: http://biomet.oxfordjournals.org/cgi/content/abstract/81/2/317 I would recommend BayesX: http://www.stat.uni-muenchen.de/~bayesx/. BayesX interfaces with R and estimates discrete (and continuous) time survival data with penalized regression methods. If

[R] sleep timer resolution on OSX

I am doing some timing experiments. To test looping performance, I used the Sys.sleep function. I noticed something in the docs that is just a little misleading: " The resolution of the time interval is system-dependent, but will normally be down to 0.02 secs or better. (On modern Unix