[R] Big data and column correspondence problem

2011-07-25 Thread murilofm
Greetings, I've been struggling for some time with a problem concerning a big database that i have to deal with. I'll try to exemplify my problem since the database is really big. Suppose I have the following data: AA = c(4,4,4,2,2,6,8,9) A1 = c(3,3,5,5,5,7,11,12) A2 = c(3,3,5,5,5,7,11,12) A = cb

Re: [R] Trouble with line of best fit

2011-07-25 Thread DimmestLemming
Daniel Malter wrote: > > Second, given the figure, a linear specification is obviously a > misspecification of your model, unless you account for autocorrelation. > I've decided to use this as a learning opportunity. I looked up autocorrelation: It does not apply in any way. By your standards,

Re: [R] grey colored lines and overwriting labels i qqplot2

2011-07-25 Thread Sigrid
Thank you Brian. Sorry for being such a noob. I am not a programmer and just learning R by myself. This is was I typed, but ended up with a couple error messages. > df <-structure(list(year = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, + 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,

Re: [R] Recoding Multiple Variables in a Data Frame in One Step

2011-07-25 Thread Peter Ehlers
On 2011-07-25 15:48, William Dunlap wrote: -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of David Winsemius Sent: Monday, July 25, 2011 3:39 PM To: Anthony Damico Cc: r-help@r-project.org Subject: Re: [R] Recoding Multiple Variables

Re: [R] Selecting unique values

2011-07-25 Thread Steven Kennedy
use unique() df <- data.frame(Species_name=c(rep('Abies concolor',4),rep('Accipiter cooperi',7)), Longitude=c(-106.601, -106.493, -106.489, -106.496, -119.688, -119.792, -118.797,-77.38333,-77.38333,-75.99153,-75.99153), Latitude=c(35.868, 35.9682, 35.892, 35.85

Re: [R] Accessing the index of factor in by() function

2011-07-25 Thread Ista Zahn
Hi Merik, by() works most easily with data.frames. Is this what you are after? my.plot <- function(dat) { print(dat$value); print(dat$month[dat$id==dat$value]) } by(dat.tmp, id, my.plot) Best, Ista On Mon, Jul 25, 2011 at 9:19 PM, Merik Nanish wrote: > Hello, > > Here are three vectors to give

Re: [R] library 'ts' not available?

2011-07-25 Thread cplusplus programmer
guess you're rightI need to look at more recent code examplesthx! On Mon, Jul 25, 2011 at 11:01 PM, Dennis Murphy wrote: > I don't recall exactly (ca 2004?) when the base R package structure > was reorganized, but some time ago a number of packages were merged > into a much smaller subse

[R] Creating png of layered legend

2011-07-25 Thread Amelia McNamara
I am trying to create a plot that has multiple plot characters for each point (e.g. a point within a triangle, a triangle within a square, etc). The workaround I have found to do this is by plotting twice, as in this example: x <- c(1.1, 2.3, 4.6) y <- c(2.0, 1.6, 3.2) plot(x, y) points(x,y, pch=2

Re: [R] Trouble with line of best fit

2011-07-25 Thread DimmestLemming
Thanks! It was easy to solve, then... I just hadn't considered that. For the record: I've never actually taken statistics, so I'm learning it at the same time as R. So yes, this is a private project. -- View this message in context: http://r.789695.n4.nabble.com/Trouble-with-line-of-best-fit-tp3

[R] Selecting unique values

2011-07-25 Thread Kumar Mainali
Greetings I have a dataset with occurrence records of multiple species. I need to get rid of multiple listings of the same occurrence point for a species (as you see below in red and blue typeface). How do I create a dataset only with unique set of longitude and latitude for each species? Thanks i

[R] trouble with uninstall

2011-07-25 Thread John Hangen
Hello, Today I installed R for Windows version 2.13.1, but realized that it's not something that I need right now.I attempted to uninstall it, but am getting the following error message: "This Installation can only be uninstalled on 64-bit Windows." I am running Windows XP Professional, Ver

[R] Dealing with many parameters in a function

2011-07-25 Thread Abraham Mathew
I'm creating a function in R. However, I have a large number of function parameters, and need to find an efficient solution for running the function with all the parameters. So in the following function, I have about 20 parameters that I assign to the function, with almost all the values being diff

[R] lme convergence error

2011-07-25 Thread Ben Grannan
Hello, I am working from a linux 64 machine on a server with R-2.12 (I can't update to 2.13). I am iterating through many linear mixed models for longitudinal data and I occasionally receive the following convergence error: > BI.lme <- lme(cd4 ~ time + genBI + genBI:time + C1 + C2 + C11 + C12, ran

[R] error in survival analysis

2011-07-25 Thread Rasika Gawde
This is a simple R program that I have been trying to run. I keep running into the "singular matrix" error. I end up with no sensible results. Can anyone suggest any changes or a way around this? I am a total rookie when working with R. Thanks, Rasika > library(survival) Loading required packa

[R] (no subject)

2011-07-25 Thread John Hangen
John W. Hangen, GISP GIS Coordinator City of Milford, CT * Please consider the environment before printing this email [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PL

[R] Is there an R program that produces optimal solution/mix of multiple samples' varying volumes and values

2011-07-25 Thread luke
Sorry about the lengthy subject line. Anyone know of an R' program that can look at several sources' varying available volumes/amounts and their individual set of values, compared to a "target" range/curve for these values, to find the optimal mixture(s) of these possible sources for the desired cu

[R] Rescaling columns in a multi-plot layout

2011-07-25 Thread Manojit Roy
Dear all, I am trying to create a 6-plot layout - 3 rows and 2 columns - so that only the top two plots have variable widths, all else with their default setting. Using "layout(matrix(c(1,2,3,4,5,6),3,2,byrow=T),widths=c(5,2))" rescales column #2 of all three rows, whereas I would like to res

[R] Debugging multiple imputation in mice

2011-07-25 Thread Harish Narayanan
Hello all, I am trying to impute some missing data using the mice package. The data set I am working with contains 125 variables (190 observations), involving both categorical and continuous data. Some of these variables are missing up to 30% of their data. I am running into a peculiar problem wh

Re: [R] squared "pie chart" - is there such a thing?

2011-07-25 Thread Naomi Robbins
Andy, Mosaic plots were suggested by Sarah Goslee on July 21st. Unfortunately, I broke the chain of messages when I was away from home, in a rush and didn't know how to respond to messages from a digest. Naomi -- Naomi B. Robbins NBR 11 Christine Court Wayne, NJ 07470 Phone: (973) 694-6009

[R] Trouble with line of best fit

2011-07-25 Thread DimmestLemming
I don't usually do much with graphs in R, and this is my first time adding a line of best fit. Hopefully this is an easy problem to solve. I'm looking at a variable called soloKills along the range 5:28. Here are all my commands, in script form: range=5:28 graph=soloKills title="Solo kill/death d

[R] Life Cycle Assessment with R.

2011-07-25 Thread Jose Bustos Melo
Hello everyone,  There's something really important about climate change and how many institutions around the globe are looking for softwares solutions in order to achieve they (and everyone) needs to improve life conditions in all the planet. Currently, they're many comercial softwares working w

Re: [R] Fwd: Rgraphviz installation problem

2011-07-25 Thread Ram H. Sharma
Hi Martin Thank you for help, that helped me to come out the problem ! On Mon, Jul 25, 2011 at 11:46 AM, Martin Morgan wrote: > On 07/25/2011 07:58 AM, Ram H. Sharma wrote: > >> Hi All >> >> How can I adjust this path. My current path of Graphviz installation is: >> >> C:\Program Files (x86)\G

[R] Installing CAIC

2011-07-25 Thread eilunedpearce
Hi, I'm trying to install CAIC directly into the newest version of R using the code on the R-Forge CAIC website and I get an error message: install.packages("CAIC", repos="http://R-Forge.R-project.org";) Warning message: In getDependencies(pkgs, dependencies, available, lib) : package ‘CAIC’ i

[R] Accessing the index of factor in by() function

2011-07-25 Thread Merik Nanish
Hello, Here are three vectors to give context to my question below: *id<- c(1,1,1,1,1,2,2,2,3,3,3)) month <- c(1, 1, 2, 3, 6, 2, 3, 6, 1, 3, 5) value <- c(10, 12, 11, 14, 16, 12, 10, 8, 14, 11, 15)* and I want to plot "value" over "month" separately for each "id". Before I can do that, I nee

Re: [R] library 'ts' not available?

2011-07-25 Thread Dennis Murphy
I don't recall exactly (ca 2004?) when the base R package structure was reorganized, but some time ago a number of packages were merged into a much smaller subset. Several of these were merged into the stats package, including (but not exclusive to) ts, mva, nls and ctest. Code using such packages

Re: [R] Recoding Multiple Variables in a Data Frame in One Step

2011-07-25 Thread David Winsemius
On Jul 25, 2011, at 6:48 PM, William Dunlap wrote: -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org ] On Behalf Of David Winsemius Sent: Monday, July 25, 2011 3:39 PM To: Anthony Damico Cc: r-help@r-project.org Subject: Re: [R] Recoding Multip

[R] Package or procedure recommendations for analysis of repeated cross-sections?

2011-07-25 Thread andrewH
I have a survey data set of 6 years and about 1500 persons surveyed per year, with roughly 200 questions per survey. The samples are drawn independently without replacement and are intended to represent the nation (USA). I would like to create something like a synthetic panel, dividing the respond

Re: [R] Simple example of using a closure in R to manage bank accounts?

2011-07-25 Thread Michael Hannon
Yes, I think that's the same example that I saw previously. Thanks, Gabor. -- Mike - Original Message > From: Gabor Grothendieck > To: Michael Hannon > Cc: R Help > Sent: Mon, July 25, 2011 3:40:58 PM > Subject: Re: [R] Simple example of using a closure in R to manage bank >accou

Re: [R] Simple example of using a closure in R to manage bank accounts?

2011-07-25 Thread Michael Hannon
Hi, Leo. Yes, that's exactly the idea. Thanks, -- Mike > >From: Leo Guelman >To: Michael Hannon >Cc: R Help >Sent: Mon, July 25, 2011 3:16:13 PM >Subject: Re: [R] Simple example of using a closure in R to manage bank accounts? > > >something like this maybe? > >account <- function (bala

Re: [R] convert TS dataframe to evenly spaced intervals?

2011-07-25 Thread David Winsemius
On Jul 22, 2011, at 7:01 PM, Derrick Lin wrote: Hi R-help, I have a dataframe consisting of a time-series [t, v]. The timestamps aren't at all evenly spaced. The values are continuous. I've been able to graph this as a step function (which is what it should be) in ggplot2, using the 'ste

Re: [R] library 'ts' not available?

2011-07-25 Thread David Winsemius
On Jul 22, 2011, at 5:03 PM, cplusplus programmer wrote: hi, When I type library(ts) (I am working on an ubuntu machine), I get the following error: library(ts) Error in library(ts) : there is no package called 'ts' You would have needed to have earlier installed that package before tryi

Re: [R] standard error of exp(coef) from coxph

2011-07-25 Thread David Winsemius
On Jul 23, 2011, at 1:41 AM, Ehsan Karim wrote: Dear List, Must be a silly question, but I was wondering whether there is a direct way of calculating "standard error of a HR or exp(coef)" from coxph objects x <- coxph(Surv(time, status) ~ age + inst, lung)> xcoef exp(coef) se(coe

Re: [R] Find pattern in matrix

2011-07-25 Thread David Winsemius
On Jul 21, 2011, at 6:12 PM, Q wrote: Hello all, I'm trying to find all combinations of 4 numbers that satisfy 4 criteria, inside of a matrix (62 x 25). I've found a way to do this using for loops, but it is extremely slow because it involves checking every possible combination of number

Re: [R] two sample histogram

2011-07-25 Thread Joshua Wiley
Hi, Just a few general notes. First, please read the posting guide. It requests that you: A) state your version of R (recommended by providing the results of sessionInfo() ) B) provide a minimal, reproducible example C) post in plain text not HTML I realize that you are new, but following the

Re: [R] Recoding Multiple Variables in a Data Frame in One Step

2011-07-25 Thread William Dunlap
> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On > Behalf Of William Dunlap > Sent: Monday, July 25, 2011 3:49 PM > To: David Winsemius; Anthony Damico > Cc: r-help@r-project.org > Subject: Re: [R] Recoding Multiple Variables in a Data Fram

Re: [R] Recoding Multiple Variables in a Data Frame in One Step

2011-07-25 Thread William Dunlap
> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On > Behalf Of David Winsemius > Sent: Monday, July 25, 2011 3:39 PM > To: Anthony Damico > Cc: r-help@r-project.org > Subject: Re: [R] Recoding Multiple Variables in a Data Frame in One Step >

Re: [R] Simple example of using a closure in R to manage bank accounts?

2011-07-25 Thread Gabor Grothendieck
On Mon, Jul 25, 2011 at 6:00 PM, Michael Hannon wrote: > Greetings.  I once ran across a simple (toy) example of using a closure in R > to > manage bank accounts.  I've got a use for it now but can no longer find it.   > If > you have it (or a similar example), will you please send it to me? > (U

Re: [R] Recoding Multiple Variables in a Data Frame in One Step

2011-07-25 Thread David Winsemius
On Jul 21, 2011, at 8:06 PM, Anthony Damico wrote: Hi, I can't for the life of me find how to do this in base R, but I'd be surprised if it's not possible. I'm just trying to replace multiple columns at once in a data frame. #load example data data(api) #this displays the three columns and

Re: [R] Simple example of using a closure in R to manage bank accounts?

2011-07-25 Thread Leo Guelman
something like this maybe? account <- function (balance=0) { function (d = 0,w = 0) { newbal <- balance + d - w balance <<- newbal newbal } } John <- account(100) John John(d=100,w=50) John() Leo <- account(1000) Leo Leo(d=1000,w=50) Leo() Leo(d=100,w

Re: [R] Help with 3-d barplot in R

2011-07-25 Thread Rolf Turner
Right on, Red Freak!!! :-) cheers, Rolf On 25/07/11 20:23, Joshua Wiley wrote: On Mon, Jul 25, 2011 at 12:49 AM, ATANU wrote: i am trying to make 3-d barplots,pie-charts in R,just like Excel. i have used rgl , but that does not produce beautiful graphs like excel(i dont need to

[R] Simple example of using a closure in R to manage bank accounts?

2011-07-25 Thread Michael Hannon
Greetings. I once ran across a simple (toy) example of using a closure in R to manage bank accounts. I've got a use for it now but can no longer find it. If you have it (or a similar example), will you please send it to me? (Unfortunately, a web search that includes the terms "bank" and "clo

Re: [R] Binding multiple data frames into single data frame

2011-07-25 Thread Joshua Wiley
This and earlier posts by you make me wonder if perhaps you are doing some data management in R that would be better handled by software explicitly designed for it limitting the import into R for higher level analyses/graphing/whatever. In any case, you can easily break this big task into smaller

Re: [R] biglm() and NeweyWest()

2011-07-25 Thread Thomas Lumley
On Tue, Jul 26, 2011 at 12:29 AM, Philipp Grueber wrote: > Dear all, > > I am working on a large dataset and need to use biglm() to perform OLS > regressions. I have detected significant ARCH effects which I try to account > for using the Newey-West correction. > > So far, I have worked with Newey

Re: [R] ggplot question: changing the label for the Y axis on a histogram

2011-07-25 Thread Steven Rytina, Prof.
Thanks, upgrading the software did the trick. -Original Message- From: istaz...@gmail.com [mailto:istaz...@gmail.com] On Behalf Of Ista Zahn Sent: Monday, July 25, 2011 3:50 PM To: Steven Rytina, Prof. Cc: r-help@r-project.org Subject: Re: [R] ggplot question: changing the label for th

Re: [R] Binding multiple data frames into single data frame

2011-07-25 Thread Madana_Babu
This is working for me... however when i use multicore there are more than 60K DF's getting created and when i use do.call() function it is taking huge time. Is there any function which can perform this operation at faster rate? Regards, Madana -- View this message in context: http://r.789695.n4

Re: [R] Binding multiple data frames into single data frame

2011-07-25 Thread Madana_Babu
Hi Joshua, This worked for me. Thanks for your help. Now i have another challenge. Since multicore is creating almost 70K DF's, do.call() function taking huge amount of time to bind all the datasets and i am not able to realize the impact of multicore due to do.call() time consumption. Can you

Re: [R] ggplot question: changing the label for the Y axis on a histogram

2011-07-25 Thread Ista Zahn
Hi Steven, It works as expected for me (see session info below). What version of R and ggplot are you using? Best, Ista sessionInfo() R version 2.13.0 (2011-04-13) Platform: i386-pc-mingw32/i386 (32-bit) locale: [1] LC_COLLATE=English_United States.1252 [2] LC_CTYPE=English_United States.1252 [3

Re: [R] Trouble with line of best fit

2011-07-25 Thread Joshua Wiley
On Jul 25, 2011, at 12:31, Daniel Malter wrote: > Hi, > > I hope this is for private purposes. Otherwise, I may cite David Winsemius: > "The advancement of science would be safer if you knew what you were doing." Hmm, even if it were not private, the science of FPS video games is a dangerous p

Re: [R] Standardizing the number of records by group

2011-07-25 Thread Peter Ehlers
On 2011-07-25 12:24, Sam Albers wrote: Hello R-help, I have some data collected at regular intervals but for a varying length of time. I would like to standardize the length of time collected and I can do this by standardizing the number of records I use for my analysis. Take for example the da

Re: [R] Finding/identifying a value within a factor

2011-07-25 Thread Peter Ehlers
On 2011-07-25 10:51, "Dénes TÓTH" wrote: Hi, you provided a character vector as an example. I guess you meant something like: x<- factor(c("1","2","3","4","<1")) # You can identify those elements with an "<" by ?grep or ?grepl: indices<- grep("<",as.character(x)) # You can transform those ele

Re: [R] Trouble with line of best fit

2011-07-25 Thread Daniel Malter
Hi, I hope this is for private purposes. Otherwise, I may cite David Winsemius: "The advancement of science would be safer if you knew what you were doing." First, your regression command is inverted. You ought to regress SoloKills on range, not vice versa. abline(lm(graph~range)) #does the tric

Re: [R] ggplot question: changing the label for the Y axis on a histogram

2011-07-25 Thread Steven Rytina, Prof.
Thanks for taking a look. As you can see from my (reduced) example below, I am drawing on counts and then asking bin to aggregate. My hunch is that ylab() somehow gets over-ridden in the collaboration implicit in geom_histogram. Accordingly, any comment that illuminates

[R] Standardizing the number of records by group

2011-07-25 Thread Sam Albers
Hello R-help, I have some data collected at regular intervals but for a varying length of time. I would like to standardize the length of time collected and I can do this by standardizing the number of records I use for my analysis. Take for example the data set below: library(plyr) x <- runif(

Re: [R] ARIMA simulations

2011-07-25 Thread Don McKenzie
?arima.sim in the stats package included with the standard distribution of R. (gurus -- is this considered state of the art? thx) On Jul 25, 2011, at 12:14 PM, Luis Felipe Parra wrote: > Hello, > > I have estimated an ARIMA model and I would like to make simulations from > this estimate

[R] ARIMA simulations

2011-07-25 Thread Luis Felipe Parra
Hello, I have estimated an ARIMA model and I would like to make simulations from this estimated model 1,5 and 10 steps ahead. Does anybody know how to do this? Thank You Felipe Parra [[alternative HTML version deleted]] __ R-help@r-project.or

Re: [R] Variable scope in functions - best practices

2011-07-25 Thread Dieter Menne
Noah Silverman wrote: > > > I have several "global" variables that I want to change with a given > function. (The variable has a different value after the function is > called.) > > Berend Hasselman wrote: > > > Maybe this helps > > ?`<<-` > It helps to get the job done, but the OP aske

Re: [R] Problem with random number simulation

2011-07-25 Thread Nordlund, Dan (DSHS/RDA)
> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- > project.org] On Behalf Of Nordlund, Dan (DSHS/RDA) > Sent: Monday, July 25, 2011 10:39 AM > To: r-help@r-project.org > Subject: Re: [R] Problem with random number simulation > <<>> > > I don't use Nabble

Re: [R] [R-pkgs] rdyncall 0.7.3

2011-07-25 Thread Dieter Menne
Daniel Adler wrote: > > Initial Announcement: Package rdyncall released on CRAN. (Version 0.7.3) > > The package was presented at the Use!R 2009 with the title > 'An improved Foreign Function Interface for R' and is now available on > CRAN > and considered stable for a large range of R platfor

Re: [R] Extrat rows from a dataframe based on second data frame

2011-07-25 Thread Dieter Menne
Bansal, Vikas wrote: > > > I have two data frames > > df > > ChrPos CaseA CaseCCaseG CaseT > 10 135349878 0.00 3.428571 0.00 20.571429 > 10 135349880 0.00 21.33 0.00 2.67 > 10 135349883 21.00 0.00 3.00 0.000

Re: [R] ggplot question: changing the label for the Y axis on a histogram

2011-07-25 Thread Ista Zahn
Hi Steven, The basic example works: qplot(rnorm(100)) + ylab("Bin Counts") Please post a reproducible example that illustrates the problem you are having. Best, Ista On Mon, Jul 25, 2011 at 1:40 PM, Steven Rytina, Prof. wrote: > Some help with how to re-label the vertical axis in a histogram  

Re: [R] Finding/identifying a value within a factor

2011-07-25 Thread Dénes TÓTH
Hi, you provided a character vector as an example. I guess you meant something like: x <- factor(c("1","2","3","4","<1")) # You can identify those elements with an "<" by ?grep or ?grepl: indices <- grep("<",as.character(x)) # You can transform those elements by ?as.numeric as.numeric(x[indices

Re: [R] Finding/identifying a value within a factor

2011-07-25 Thread Daniel Malter
b<-c("1","2","3","4","<1") grep('<',b) HTH, Daniel Ryan Utz-2 wrote: > > Hi all, > > I'm trying to identify a particular digit or value within a vector of > factors. Specifically, this is environmental data where in some cases the > minimum value reported is "<" a particular number (and I wan

Re: [R] Finding/identifying a value within a factor

2011-07-25 Thread Sarah Goslee
Hi Ryan, The key to this is the grep() command. It's easiest if you import your data as character rather than factor, since these aren't properly factors. You don't say how that's done, but if you are using read.table() then as.is=TRUE will prevent conversion to factor. For a character vector, h

[R] ggplot question: changing the label for the Y axis on a histogram

2011-07-25 Thread Steven Rytina, Prof.
Some help with how to re-label the vertical axis in a histogram would be appreciated. qplot(off.sc,weight=rel.freq,binwidth=.29,main="test Figure"+ylab("New from inside"))+ylab("New from outside")+ xlab("off.sc\nAg

Re: [R] Problem with random number simulation

2011-07-25 Thread Nordlund, Dan (DSHS/RDA)
> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- > project.org] On Behalf Of Petr PIKAL > Sent: Monday, July 25, 2011 9:05 AM > To: AriGold > Cc: r-help@r-project.org > Subject: Re: [R] Problem with random number simulation > > Hi > > r-help-boun...@r-pr

Re: [R] Plotting compound functions--help with defining x-axis as f(x)

2011-07-25 Thread Uwe Ligges
On 25.07.2011 02:34, Aimee Jones wrote: Thanks for your assistance, this worked perfectly, and sorry for posting in html! I shall check it's plain text for future postings.. Would it be possible for you to explain why this works? I'm unsure as to why redefining t works, when t is perfectly def

[R] Finding/identifying a value within a factor

2011-07-25 Thread Ryan Utz
Hi all, I'm trying to identify a particular digit or value within a vector of factors. Specifically, this is environmental data where in some cases the minimum value reported is "<" a particular number (and I want to manipulate only these). For example: x<-c("1","2","3","4","<1") For a dataset

Re: [R] How to find the likelihood of a null model in R

2011-07-25 Thread Rubén Roa
I know how to do what you want. However, the fact that you didn't bother to follow the rules of the list eliminates the initial, admittedly faint, enthusiasm to help you. Maybe someone else will guide you anyways. The rules are extremely simple. Read them at the end of this message. If you do yo

[R] Odp: How to merge a matrix and a dataframe with different types of columns

2011-07-25 Thread Petr PIKAL
Hi > [R] How to merge a matrix and a dataframe with different types of columns > > Hi, > > > > I hope someone can help here. I have been scratching my head for many > hours, trying to find out how to merge two datasets that were created by > different means. The initial miniMarket contains

Re: [R] Wide confidence intervals or Error message in a mixed effects model (nlme)

2011-07-25 Thread Ben Bolker
Dieter Menne menne-biomed.de> writes: [snip] > I always thought that that list was about > lme4/lmer development, it still says > > "notably lmer() related" > > Does Douglas Bate really want questions of this type, and on nlme/lme, on > that list? For what it's worth, in practice the list

Re: [R] Wide confidence intervals or Error message in a mixed effects model (nlme)

2011-07-25 Thread Dieter Menne
Bert Gunter wrote: > > I would also address future questions of this sort to the > R-sig-mixed-models list, as this is not really an r-help kind of question. > > That's a bit confusing, Bert. I always thought that that list was about lme4/lmer development, it still says "notably lmer() relate

Re: [R] Problem with random number simulation

2011-07-25 Thread Petr PIKAL
Hi r-help-boun...@r-project.org napsal dne 25.07.2011 05:23:51: > Re: [R] Problem with random number simulation > > Hi, > > I was just hoping someone could try running the program to see if they get > any output. This is independent research. I just need to run this > simulation to see how a

Re: [R] Lattice: distance of tick labels from axis line

2011-07-25 Thread Peter Ehlers
On 2011-07-25 07:54, marcel wrote: Hi, I am doing fine editing on a lattice plot, now have the Y-axis title in the correct position, but the tick labels are too far from the axis line. I looked at the help documentation but could not find how to change this. This seems to be easy to do in basic p

[R] Extrat rows from a dataframe based on second data frame

2011-07-25 Thread Bansal, Vikas
Dear all, I have also attached my question in txt file because of bad column spacing. I have two data frames df ChrPos CaseA CaseCCaseG CaseT 10 135349878 0.00 3.428571 0.00 20.571429 10 135349880 0.00 21.33 0.00 2.67 10 13

Re: [R] Fwd: Rgraphviz installation problem

2011-07-25 Thread Martin Morgan
On 07/25/2011 07:58 AM, Ram H. Sharma wrote: Hi All How can I adjust this path. My current path of Graphviz installation is: C:\Program Files (x86)\Graphviz2.20\bin and Rgraphviz is C:\R\library\Rgraphviz\libs\i386 Please help me. I am new to Rgraphviz. PATH is an environment variable. Fro

Re: [R] Binning numbers into integer-valued intervals (or: a version of cut or cut2 that makes sense)

2011-07-25 Thread William Dunlap
> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On > Behalf Of Karl Ove > Hufthammer > Sent: Monday, July 25, 2011 7:01 AM > To: r-h...@stat.math.ethz.ch > Subject: [R] Binning numbers into integer-valued intervals (or: a version of > cut or

Re: [R] Wide confidence intervals or Error message in a mixed effects model (nlme)

2011-07-25 Thread Bert Gunter
Inline below. -- Bert On Mon, Jul 25, 2011 at 7:42 AM, Dieter Menne wrote: > > Menelaos Stavrinides-2 wrote: >> >> I am analyzing a dataset on the effects of six pesticides on population >> growth rate of a predatory mite. The response variable is the population >> growth rate of the mite (ranges

Re: [R] squared "pie chart" - is there such a thing?

2011-07-25 Thread Liaw, Andy
Has anyone suggested mosaic displays? That's the closest I can think of as a "square pie chart"... > -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Naomi Robbins > Sent: Sunday, July 24, 2011 7:09 AM > To: Thomas Levine > Cc

[R] Lattice: distance of tick labels from axis line

2011-07-25 Thread marcel
Hi, I am doing fine editing on a lattice plot, now have the Y-axis title in the correct position, but the tick labels are too far from the axis line. I looked at the help documentation but could not find how to change this. This seems to be easy to do in basic plots with "at=c(value1, value2, value

[R] LogLik warning() corresponding probabilities for node ... were 0 or NaN.

2011-07-25 Thread Liam Woulfe
Hello R Helpers, I am trying to find the log liklihood between the bn.fit of one set of data and another. At small sample sizes I get the warning message In FUN(c("...", "...", "...", "...", ... : 1 observations were dropped because the corresponding probabilities for node ... were 0 o

[R] Fwd: Rgraphviz installation problem

2011-07-25 Thread Ram H. Sharma
Hi All How can I adjust this path. My current path of Graphviz installation is: C:\Program Files (x86)\Graphviz2.20\bin and Rgraphviz is C:\R\library\Rgraphviz\libs\i386 Please help me. I am new to Rgraphviz. Ram "your Windows system 'PATH' variable needs to be adjusted once you've instal

[R] error in optimization when I include constant term in Klein and Spady (np package)

2011-07-25 Thread Dimitris.Kapetanakis
Dear all, I am trying to use the np package for the estimation of a model with Klein and Spady's semiparametric estimator. When though, I include a constant term ( a column with 1s in X) then the following message appear: Multistart 1 of 3...Error in optim(c(beta, h), fn = optim.fn, gr = NULL, m

Re: [R] compute the mean of two (or more) correlations

2011-07-25 Thread Joshua Wiley
Hi Martin, Off hand I do not know of a pre-defined function to do it, but the "z-transformation" is just the inverse hyperbolic function, the mean is just the mean, and the back transformation is the hypoerbolic function so... x <- c(.5, .4) tanh(mean(atanh(x))) should

Re: [R] Wide confidence intervals or Error message in a mixed effects model (nlme)

2011-07-25 Thread Dieter Menne
Menelaos Stavrinides-2 wrote: > > I am analyzing a dataset on the effects of six pesticides on population > growth rate of a predatory mite. The response variable is the population > growth rate of the mite (ranges from negative to positive) and the > exploratory variable is a categorical variabl

Re: [R] do.call in "with" construction

2011-07-25 Thread Joshua Wiley
Dear Johannes, because myFun has an environment, and that the global environment. Look at: with(myEnv, { myFun2 <- function(symbols){do.call("print", lapply(symbols, FUN=as.name))} do.call("myFun2", list(symbols.env)) do.call("myFun", list(symbols.env)) }) Also, for example, compare

Re: [R] wrong name of input file and goto like function

2011-07-25 Thread MacQueen, Don
While you study the documentation as others suggested, may I suggest that you take a look at the file.exists() function. Here is an example of using file.exists(). tmpfile <- define() if (file.exists(tmpfile)) { ## read the file } else { ## tell the user to try again } On Jul 24, 2011, at

[R] compute the mean of two (or more) correlations

2011-07-25 Thread Martin Batholdy
Hi, is there a function in R that computes the mean of two (or more) correlations? (that is doing the z-transformation of the correlations, computing the mean of the z-values and then retransform it to a correlation). Or is there a fisher-table implemented? I have not found anything via help.

Re: [R] inside p value 'e'

2011-07-25 Thread Karl Ove Hufthammer
Dieter Menne wrote: >> I have like 5.075e-12 , 3.207e-05, 7.438e-07 and 9.393e-08 *** , i dont >> know what number they are > > http://en.wikipedia.org/wiki/Floating_point Or use ‘format’ in R. Example: format(5.075e-8, scientific=FALSE) -- Karl Ove Hufthammer __

[R] Binning numbers into integer-valued intervals (or: a version of cut or cut2 that makes sense)

2011-07-25 Thread Karl Ove Hufthammer
Dear list members, I’m looking for a way to divide numbers into simple (i.e., integer-valued) intervals, and thought the ‘cut’ function in ‘base’ or the ‘cut2’ function in ‘Hmisc’ would, er, cut it. However, they seem to give rather surprising results. Since I want the endpoints of the interva

Re: [R] load data

2011-07-25 Thread Dieter Menne
Crock wrote: > > i`m an absolute R-beginner and have a problem loading data from datasets. > For example, giving the command: >load ("Affairs"), i get the message, > that there is an error in readChar and there cannot be made any > connection. > I tried to find out on introduction papers but was

Re: [R] Rgraphviz installation problem

2011-07-25 Thread Martin Morgan
On 07/25/2011 06:05 AM, Ram H. Sharma wrote: Dear experts I installed Rgraphviz with the following command: source("http://www.bioconductor.org/biocLite.R";) biocLite("Rgraphviz") But had a problem when I tried to load it. This says that libcdt-4.dll is missing in your computer and th

[R] Getting more results from lda (MASS)

2011-07-25 Thread Stephen T
Hello, I am using linear discriminant analysis (lda) from the MASS library to classify data in two classes. 1. How do I get the full LDA model? The function lda reports coefficients but not the constant term? Currently I run the linear model function lm on the LDA scores and variables to find

[R] [R-pkgs] rdyncall 0.7.3

2011-07-25 Thread Daniel Adler
Initial Announcement: Package rdyncall released on CRAN. (Version 0.7.3) The package was presented at the Use!R 2009 with the title 'An improved Foreign Function Interface for R' and is now available on CRAN and considered stable for a large range of R platforms. The package provides a cross-pl

[R] How to find the likelihood of a null model in R

2011-07-25 Thread Partha Pratim PATTNAIK
Dear All, I am working on a dataset having the dependent variable as ordinal data(discrete data) and multiple independent variables. I need to find the likelihood for the NULL model.i.e the model with only the dependent variable and all other independent variables as zero. Kindly let me k

[R] smallest space analysis

2011-07-25 Thread marco
Hi everyone just a brief question: I'm trying to perform a smallest space analysis on nonmetric variables with R. It seems there are at lear a couple of packages to perform non metric MDS but still not found any info concerning SSA. Any clue? thanks in adv. Marco -- View this message in context:

[R] do.call in "with" construction

2011-07-25 Thread Johannes Egner
Dear all, I'd appreciate any help to rectify what must be a misconception of mine how environments work: ## myEnv <- new.env() myEnv$a.env <- 1 myEnv$symbols.env <- "a.env" a.global <- 2 symbols.global <- "a.global" myFun <- function(symbols){do.call("print", lapply(symb

[R] Rgraphviz installation problem

2011-07-25 Thread Ram H. Sharma
Dear experts I installed Rgraphviz with the following command: source("http://www.bioconductor.org/biocLite.R";) biocLite("Rgraphviz") But had a problem when I tried to load it. This says that libcdt-4.dll is missing in your computer and the following message in R window: > library("Rgr

[R] Cox's regression analysis with Left truncated data

2011-07-25 Thread natalie.vanzuydam
Hi, I have a fairly simple question. I would like to use the survival package to perform an analysis on data where an event can have occurred before individuals were recruited into a study. I'm not sure how to do this using the Surv() function. I would have a date of an event and then the enrol

Re: [R] Rcompression on MAC - where is it?

2011-07-25 Thread Eduardo M. A. M. Mendes
Hi there Yes, I am sure. I managed to compile rcompression and the package is now available. Many thanks Ed On Jul 24, 2011, at 7:17 PM, Henrik Bengtsson wrote: > Hi. > > On Tue, Jul 12, 2011 at 5:40 PM, Eduardo M. A. M. Mendes > wrote: >> Hello >> >> I need to run an r-file that works on

[R] load data

2011-07-25 Thread Crock
hi, i`m an absolute R-beginner and have a problem loading data from datasets. For example, giving the command: >load ("Affairs"), i get the message, that there is an error in readChar and there cannot be made any connection. I tried to find out on introduction papers but wasn`t able to find anythi

  1   2   >