Re: [R] SAS and R on multiple operating systems

Douglas Bates wrote: The first VAX at our university, a VAX-11/780, was allegedly a 1 MIPS machine but it was pretty difficult to get it to do a million of any instruction in one second. I eventually succeeded because it had a special instruction for decrement and branch on zero so if you

Re: [R] predict.lm

Hi Felipe, Make sure predictors data.frame has the EXACT same column names as the data.frame you used for creating ModeloLineal. Cheers, Tal Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me:

Re: [R] NMDS Ordination Graphics Problem

hi, maybe you should see: http://cc.oulu.fi/~jarioksa/softhelp/FAQ-vegan.html#How-to-use-different-plotting-symbols-in-ordination-graphics_003f greetings, kay -- View this message in context: http://n4.nabble.com/NMDS-Ordination-Graphics-Problem-tp1751845p1752557.html Sent from the R help

Re: [R] bootstrap confidence intervals, non iid

Kay Cichini wrote: i need to calculate ci's for each of 4 groups within a dataset, to be able to infere about differences in the variable similarity. the problem is that data within groups is dependent, as assigned by the blocking-factor site. my guess was to use a block bootstrap but

Re: [R] bootstrap confidence intervals, non iid

not yet, but it would either not solve the non iid problem, or would it? -- View this message in context: http://n4.nabble.com/bootstrap-confidence-intervals-non-iid-tp1751619p1752583.html Sent from the R help mailing list archive at Nabble.com. __

[R] GridR

Dear All I am doing one distributed data mining program. So I selected GridR package for the distributed programming. *grid.init(service=local,debug=FALSE, localTmpDir=GridRTmp/) grid.apply(x,UCS, wait=TRUE )* UCS is a function. When I execute this statement, I have some errors like : *cannot

[R] about ARMA(p,q) SCAN method: SAS vs. R

Hi all, I am modifying a program I wrote before to perform smallest canonical (SCAN) correlation method for identification of ARMA(p,q) orders in Time Series, but when I compared the output with SAS, there are some differences. My SCAN R code can be downloaded in the following URL:

Re: [R] Palette color order in bwplot (lattice violin plot) vs. boxplot

On 04/04/2010 06:21, Deepayan Sarkar wrote: On Thu, Apr 1, 2010 at 4:10 AM, Luigi Ponti lpo...@inbox.com wrote: Hello, I am trying to give different colors to boxes in a violin plot obtained via bwplot from lattice package using a color palette from RColorBrewer: [...] However, when I

Re: [R] GridR

*Always* report sessionInfo() and traceback() when reporting problems. /Henrik On Tue, Apr 6, 2010 at 11:11 AM, wesley mathew wesleycmat...@gmail.com wrote:  Dear All I am doing one distributed data mining program. So I selected GridR package for the distributed programming.

Re: [R] Error with read.csv.sql on processing large file

1. Are you able to determine if its the filtering step that is the problem? Also is it the size of the file or is it a particular record that is causing the problem? If you read it in in two pieces do both pieces work or just one piece? It might be a problem in the filtering step (in which

Re: [R] logistic regression in an incomplete dataset

Hi Bert, Thanks for your reply. I AM making an assumption of MAR data, because informative missingness (I assume you mean NMAR) is too hard to deal with I have quite a few covariates (so the observed is likely to predict the missing and mitigate against informative missingness) the

Re: [R] logistic regression in an incomplete dataset

Dear Thomas, Thanks for your reply. Yes you are quite right (your example) complete case does not require MCAR, however as well as being a bit less robust than ML it is throwing away data. Missing Data in Clinical Studies, Geert Molenberghs, Michael Kenward, have a nice section in chapter 3

[R] is.na- doesn't seem to work with labelled variables?

Dear All, I seem entirely unable to solve what should be a very simple problem. I have imported a SPSS dataset into R using spss.get from Frank Harrell's Hmisc library. The variables are imported classed as 'labelled': missing values are coded as either the SPSS missing value code, 8 or 88.

[R] respecting original matrix dimensions

Hi All, I'm hoping someone else can help me out with this. I am doing some matrix algebra using sub-parts of matrices, and sometimes I need only a single row/column of the original matrix. However, whenever I pull out only a single row/column, R returns a row vector, but often this will break my

Re: [R] is.na- doesn't seem to work with labelled variables?

a - c(1, 8, 9, 10, 8) a[which(a==8)] - NA Nikhil On Apr 6, 2010, at 7:31 AM, David Foreman wrote: Dear All, I seem entirely unable to solve what should be a very simple problem. I have imported a SPSS dataset into R using spss.get from Frank Harrell's Hmisc library. The variables are

Re: [R] is.na- doesn't seem to work with labelled variables?

On Apr 6, 2010, at 6:31 AM, David Foreman wrote: Dear All, I seem entirely unable to solve what should be a very simple problem. I have imported a SPSS dataset into R using spss.get from Frank Harrell's Hmisc library. The variables are imported classed as 'labelled': missing values are

Re: [R] respecting original matrix dimensions

On Apr 6, 2010, at 6:39 AM, Robert M. Flight wrote: Hi All, I'm hoping someone else can help me out with this. I am doing some matrix algebra using sub-parts of matrices, and sometimes I need only a single row/column of the original matrix. However, whenever I pull out only a single

[R] checking bivariate normality

x - iris$Sepal.Length[1:50]/iris$Sepal.Width[1:50] y - iris$Petal.Length[1:50]/iris$Petal.Width[1:50] I want to check whether (x,y) follows a bivariate normal distribution or not, using density plot or scatter plot. Is it possible to plot a bivariate density in R. I cant find any. Arindam

Re: [R] respecting original matrix dimensions

yes, there is -- have a look at the online help file ?[, especially argument 'drop'. Try also this: crud - matrix(c(1,2,3,4,5,6,7,8,9,10,11,12),nrow=3,ncol=4) crud[, 1, drop = FALSE] I hope it helps. Best, Dimitris On 4/6/2010 1:39 PM, Robert M. Flight wrote: Hi All, I'm hoping someone

Re: [R] R 2.6 Support Question

Many others have made good comments regarding your questions, my suggestion would be this: You may or may not know that R is freely available and that you can install more than one version of R on the same PC, so I would suggest installing both R 2.6 and the latest version of R (2.10.x) and try

[R] Extracting formulae from expression() / deriv()

I am attempting to extract the derivative/ gradient from this expression df1p - deriv(f1, P) df1p expression({ .value - s - c - a * P .grad - array(0, c(length(.value), 1L), list(NULL, c(P))) .grad[, P] - -a attr(.value, gradient) - .grad .value }) So in this case I want

Re: [R] changing column names in a dataframe

jda wrote: I have imported data from an Excel spreadsheet. Columns in that spreadsheet are named name, x, and y, and several sets of those columns appear in the worksheet. For example: name x y name x y test1 1 3 test2 4 4 test1 2 2 test2 5 5 test1 3 1 test2 6 6 When I

[R] svm of e1071 package

Hello List, I am having a great trouble using svm function in e1071 package. I have 4gb of data that i want to use to train svm. I am using Amazon cloud, my Amazon Machine Image(AMI) has 34.2 GB of memory. my R process was killed several times when i tried to use 4GB of data for svm. Now I am

[R] median of two groups

Dear all, What is the right test to test whether the median of two groups are statistically significant? Is it the wilcox.test, mood.test or the ks.test? In the text book I have got there is explanation for the Wilcoxon (Mann Whitney) test which tests ob the two variable are from the same

Re: [R] Extracting formulae from expression() / deriv()

Try this: attributes(eval(df1p))$gradient[,'P'] On Tue, Apr 6, 2010 at 9:09 AM, nickymcp nickymcpher...@gmail.com wrote: I am attempting to extract the derivative/ gradient from this expression df1p - deriv(f1, P) df1p expression({    .value - s - c - a * P    .grad - array(0,

Re: [R] Extracting formulae from expression() / deriv()

attributes(eval(df1p))$gradient[,'P'] Error in eval(expr, envir, enclos) : object 's' not found -- View this message in context: http://n4.nabble.com/Extracting-formulae-from-expression-deriv-tp1752738p1752804.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] Extracting formulae from expression() / deriv()

Post your f1 function to us On Tue, Apr 6, 2010 at 10:23 AM, nickymcp nickymcpher...@gmail.com wrote: attributes(eval(df1p))$gradient[,'P'] Error in eval(expr, envir, enclos) : object 's' not found -- View this message in context:

[R] lattice package: line end style

First, apologies for no example data but I don't think it's needed in this case, Q: Can (and if so how ) the line end style be changed for 'cloud' plots? I've tried par(lend=2), trellis.par.set(add.line = list(lend=2)) and much googling but to no avail Thanks in advance Dan P.S. the reason

Re: [R] Extracting formulae from expression() / deriv()

f1 ~s - c - a * P -- View this message in context: http://n4.nabble.com/Extracting-formulae-from-expression-deriv-tp1752738p1752816.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

Re: [R] Extracting formulae from expression() / deriv()

Try fixing the constants s - 2 c - 3 a - 2 P - 1 attributes(eval(df1p))$gradient[,'P'] On Tue, Apr 6, 2010 at 10:30 AM, nickymcp nickymcpher...@gmail.com wrote: f1 ~s - c - a * P -- View this message in context:

Re: [R] Data manipulation problem

OK... next question.. Which is still a data manipulation problem so I believe the heading is still OK. ##So now I read my population data from excel. pop-read.csv(pop.csv) typeof(pop) ## yields a list where I have age-specific population rows and a yearly column population, where the years are

Re: [R] Extracting formulae from expression() / deriv()

Thanks for your help - this works and I get s - 2 c - 3 a - 2 P - 1 attributes(eval(df1p))$gradient[,'P'] P -2 What I really wanted was an expression that I can use later - do you think there's any way to do this? -- View this message in context:

Re: [R] checking bivariate normality

On Apr 6, 2010, at 7:48 AM, arindam fadikar wrote: x - iris$Sepal.Length[1:50]/iris$Sepal.Width[1:50] y - iris$Petal.Length[1:50]/iris$Petal.Width[1:50] I want to check whether (x,y) follows a bivariate normal distribution or not, using density plot or scatter plot. Is it possible to plot

Re: [R] Extracting formulae from expression() / deriv()

df1p is a expression, and you can use later. On Tue, Apr 6, 2010 at 11:00 AM, nickymcp nickymcpher...@gmail.com wrote: Thanks for your help - this works and I get s - 2 c - 3 a - 2 P - 1 attributes(eval(df1p))$gradient[,'P']  P -2 What I really wanted was an expression that I can use

Re: [R] SAS and R on multiple operating systems

Roger DeAngelis(xlr82sas) wrote: Hi, About the forest plot. Some Phrarma companies demand the report and graphics follow very restrictive layouts. Thank goodness that the FDA does not require that. Frank SAS allows uses to use one template for graphs and tables. Margins have

[R] [R-pkgs] betareg 2.2-2: Beta regression

Dear useRs, version 2.2-2 of the betareg package has just been released on CRAN http://CRAN.R-project.org/package=betareg accompanied by an article in the Journal of Statistical Software http://www.jstatsoft.org/v34/i02/ The package provides beta regression for data in the unit interval

[R] [R-pkgs] Formula 1.0-0: Model formulas with multiple parts and responses

Dear useRs, version 1.0-0 of the Formula package has just been released on CRAN http://CRAN.R-project.org/package=Formula accompanied by an article in the Journal of Statistical Software http://www.jstatsoft.org/v34/i01/ It provides simple infrastructre for processing formulas like y

Re: [R] SAS and R on multiple operating systems

 Also I have seen 5,000 page listings in SAS. Is this a pro or a con? Hadley -- Assistant Professor / Dobelman Family Junior Chair Department of Statistics / Rice University http://had.co.nz/ __ R-help@r-project.org mailing list

[R] approx function wierd result

Dear R-list members, I am calculating the linear extrapolation for a data set, using the function found in Hmisc. x=c(0.0265,-0.0003,0.0142,0.0263,0.0634,0.1145,0.2504) y=c(58,107,152,239,362,512,724) x1=0.0393216 approxExtrap(x,y,x1, method=linear) approx(x,y,x1) #to see what is happening:

Re: [R] c(), or cat(), or paste(), all cause unwanted reordering

On Thu, 2010-03-25 at 11:54 -0800, Jeff Brown wrote: Wow, you guys are awesome. Thanks! Thanks for the cat() question Jeff and to all guRus out there for the replies. This is something I was looking for the last hour. Nikos __ R-help@r-project.org

Re: [R] Extracting formulae from expression() / deriv()

On Tue, Apr 6, 2010 at 9:09 AM, nickymcp nickymcpherson at gmail.com https://stat.ethz.ch/mailman/listinfo/r-help wrote: * ** I am attempting to extract the derivative/ gradient from this expression ** ** df1p - deriv(f1, P) ** df1p ** expression({ **.value - s - c - a * P **.grad -

Re: [R] approx function wierd result

On Apr 6, 2010, at 10:42 AM, marta rufino wrote: Dear R-list members, I am calculating the linear extrapolation for a data set, using the function found in Hmisc. x=c(0.0265,-0.0003,0.0142,0.0263,0.0634,0.1145,0.2504) y=c(58,107,152,239,362,512,724) x1=0.0393216 approxExtrap(x,y,x1,

Re: [R] glmpath in R

Hi Claire, I'm replying and CC-ing to the R-help list to get more eyes on your question since others will likely have more/better advice, and perhaps someone else in the future will have a similar question, and might find this thread handy. I've removed your specific research aim since that

Re: [R] memory error

someone else on this list may be able to give you a ballpark estimate of how much RAM this merge would require. I don't have an absolute estimate, but try data.table::merge, as it needs less working memory than base::merge. 20 million rows of 5 columns isn't beyond 32bit : (1*4 +

[R] Caret package and lasso

Dear all, I have used following code but everytime I encounter a problem of not having coefficients for all the variables in the predictor set. # code rm(list=ls()) library(caret) # generating response and design matrix X-matrix(rnorm(50*100),nrow=50) y-rnorm(50*1) # Applying caret package

Re: [R] svm of e1071 package

Hi, On Tue, Apr 6, 2010 at 8:07 AM, Shyamasree Saha [shs] s...@aber.ac.uk wrote: Hello List, I am having a great trouble using svm function in e1071 package. I have 4gb of data that i want to use to train svm. I am using Amazon cloud, my Amazon Machine Image(AMI) has 34.2 GB of memory. my

Re: [R] approx function wierd result

Thank David, That is completely correct! I have no ideia how come matlab and Octave do not do this... maybe the interp1 function has some kind of trick. .. maybe the algorythm looks for the closest ranked values and not by its true value. What ever... Just sort it out in a not prety-but

Re: [R] approx function wierd result

ok, not really sorted :-(does not work in other values,...and I am not really sure how to come over it x=c(0.0191,0.0088,0.0150,0.0105,0.0430,0.0697,0.1999) y=c(58,107,152,239,362,512,724) x1=0.02220172 kk=y[x==min(x)] approxExtrap(x[ykk],y[ykk],x1, method=linear) approx(x[ykk],y[ykk],x1) plot(x,

Re: [R] approx function wierd result

On Apr 6, 2010, at 12:13 PM, marta rufino wrote: ok, not really sorted :-(does not work in other values,...and I am not really sure how to come over it x=c(0.0191,0.0088,0.0150,0.0105,0.0430,0.0697,0.1999) y=c(58,107,152,239,362,512,724) x1=0.02220172 kk=y[x==min(x)]

[R] Out of sample forecast

Hello I am trying to use predict, but I am having trouble getting out of sample predictions. I am getting the same output if I use the following three commands: predict(ModeloLineal,predictors[721:768,]) predict(ModeloLineal,predictors[1:768,]) predict(ModeloLineal) where ModeloLineal is the

[R] help in function in R akin to macro in SAS

Dear Whom it may concern, I need help to figure the macro function in R: I need to plot the different data sets by a plotxyf function, I want the title to be different for different data set. # get the data set final.xyf- xyf(data=as.matrix(my.final), Y=classvec2classmat(final$outcome),

Re: [R] Out of sample forecast

On Apr 6, 2010, at 12:25 PM, Luis Felipe Parra wrote: Hello I am trying to use predict, but I am having trouble getting out of sample predictions. I am getting the same output if I use the following three commands: predict(ModeloLineal,predictors[721:768,])

Re: [R] help in function in R akin to macro in SAS

On Tue, 2010-04-06 at 09:33 -0700, Changbin Du wrote: #plot function plotxyf-function (data) { plot(data, type=mapping, labels=final$target, col=final$outcome+1,, main=Supervised: Mapping plot for data) Change this to: main=paste(Supervised: Mapping plot for, data)

Re: [R] glmpath in R

Steve Lianoglou mailinglist.honeypot at gmail.com writes: Hi Claire, I'm replying and CC-ing to the R-help list to get more eyes on your question since others will likely have more/better advice, and perhaps someone else in the future will have a similar question, and might find this

Re: [R] median of two groups

None of them. - mood.test() looks promising until you read the help page and see that it does not do Mood's test for equality of quantiles, it does Mood's test for equality of scale parameters. - wilcox.test() is not a test for equal medians - ks.test() is not a test for equal medians.

Re: [R] help in function in R akin to macro in SAS

A small correction. On Tue, Apr 6, 2010 at 9:43 AM, Stuart Luppescu s...@ccsr.uchicago.eduwrote: On Tue, 2010-04-06 at 09:33 -0700, Changbin Du wrote: #plot function plotxyf-function (data) { plot(data, type=mapping, labels=final$target, col=final$outcome+1,, main=Supervised: Mapping

Re: [R] help in function in R akin to macro in SAS

On Tue, 6 Apr 2010, Stuart Luppescu wrote: On Tue, 2010-04-06 at 09:33 -0700, Changbin Du wrote: #plot function plotxyf-function (data) { plot(data, type=mapping, labels=final$target, col=final$outcome+1,, main=Supervised: Mapping plot for data) Change this to: main=paste(Supervised:

Re: [R] SAS and R on multiple operating systems

Hi, First a quick note about FDA and Style Guides. Frank is correct the FDA does not require strict formatting. I once worked for small startup and they sent EXCEL sheets to the FDA. However big pharma usually has very strict formatting requirements. I once produced a table where one of the

Re: [R] c(), or cat(), or paste(), all cause unwanted reordering

Jeff Brown wrote: Wow, you guys are awesome. Thanks! Nikos Alexandris wrote: Thanks for the cat() question Jeff and to all guRus out there for the replies. This is something I was looking for the last hour. I can't seem to make this run: I have a function ( Column.of.Matrix.1 ,

Re: [R] svm of e1071 package

I think the problem is that you have R configured as 32-bits. If that is the case, then you will only have access to 4 gigs of RAM (see http://www.brianmadden.com/blogs/brianmadden/archive/2004/02/19/the-4gb-windows-memory-limit-what-does-it-really-mean.aspx). Try booting up an ubuntu instance in

Re: [R] library sets: A EMPTY does not work; gset_intersection(A,EMPTY) works

Hi Peter, This looks like another one: gset_is_equal(X,X) and X == X evaluate to FALSE ? uv - c('a', 'b', 'c'); s - gset(uv, runif(length(uv))) s {a [0.0811552], b [0.3552998], c [0.996772]} gset_is_equal(s, s) [1] FALSE s == s [1] FALSE class(s) [1] gset cset Best regards, Ryszard

[R] plot.png

I've shared a document with you: plot.png http://docs.google.com/leaf?id=0B_J5_yv8GrERYWRhZjU4NDItMmMyOS00ZGIxLWE2ZDUtMDFiYTY3MzliMzhkhl=eninvite=CI2a3aYC It's not an attachment -- it's stored online at Google Docs. To open this document, just click the link above. [[alternative HTML

Re: [R] c(), or cat(), or paste(), all cause unwanted reordering

On Tue, 2010-04-06 at 19:29 +0200, Nikos Alexandris wrote: Jeff Brown wrote: Wow, you guys are awesome. Thanks! Nikos Alexandris wrote: Thanks for the cat() question Jeff and to all guRus out there for the replies. This is something I was looking for the last hour. I can't seem to

Re: [R] SAS and R on multiple operating systems

On Tue, Apr 6, 2010 at 6:27 PM, Roger DeAngelis(xlr82sas) rdean...@amgen.com wrote: Hi,  First a quick note about FDA and Style Guides.  Frank is correct the FDA does not require strict formatting. I once worked for small startup and they sent EXCEL sheets to the FDA. However big pharma

Re: [R] SAS and R on multiple operating systems

Hi, Just a couple of notes to reply to others questions. 1. unx is a file handle that points to my UNIX SAS work directory(like an R workspace). This directory is not mounted under windows. most of the SAS objects, ie SAS datasets in the 64bit Unix directory are available to my windows SAS

[R] plot.png

I've shared a document with you: plot.png http://docs.google.com/leaf?id=0B_J5_yv8GrERYWRhZjU4NDItMmMyOS00ZGIxLWE2ZDUtMDFiYTY3MzliMzhkhl=eninvite=CIa4nqkK It's not an attachment -- it's stored online at Google Docs. To open this document, just click the link above. Hi there, Using the code

Re: [R] SAS and R on multiple operating systems

Hi, I agree with Barry, but there are some of us out there that are trying to change the image of big Pharma. R helps statisticians do real science. I have complained for about 10 years about EXCEL and WORD driving the delivery of inputs and outputs. Years ago, I did a short gig with

[R] Need help on predictive modeling of count data

Hi, I have a data set of a number of compounds and their C13 chemical shifts data. The classes of the compounds are either strong or weak. There are 100 feature variables which are different C13 bins. In each C13 bin the values are counts (number of times a chemical shift is present in

Re: [R] Scope and assignment: baffling

Thanks, Bill! That was deep, and took me a long time to work through, but I get it now. And Gabor -- r-proto is great! -- View this message in context: http://n4.nabble.com/Scope-and-assignment-baffling-tp1747582p1753321.html Sent from the R help mailing list archive at Nabble.com.

[R] Adding labels on maps (using sp, maptools, or something else suggested)

Hello, I would like to place region names (or abbreviations if space is too limited) on a map. I have tried the following, copying and adapting various examples I have seen but without success. The first example is self-contained, but the second one requires the .shp file, which can retrieved

Re: [R] Data manipulation problem

On Apr 6, 2010, at 9:56 AM, moleps islon wrote: OK... next question.. Which is still a data manipulation problem so I believe the heading is still OK. ##So now I read my population data from excel. No, you read it from a text file and providing the first ten lines of that text file should

Re: [R] Data manipulation problem

On Apr 6, 2010, at 3:30 PM, David Winsemius wrote: On Apr 6, 2010, at 9:56 AM, moleps islon wrote: OK... next question.. Which is still a data manipulation problem so I believe the heading is still OK. ##So now I read my population data from excel. No, you read it from a text file and

[R] Merging data frames on two conditions

Hi Guys I have two data frames which I would like to merge on two conditions. I am doing the following (abstract form) new.data.frame - merge(df1,df2, by=c(Col1,Col2)) It is giving me a null result. Basically I need to apply two conditions. I also tried sqldf but it is running forever. Will

[R] Error bar

Hi there, Using the code below I generated the plot attached. The error bars at day = -3, 2 and 4 appear are larger. I was thinking if there is a way I could make all the error bars to be of the same size. I don't know if that makes sense. If not, then, is there a way I can plot only these 3

Re: [R] NMDS Ordination Graphics Problem

Hi Trey, On Mon, Apr 5, 2010 at 1:58 PM, Trey trey3...@hotmail.com wrote: Dr. Stevens, Hi, my name is Trey Scott, and I'm a grad student of Brian McCarthy's.  He referred me to you because of your expertise in handling complex R problems. We were hoping you could help us solve a nagging

[R] Need Help with Manipulating Columns

I have uploaded the property data from Lee County Property Appraiser into R, but my problem is that I am unable manipulate and use many of the columns because they are classified as factors. I think this is because the values represent selling prices and therefore have $ in front of them. I

Re: [R] Merging data frames on two conditions

On Apr 6, 2010, at 3:54 PM, Abhishek Pratap wrote: Hi Guys I have two data frames which I would like to merge on two conditions. I am doing the following (abstract form) new.data.frame - merge(df1,df2, by=c(Col1,Col2)) What does str(df1) ; str(df2) ... show? It is giving me a null

Re: [R] Merging data frames on two conditions

Hi David Here it is. You can ignore the bio jargon if it sounds confusing. The corresponding data type of column (SNP, chr) on which I am applying merge is same. merge(data_lane6_snps, data_lane6_snps_rsid , by = c(SNP,chr)) str(data_lane6_snps) 'data.frame': 7724462 obs. of 10 variables:

Re: [R] Merging data frames on two conditions

And I should also add that if I merge only on one column it works fine but the result is not what I want. merge(data_lane6_snps, data_lane6_snps_rsid , by = c(SNP) : works as expected. Is the chr column being a factor creating probs here ? -A On Tue, Apr 6, 2010 at 4:03 PM, Abhishek Pratap

Re: [R] Need Help with Manipulating Columns

possible hint: as.numeric(gsub(\\$, , c($1.23, $1.32))) cavalier33901 wrote: I have uploaded the property data from Lee County Property Appraiser into R, but my problem is that I am unable manipulate and use many of the columns because they are classified as factors. I think this is because

Re: [R] Need Help with Manipulating Columns

On 06/04/2010 3:40 PM, cavalier33901 wrote: I have uploaded the property data from Lee County Property Appraiser into R, but my problem is that I am unable manipulate and use many of the columns because they are classified as factors. I think this is because the values represent selling prices

Re: [R] Need Help with Manipulating Columns

Erik Iverson wrote: possible hint: as.numeric(gsub(\\$, , c($1.23, $1.32))) and convert to character first before as.numeric if they were factors to begin with... cavalier33901 wrote: I have uploaded the property data from Lee County Property Appraiser into R, but my problem is that

Re: [R] svm of e1071 package

Hi Shyama, Don't forget to CC the r-help list in your discussions so that there are more eyes on this problem, and others might potentially benefit from discussion. Comments in line. On Tue, Apr 6, 2010 at 4:06 PM, Shyamasree Saha [shs] s...@aber.ac.uk wrote: Dear Steve, Thanks a lot for

Re: [R] Merging data frames on two conditions

On Apr 6, 2010, at 4:03 PM, Abhishek Pratap wrote: Hi David Here it is. You can ignore the bio jargon if it sounds confusing. Sometimes it is essential to have domain details. The corresponding data type of column (SNP, chr) on which I am applying merge is same. merge(data_lane6_snps,

Re: [R] Merging data frames on two conditions

Hi David I can understand looking the SNP data values it can be felt that they are different values and hence no result in merge. However the columns still have ~700K SNPs common. What I am looking for is a merge where the SNP and Chr matches. If I match only the SNP column I get partially

Re: [R] Merging data frames on two conditions

Just so you know length(intersect(data_lane6_snps$SNP, data_lane6_snps_rsid$SNP)) 796120 I just need to include the chr condition now where I am stuck. -Abhi On Tue, Apr 6, 2010 at 4:51 PM, Abhishek Pratap abhishek@gmail.comwrote: Hi David I can understand looking the SNP data values

Re: [R] lattice package: line end style

On 2010-04-06 7:28, Daniel Alcock wrote: First, apologies for no example data but I don't think it's needed in this case, Q: Can (and if so how ) the line end style be changed for 'cloud' plots? I've tried par(lend=2), trellis.par.set(add.line = list(lend=2)) and much googling but to no avail

Re: [R] Merging data frames on two conditions

OK, not the SNP's. So look at the chr's. I will bet that you get 0 when you try : length(intersect(data_lane6_snps$chr, data_lane6_snps_rsid$chr)) ... since one is using a format of chrNN and the other is using just NN. You need to get the chromosome naming convention straightened out.

Re: [R] checking bivariate normality

On Tue, Apr 6, 2010 at 12:48 PM, arindam fadikar arindam.fadi...@gmail.com wrote: x - iris$Sepal.Length[1:50]/iris$Sepal.Width[1:50] y - iris$Petal.Length[1:50]/iris$Petal.Width[1:50] I want to check whether (x,y) follows a bivariate normal distribution or not,  using density plot or scatter

Re: [R] Merging data frames on two conditions

You got the error. It is different naming convention of chr. I should be able to fix that pretty easily. In case the problem persists, I will contact the list. Thanks! -Abhi On Tue, Apr 6, 2010 at 5:01 PM, David Winsemius dwinsem...@comcast.netwrote: OK, not the SNP's. So look at the chr's. I

Re: [R] checking bivariate normality

Below. -- Bert Bert Gunter Genentech Nonclinical Statistics -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Liviu Andronic Sent: Tuesday, April 06, 2010 2:06 PM To: arindam fadikar Cc: r-help@r-project.org Subject: Re: [R]

[R] How to import data from teradata

Hello All, I am very new to the R and i want to know how to connect to the teradata in R and import data from there. Currently I am using SAS to import data and i use the following syntax to connect to teradata. *Proc Sql; connect to odbc (dsn=teradata username = '**' password =

[R] appending an R-object to a list

How do I append an R-object to a list? I want to start with an empty list, and append R-objects one by one. Does this start with a command like mylist - NULL ?? I have read a few answers on R-help to questions like this, but they all seem to be well off the point. Sometimes it's assumed that the

[R] estimating the starting value within a ODE using nls and lsoda

All- I am interested in estimating a parameter that is the starting value for an ODE model. That is, in the typical combined fitting procedure using nls and lsoda (alternatively rk4), I first defined the ODE model: minmod - function(t, y, parms) { G - y[1] X - y[2]

Re: [R] Need Help with Manipulating Columns

Nice, thanks A lot. I got it figured out -- View this message in context: http://n4.nabble.com/Need-Help-with-Manipulating-Columns-tp1753364p1753471.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

Re: [R] read numeric values with thousands seperator from csv file

Thanks I tried it with your values and it works, but when I try it with my values, plugging in just (not sure if i'm doing it right): as.numeric(gsub(\\$, , (just))) This is what I get: [1] NA Warning message: NAs introduced by coercion Btw here are the values im wokring with under

Re: [R] appending an R-object to a list

On Tue, Apr 6, 2010 at 10:28 PM, David.Epstein david.epst...@warwick.ac.uk wrote: How do I append an R-object to a list? I want to start with an empty list, and append R-objects one by one. Does this start with a command like mylist - NULL ?? I have read a few answers on R-help to

Re: [R] Odd results with %% and conserving memory

Thanks for the response. I read through the Inferno and there's a lot of useful stuff in there. However I think you may have missed something. Having some experience with other languages, I know that floating point numbers aren't represented exactly, the issue is that R is off by much more than

Re: [R] How to import data from teradata

On Apr 6, 2010, at 4:26 PM, sumit gupta wrote: Hello All, I am very new to the R and i want to know how to connect to the teradata in R and import data from there. Currently I am using SAS to import data and i use the following syntax to connect to teradata. *Proc Sql; connect to

Re: [R] Odd results with %% and conserving memory

Ummm... I think the problem is with _your_ inadequate understanding. Consider: identical( (4.1 * 100 ) %/% 1 , 409 ) 4.1*100 410 [1] TRUE So, due to finite precision, 4.1*100 is slightly less than 410. Whence the largest integer 1 can be divided into it is 409. i.e. (4.1*100) %/% 1 [1] 409

  1   2   >