Re: [R] C module causing rounding errors?

Insufficient data. Please read the posting guide, particularly the part about reproducible code. It might have something to do with "double" vs "float". "Joseph N. Paulson" wrote: >Hi all! > >I am currently writing a C-module for a for loop in which I permute columns >in a matrix (bootstrappin

Re: [R] Histogram Principal component analysis in R

Hello Dhiman , I have never tried doing it on such data, so I am not sure how this methodology applies or what covets it may hold. However, R wise, it sounds to me like simply reusing the code here: http://www.statmethods.net/advstats/factor.html Would do the trick. p.s: you wrote "It would be g

Re: [R] Brewer.pal error

Hi there, For me the package is working. Try to install.packages("RColorBrewer") It. And then library("RColorBrewer") Then something like: brewer.pal(7,"Greens") Write in the mailing list what errors you got (if any) at any of these commands. Tal Contact Details:-

Re: [R] Data Frame Manipulation using function

On Jul 8, 2010, at 10:33 PM, Erik Iverson wrote: I have a data frame: id url urlType 1 1 www.yahoo.com www.yahoo.com>1 2 2 www.google.com/?search=

Re: [R] Data Frame Manipulation using function

I have a data frame: id url urlType 1 1 www.yahoo.com 1 2 2 www.google.com/?search= 2 3 3

Re: [R] Data Frame Manipulation using function

On Jul 8, 2010, at 10:09 PM, harsh yadav wrote: Hi, Here is a somewhat detailed explanation of what I want to achieve: I have a data frame: id url urlType 1 1 www.yahoo.com1 2 2 www.google.com/?search= 2 3

Re: [R] Data Frame Manipulation using function

Hi, Here is a somewhat detailed explanation of what I want to achieve: I have a data frame: id url urlType 1 1 www.yahoo.com1 2 2 www.google.com/?search= 2 3 3 www.google.com

[R] C module causing rounding errors?

Hi all! I am currently writing a C-module for a for loop in which I permute columns in a matrix (bootstrapping) and I send a couple of variables into C initially. All of it is working, except the initial values I send to R are rounded/truncated (I believe rounded). I am using a 32 bit machine to

Re: [R] Data Frame Manipulation using function

It will be a lot easier to help you if you follow the posting guide and PLEASE do read the posting guide and provide commented, minimal, self-contained, reproducible code. You gave your function definition, which is good. Use ?dput to give us a small data.frame that can accurately show what y

[R] Data Frame Manipulation using function

Hi all, I have a data frame for which I want to limit the output by checking whether row values for specific column meets particular conditions. Here are the more specific details: I have a function that checks whether an input string exists in a defined list:- checkBaseLine <- function(s){ fo

Re: [R] Time value not sorting properly

Jared, I am not sure how you converted your 'time' variable from a factor to numeric, but you probably actually want to convert it to one of the 'time' classes. To learn more about them in R, see ?DateTimeClasses Another nice feature of these special time classes is that they can handle year, mon

Re: [R] xyplot help

Thank you for pointing out panel.superpose. I got it working with panel.xyplot as below: xyplot(y~x,groups=type,plotdata, auto.key=list(space='bottom',points=T) , panel = function(x,y,subscripts, groups) { panel.xyplot(x,y, col=groups[subscripts]) x<-sort(x) panel.lines(x,-1.324+0.11

[R] Multivariate Kalman Filter with inequality constraints

I'm trying to do a multivariate Kalman Filter with inequality constraints. Specifically, I'm trying to force all betas to be positive in a y = beta * x equation. Any ideas? Pierre Lapointe [[alternative HTML version deleted]] __ R-help@r-proj

[R] Time value not sorting properly

I have a dataframe of animal locations that I need to have in incremental order so that I can calculate the distance traveled between each time step. However, I have identified a few values that don't seem to sort properly. For instance, the last value in the table below should be the first value a

Re: [R] xyplot help

Read the "panel" part of ?xyplot carefully and note ?panel.superpose. Typically, when you have a groups argument your panel function must be of the form panel = function(x,y,...){ panel.superpose(x,y,...) etc. } However, I think it might all happen automagically if you do the cal

Re: [R] Producing residual plots by time for lme object

The plot also works if you replace "na.exclude" with "na.omit". On Wed, 7 Jul 2010, Deepayan Sarkar wrote: DS > On Fri, Jul 2, 2010 at 11:40 AM, Gregory Gentlemen DS > wrote: DS > > Fellow R-users, DS > > DS > > I have a longitudinal data set with missing values in it. I would like to produce

[R] xyplot help

Hi, I am learning xyplot. I have an example dataset attached. plotdata<-read.table("plotdata.txt",sep='\t',header=T,row.names=1) head(plotdata,n=4) y x type 1 -4.309601 -0.7448405A 2 -4.715421 0.7875994A 3 -2.310638 0.5455310A 4 -2.685803 10.4116868A xyplot(y

[R] New R-SIG for Discrete Choice Modelling

Hello all, I'd like to announce the availability of a mailing list for a newly-formed SIG (Special Interest Group) dedicated to using R for Discrete Choice Modelling. This list is intended for discussion of issues revolving around the design and analysis of Discrete Choice (aka Stated Choice, St

Re: [R] Column header strategy

On Thu, Jul 8, 2010 at 3:14 PM, Addi Wei wrote: > > Hopefully simple question: What is the best way to name, and treat factor > columns for data that has lots of columns? > > This is my column list: > id pID50 D.1 D.2 D.3 D.4 D.5 , etc. all the way to D.185 > It would be much better to use a fac

Re: [R] Column header strategy

On 2010-07-08 13:14, Addi Wei wrote: Hopefully simple question: What is the best way to name, and treat factor columns for data that has lots of columns? This is my column list: id pID50 D.1 D.2 D.3 D.4 D.5 , etc. all the way to D.185 I was under the impression from several R examples in pls

Re: [R] Column header strategy

On Thu, Jul 8, 2010 at 12:14 PM, Addi Wei wrote: > > Hopefully simple question:  What is the best way to name, and treat factor > columns for data that has lots of columns? > > This is my column list: > id pID50 D.1 D.2 D.3 D.4 D.5 , etc. all the way to D.185 > > I was under the impression from se

Re: [R] Loading .Rdata within an R function

On 08/07/2010 6:21 PM, Duncan Murdoch wrote: On 08/07/2010 3:21 PM, Giles Crane wrote: Colleagues: I am having trouble loading data from within .Rdata file within the environment of a function. That is, the following always loads to the global environment: f1 <- function(){

Re: [R] Loading .Rdata within an R function

On 08/07/2010 3:21 PM, Giles Crane wrote: Colleagues: I am having trouble loading data from within .Rdata file within the environment of a function. That is, the following always loads to the global environment: f1 <- function(){ load("mydata.Rdata") # compute w

[R] Histogram Principal component analysis in R

Hi, I am trying to do a Principal component analysis on histogram data. Basically, I have a group of subjects and for each of them, I have a column of bin-counts (vis-a-vis intervals) and a corresponding column of frequencies (or normalized frequencies). The bin counts are the same for all the sub

Re: [R] Loading .Rdata within an R function

On Jul 8, 2010, at 3:21 PM, Giles Crane wrote: Colleagues: I am having trouble loading data from within .Rdata file within the environment of a function. That is, the following always loads to the global environment: f1 <- function(){ load("mydata.Rdata") # compute

Re: [R] Column header strategy

On Jul 8, 2010, at 3:14 PM, Addi Wei wrote: Hopefully simple question: What is the best way to name, and treat factor columns for data that has lots of columns? This is my column list: id pID50 D.1 D.2 D.3 D.4 D.5 , etc. all the way to D.185 I was under the impression from several R exam

Re: [R] What is the name of function for calculating mod?

?Arithmetic , specifically %% Soyeon Kim wrote: Dear All, I try to find the name of function to calculate mod(or Modular). Do you anyone know the name of function? Thank you so much, Soyeon Kim __ R-help@r-project.org mailing list https://stat.ethz.

[R] What is the name of function for calculating mod?

Dear All, I try to find the name of function to calculate mod(or Modular). Do you anyone know the name of function? Thank you so much, Soyeon Kim __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posti

[R] Brewer.pal error

Hi all, I have 35 years, I have been working since I exit from University. I have returned to make a phd...I´m fortunatelly I´m working so unfortunatelly I have few time to study (I have also a baby with 16 months). I´m trying to use R but sometimes it gets very difficult. I´m trying to put this

[R] ttrTests Error

spData <- as.vector(getYahooData("SPY", start="19900101",end="20081231")[,"Close"]) > cr <- cReturns(spData, ttr = "sma", params=c(20)) Error in ind[t - k] <- pos[t - k + 1] - pos[t - k] : replacement has length zero I am getting the above error when running the cReturns function. Any ideas on t

[R] Loading .Rdata within an R function

Colleagues: I am having trouble loading data from within .Rdata file within the environment of a function. That is, the following always loads to the global environment: f1 <- function(){ load("mydata.Rdata") # compute with objects from mydata.Rdata.

Re: [R] R2 function from PLS to use a model on test data

I figured out the answer. Since I have both the predicted and actual values, I simply need to call the lm function, and summary to see R^2 which matches what MOE provided. myline.fit <-lm( y~x) summary(myline.fit) -- View this message in context: http://r.789695.n4.nabble.com/R2-function-fro

[R] Column header strategy

Hopefully simple question: What is the best way to name, and treat factor columns for data that has lots of columns? This is my column list: id pID50 D.1 D.2 D.3 D.4 D.5 , etc. all the way to D.185 I was under the impression from several R examples in pls that if you name your columns like ab

Re: [R] Kite diagrams

I asked the same question on R-sig-eco, and Ben Bolker provided this solution, which as I assume this should show up in a search I copy here. However, if someone can come up with a single function, that would be good. --

Re: [R] Function for gruping similar variables?

Does anyone have an idea on this? On 6 Lip, 11:32, Timo W wrote: > Hi, > > I have a matrix of results of multiple 2x2 chi^2 tests, non- > significant tests are marked as TRUE. Is there a function for grouping > those variables in a similar way LSD.test from agricolae library does? > I reviewed LS

Re: [R] Function for gruping similar variables?

Does anyone have an idea on this? On 6 Lip, 11:32, Timo W wrote: > Hi, > > I have a matrix of results of multiple 2x2 chi^2 tests, non- > significant tests are marked as TRUE. Is there a function for grouping > those variables in a similar way LSD.test from agricolae library does? > I reviewed LS

Re: [R] Function for gruping similar variables?

Does anyone have an idea on this? On 6 Lip, 11:32, Timo W wrote: > Hi, > > I have a matrix of results of multiple 2x2 chi^2 tests, non- > significant tests are marked as TRUE. Is there a function for grouping > those variables in a similar way LSD.test from agricolae library does? > I reviewed LS

Re: [R] plotting lines when data has missing/NA values

On 08-Jul-10 19:52:36, Justin Fincher wrote: > Howdy, >I have done many searches and can't seem to find a way around this. > I am reading in a .csv file where each row is a dataset and each > column represents a position. The values are sparse (there are 2003 > positions but usually only 100-

Re: [R] combining P values from different tests

On 08/07/2010 2:45 PM, Jonathan Baron wrote: On 07/08/10 19:25, Knut Krueger wrote: > Hi to all, > is there any function available to combine those p values; > > > Kind Regards Knut Here's a couple of functions I've used for one thing that you might mean. The idea is that you have two p-valu

Re: [R] Error in which()

On Jul 8, 2010, at 3:23 PM, Muhammad Rahiz wrote: Hi all, I'm trying to filter data into respective numbers. For example, if the data ranges from 0 to <0.1, group the data. And so on for the rest of the data. There are inconsistencies in the output. For example, b1[[3]] lumps all the 0.2

[R] plotting lines when data has missing/NA values

Howdy, I have done many searches and can't seem to find a way around this. I am reading in a .csv file where each row is a dataset and each column represents a position. The values are sparse (there are 2003 positions but usually only 100-200 values) and the idea is to be able to plot each dat

Re: [R] Error in which()

Congratulations, you've run afoul of: 7.31 Why doesn't R think these numbers are equal? > x0[7] [1] 0.6 > > x0[7] == .6 [1] FALSE > Sarah On Thu, Jul 8, 2010 at 3:23 PM, Muhammad Rahiz wrote: > Hi all, > > I'm trying to filter data into respective numbers. For example, if the data > ranges fr

[R] Error in which()

Hi all, I'm trying to filter data into respective numbers. For example, if the data ranges from 0 to <0.1, group the data. And so on for the rest of the data. There are inconsistencies in the output. For example, b1[[3]] lumps all the 0.2s and 0.3s together while 0.6s are not in the output.

Re: [R] combining P values from different tests

On Thu, 8 Jul 2010, David Winsemius wrote: On Jul 8, 2010, at 2:17 PM, Stephan Kolassa wrote: Knut Krueger schrieb: Hi to all, is there any function available to combine those p values; Hi Knut, ?"+" That is certainly in the running for shortest possible fortune. Added to the devel-

Re: [R] combining P values from different tests

On 07/08/10 19:25, Knut Krueger wrote: > Hi to all, > is there any function available to combine those p values; > > > Kind Regards Knut Here's a couple of functions I've used for one thing that you might mean. The idea is that you have two p-values for different tests of the same hypothesis, a

Re: [R] combining P values from different tests

On Jul 8, 2010, at 2:17 PM, Stephan Kolassa wrote: Knut Krueger schrieb: Hi to all, is there any function available to combine those p values; Hi Knut, ?"+" That is certainly in the running for shortest possible fortune. > PLEASE do read the posting guide > http://www.R-project.org/po

[R] Query about using timestamps returned by SQL as 'factor' for split

I have a simple query as follows: "SELECT m_id,sale_date,YEAR(sale_date),WEEK(sale_date),return_type,DATEDIFF(return_date,sale_date) AS elapsed_time FROM risk_input" I can get, and view, all the data that that query returns. The question is, sale_date is a timestamp, and I need to call split to

[R] Psudeo R^2 (or other effect size) in spatial gls regressions

Dear all, I have been using the function gls in the package nlme in R to fit some spatial regressions (as described in Dormann et al.). However, I have been struggling trying to find a way to calculate a measure of effect size from these models, so I wanted to know if any of you had an idea on

Re: [R] combining P values from different tests

On 08/07/2010 1:25 PM, Knut Krueger wrote: Hi to all, is there any function available to combine those p values; No, not in general, because the proper way to do it depends on the joint distribution of the inputs, and what you want to use the combined value to test. Duncan Murdoch

Re: [R] combining P values from different tests

Hi Knut, ?"+" > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. Best, Stephan Knut Krueger schrieb: Hi to all, is there any function available to combine those p values; Kind Regards Knut _

Re: [R] Using nlm or optim

Hi All, Thanks for all of your help. I hadn't posted error messages, because I had tried many different formulations of those functions and each one gave a different error message. However, using your suggestions, I managed to make it work. You were right in that I only needed to call the constan

[R] combining P values from different tests

Hi to all, is there any function available to combine those p values; Kind Regards Knut __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and

Re: [R] Custom multi-line header for pair plots

Have a look at the text.panel parameter in ?pairs Cheers, 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 (

Re: [R] How can i draw a graph with high and low data points

Hi Nathaniel, Here are a few links for (short) articles that can help you get comfortable in doing such tasks: http://www.statmethods.net/graphs/line.html http://www.statmethods.net/graphs/scatterplot.html Here is a tiny example: x <- 1:100 y <- rnorm(100) plot(y~x) lines(loess(y~x)$fitted ~ x

Re: [R] download gplots and caTools

Hi Tra, Try another mirror, and worst case - just download it from the CRAN website Cheers, Tal Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talgalili.com (Hebrew) | www.biostatistics.co.il

Re: [R] glmnet - choosing the number of features

Hi Brian, On Wed, Jul 7, 2010 at 10:54 PM, Brian Tsai wrote: > Hi, > > I am trying to use the glmnet package to do some simple feature selection. >  However,  I would ideally like to be able to specify the number of features > to return (the glmnet package, as far as I can tell, only allows > spe

Re: [R] Vector Assignment -- Why won't it let me do it?

That was it! You honestly made my week, thank you for that. -- View this message in context: http://r.789695.n4.nabble.com/Vector-Assignment-Why-won-t-it-let-me-do-it-tp2281698p2282562.html Sent from the R help mailing list archive at Nabble.com. __

[R] How can i draw a graph with high and low data points

Hi, This is my first time so bear with me if I do anything silly. I'll learn! I'm a relatively new user of R as well so will appreciate any help. I have a data set in the form Trial.Group MeanHeart Rate Upper confidence level Lower confidence level 333subj

Re: [R] Calling Gnuplot from R

See this page for an alternative approach. It shows how you can use R with gnuplot. http://orgmode.org/worg/org-contrib/babel/examples/data-collection-analysis.php Octave is also supported by org-babel. Matt Shotwell wrote: I recently wrote a small R function to draw simple ASCII scatterpl

[R] download gplots and caTools

Dear list, I am using R.2.11.1 version on a PC. I downloaded successfully gplots_2.8.0 but I could not find where to download caTools > library("gplots") Loading required package: caTools Error: package 'caTools' could not be loaded In addition: Warning message: In library(pkg, character.only =

Re: [R] Calling Gnuplot from R

I recently wrote a small R function to draw simple ASCII scatterplots. http://biostatmatt.com/archives/491 Bill Harris commented that the plots reminded him of the "dumb" terminal of Gnuplot. I think it would be really neat to have an R graphics driver to Gnuplot in order to generate more complete

Re: [R] How do I test against a simple null that tworegressions coefficients are equal?

Bert, Thanks for replying. Actually, Chuck and Mark's suggestions already helped, especially regarding to the problem you call type III error. Best, Jia On Thu, Jul 8, 2010 at 11:06 AM, Bert Gunter wrote: > Inline below. > > Bert Gunter > Genentech Nonclinical Statistics > >> -Original Mess

Re: [R] How do I test against a simple null that two regressions coefficients are equal?

OK. Thanks again. I will read the references more. Best, Jia On Thu, Jul 8, 2010 at 10:51 AM, wrote: > hi: no. it's not the same. if you read the paper that I referenced last > night, that explains how to do the following  test : > > Ho: R2 = R1 > > H1: R2 != R1 > > that's a different test fro

Re: [R] Calling Gnuplot from R

If you use Emacs, you can use org-mode with org-babel to facilitate this... I'll refrain from asking why :). See: http://orgmode.org/worg/org-contrib/babel/index.php Christopher Desjardins wrote: Hi, I am wondering if there is a way to call Gnuplot from R and/or if anyone can recommend a packa

[R] Calling Gnuplot from R

Hi, I am wondering if there is a way to call Gnuplot from R and/or if anyone can recommend a package on CRAN capable of doing this? Thanks, Chris PS - Please cc me on the response. [[alternative HTML version deleted]] __ R-help@r-project.org mai

Re: [R] How do I test against a simple null that tworegressions coefficients are equal?

Inline below. Bert Gunter Genentech Nonclinical Statistics > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] > On Behalf Of markle...@verizon.net > Sent: Thursday, July 08, 2010 7:52 AM > To: chen_1...@fisher.osu.edu > Cc: r-help@r-project.org

Re: [R] F# vs. R

On Thu, Jul 8, 2010 at 4:44 AM, Patrick Burns wrote: > I'd like to hear answers to this as well. > A language doesn't have to be a complete > replacement to be useful. Totally agree. > F# seems to have some nice features. My intention isn't to have this explode into a language war, so apologies

Re: [R] F# vs. R

It's been a long time since I used Fortran, and I have only dabbled in F#, but I do not think translating Fortran (or R) to F# will be easy. F# is basicly a functional language (like ML) and a very differant mind set than Fortran (or R). -EdK Ed Keith e_...@yahoo.com Blog: edkeith.blogspot.

Re: [R] How do I test against a simple null that two regressions coefficients are equal?

hi: no. it's not the same. if you read the paper that I referenced last night, that explains how to do the following test : Ho: R2 = R1 H1: R2 != R1 that's a different test from what you did but i think it's what you want. On Jul 8, 2010, chen jia wrote: Thanks, Chuck.

Re: [R] How do I test against a simple null that two regressions coefficients are equal?

Thanks, Chuck. I am reading the references, which are helpful. Just to understand what I have done wrong here, I proposed an alternative testing strategy: I run regressions y = a3 + b1 * x + b2 * z + e3 and test alternative hypothesis b1 != b2 against the null hypothesis b1 = b2 in this equation

[R] Custom multi-line header for pair plots

Hi All, Is there a possible way to customize the header text appearing on the diagonal of a pair plot? My header is: Biomass PoB QoB ... etc. I want to put additional information in the diagonal boxes where header text appears such as mean, standard deviation of data. If I add these information

[R] ANOVA-Formula

Hi, I have taken one microarray experiment and trying to implement same statistical measures what they have done.I have taken datasets from GEO platform with accession number GSE1557. In the experiment,about half of double transgenic rats (dTGR) over-expressing the human renin and angiotensinogen

Re: [R] F# vs. R

True, porting old C and Fortran code to C# or F# would be a pain and probably riddled with errors but it is not too soon to start looking to see if there is a better way. There have been numerous ports of LAPACK, BLAS, etc. to C#. Maybe they could be leveraged. Maybe just allowing packages to b

Re: [R] package installation for Windows 7

On Thu, Jul 8, 2010 at 3:39 PM, Duncan Murdoch wrote: > On 08/07/2010 9:26 AM, David Bickel wrote: >> >> Yes, the User into which I logged in before launching RGui is an >> Administrator. Correct, the problem is not limited to Bioconductor packages. >> > > On Windows 7, it's not enough to have the

Re: [R] package installation for Windows 7

Thanks, Joris. Your suggestion solved the problem without requiring me to run R as an administrator. Best regards, David On 08/07/2010 8:45 AM, Joris Meys wrote: Hi, I am running Windows 7 and R 2.11.1, and everything is installing just fine for me. Did you install R in the "Program Files" f

Re: [R] Appropriateness of survdiff {survival} for non-censored data

The query: "Thus I question the appropriateness of using survdiff in my analysis; I have exact data yet I would be testing on the Kaplan-Meir estimate of these data in survdiff. Thanks for any help." My thoughts: There are two aspects of survival analysis you need to think of. The first, as y

Re: [R] package installation for Windows 7

Yes, the User into which I logged in before launching RGui is an Administrator. Correct, the problem is not limited to Bioconductor packages. David On 08/07/2010 8:55 AM, Duncan Murdoch wrote: On 08/07/2010 7:15 AM, Dave Bickel wrote: Neither biocLite nor the GUI menus can install packages o

Re: [R] Using nlm or optim

Without data I can't check, but try : mle(nll,start=list(c=0.01,z=2.1,s=200),fixed=list(V=Var,M=Mean)) With a random dataset I get : > Mean <- rnorm(136) > Var <- 1 + rnorm(136)^2 > mle(nll,start=list(c=0.01,z=2.1,s=200),fixed=list(V=Var,M=Mean)) Error in optim(start, f, method = method, hessian

Re: [R] package installation for Windows 7

Yes, the User into which I logged in before launching RGui is an Administrator. Correct, the problem is not limited to Bioconductor packages. David On 08/07/2010 8:55 AM, Duncan Murdoch wrote: On 08/07/2010 7:15 AM, Dave Bickel wrote: Neither biocLite nor the GUI menus can install packages on

Re: [R] random sample from arrays

On Jul 8, 2010, at 2:04 AM, Assa Yeroslaviz wrote: Hello R users, I'm trying to extract random samples from a big array I have. I have a data frame of over 40k lines and would like to produce around 50 random sample of around 200 lines each from this array. this is the matrix ID

Re: [R] Using nlm or optim

Hi, When you report errors please send a reproducible example that will allow us to better help you. At a minimum, you should cut and paste the actual call and the error message. I think the problem is that you should only have parameters as arguments, i.e. only c and z should be in the argum

[R] 0 exit status on packages update try

I am not entirely sure what is going on. I am sure that I am missing something. maptools won't detach... Thanks for all of your help. R 2.11.1 OS Ubuntu 10.04 R --vanilla install.packages("sp", dep=TRUE) Loading required package: reshape Loading required package: plyr Loading required packag

Re: [R] xyplot: key inside the plot region / lme: confidence bands for predicted

Kingsford Jones wrote: On Wed, Jul 7, 2010 at 8:59 AM, Michael Friendly wrote: No one replied to my second question: how to get standard errors or confidence intervals for the estimated fixed effects from lme().AFAICS, intervals() only gives CIs for coefficients. Assuming you mean

Re: [R] package installation for Windows 7

On 08/07/2010 9:26 AM, David Bickel wrote: Yes, the User into which I logged in before launching RGui is an Administrator. Correct, the problem is not limited to Bioconductor packages. On Windows 7, it's not enough to have the user be an administrator, you need to run programs specifically

Re: [R] strsplit("dia ma", "\\b") splits characterwise

On Thu, Jul 8, 2010 at 4:15 AM, Suharto Anggono Suharto Anggono wrote: > \b is word boundary. > But, unexpectedly, strsplit("dia ma", "\\b") splits character by character. > >> strsplit("dia ma", "\\b") > [[1]] > [1] "d" "i" "a" " " "m" "a" > >> strsplit("dia ma", "\\b", perl=TRUE) > [[1]] > [1] "

[R] [OOPS] Re: hi... problems about adjacency matrix

See at end ... On 08-Jul-10 12:41:03, Ted Harding wrote: > On 08-Jul-10 10:33:56, Gina Liao wrote: >> Dear all, >> Hi, I have the problems about converting the matrix to adjacency >> matrix.Here's my example, >> ab cd >> e

Re: [R] random sample from arrays

Don't know what exactly you're trying to do, but you make a matrix with 11 columns and 50 rows, then treat it as a vector. On top of that, you try to fill 50 rows/columns with 50 values. Off course that doesn't work. Did you check the warning messages when running the code? Either do : for(i in

Re: [R] Adding two files into one and vlookup

Hi: > merge(file2, file1, by = 'Date') DatePrice.x Price.y 1 03/07/2010 5.312006 56.92826 2 04/07/2010 673.070592 39.27409 3 05/07/2010 442.467939 42.59834 4 06/07/2010 851.915899 70.68512 5 07/07/2010 581.859242 10.92505 HTH, Dennis On Wed, Jul 7, 2010 at 8:07 AM, Raghu wrote:

Re: [R] package installation for Windows 7

As far as my experience goes, there is no need whatsoever to run R as an administrator for installation of any package, including BioConductor, provided you stay away from the "Program Files" folder. Cheers Joris On Thu, Jul 8, 2010 at 2:55 PM, Duncan Murdoch wrote: > On 08/07/2010 7:15 AM, Dave

Re: [R] strsplit("dia ma", "\\b") splits characterwise

l guess this is expected behaviour, although counterintuitive. \b represents an empty string indicating a word boundary, but is coerced to character and thus simply the empty string. This means the output you get is the same as > strsplit("dia ma", "",perl=T) [[1]] [1] "d" "i" "a" " " "m" "a" I'd

Re: [R] forcing a zero level in contr.sum

Peter, I am asked to show seasonal effects vs a general averaged level without any seasonal effects. That's why contr.sum and not treatment. The estimates show that some coefs are not stat different from 0, which is also the enforced total average. So, my understanding of model selection leads m

Re: [R] hi... problems about adjacency matrix

Hi Gina, On Thu, Jul 8, 2010 at 6:33 AM, Gina Liao wrote: > > Dear all, > Hi, I have the problems about converting the matrix to adjacency > matrix.Here's my example, >          a            b               c                d                 e     >           fa 1.000 0.4048823682 0.1228531

Re: [R] can't open file

On Thu, Jul 8, 2010 at 8:47 AM, Sebastian Kruk wrote: >>> I tried: >>> >>> L <- readLines(con <- file("log2.log", encoding="UCS-2LE") >>> z <- read.zoo(textConnection(L), index = 1:2,FUN = function(x) >>> paste(x[,1], x[,2])) >>> >>> Error: >>> Error en x[, 1] : número incorreto de dimensiones >>>

Re: [R] package installation for Windows 7

On 08/07/2010 7:15 AM, Dave Bickel wrote: Neither biocLite nor the GUI menus can install packages on my system. This is probably a permissions problem. Are you running R as an administrator when you try the install? If not, you won't be able to install to the default location, but you sh

Re: [R] how to determine order of arima bt acf and pacf

Did you read these already? http://www.statoek.wiso.uni-goettingen.de/veranstaltungen/zeitreihen/sommer03/ts_r_intro.pdf http://www.barigozzi.eu/ARIMA.pdf Cheers Joris On Thu, Jul 8, 2010 at 11:45 AM, vijaysheegi wrote: > > Hi R community, > I am new to R.Have  some doubts on ACF and pacf.Appy

Re: [R] can't open file

>> I tried: >> >> L <- readLines(con <- file("log2.log", encoding="UCS-2LE") >> z <- read.zoo(textConnection(L), index = 1:2,FUN = function(x) >> paste(x[,1], x[,2])) >> >> Error: >> Error en x[, 1] : número incorreto de dimensiones >> >Please provide a reproducible example like this: >library(zo

Re: [R] package installation for Windows 7

Hi, I am running Windows 7 and R 2.11.1, and everything is installing just fine for me. Did you install R in the "Program Files" folder? If so, uninstall and try to re-install R in another folder (e.g. c:\R\R2.11.1\ like on my computer). I noticed in the past that the access control of Windows tre

Re: [R] hi... problems about adjacency matrix

On 08-Jul-10 10:33:56, Gina Liao wrote: > Dear all, > Hi, I have the problems about converting the matrix to adjacency > matrix.Here's my example, > ab cd > e fa 1.000 0.4048823682 0.1228531 0.49046991 > 0.494515886

Re: [R] how to process this in R

I don't have R on the netbook I'm writing this from, but I *think* the following will work: > 1. how can i only read in the files with the string patterns ggg or fff as > part of the file names? >   for instance, I only need the file names with the ggg or fff in it >      x_ggg_y_1.txt >  

Re: [R] x-axis labels for barplot of zoo object

On Thu, Jul 8, 2010 at 12:42 PM, Gabor Grothendieck wrote: > On Thu, Jul 8, 2010 at 4:17 AM, Rainer M Krug wrote: > > Hi > > > > Let's assume, I have a dataset of 1000 datapoints, which represent daily > > recordings of a measurement. > > They are stored as a ts object (see example below) > > >

  1   2   >