Re: [R] Problem in r help me

From: Ghada Almousa [mailto:ghada.f...@gmail.com] Sent: Monday, May 04, 2015 9:20 PM To: PIKAL Petr Subject: Re: Problem in r help me I used internet search to get info about methods but I can't find any info about the comparison between those methods The comparison in the number of iterations ,

Re: [R] GSVD and qr questions

off line I looked at the matlab page and the wikipedia page. Wikipedia says there are two different things with the gsvd name. So I guess we need to see a tiny example of your starting matrix and the decomposition that the matlab function you call returns to you. I guess you will have to send the

Re: [R] GSVD and qr questions

Thanks. But to be honest I am outside of what I really know. so I don’t know if what you describe is what I want I know what the Matlab function I have calls and what it does, and I am trying to mimic its functionality. So I am trying to see if there is equivalent to the Matlab qr(A,0) call,

Re: [R] GSVD and qr questions

Roy, Look at the qr.R and qr.Q functions. With the complete=TRUE argument they bring their matrices up to full size, by constructing an arbitrary orthogonal completion. Does this provide what you are looking for? Rich On Mon, May 4, 2015 at 11:56 PM, Roy Mendelssohn - NOAA Federal wrote: > Hi

[R] using filter() to sum up

Hi everybody, I am trying to create a code for the formula in the attachment. I first tried following code: ltau <- m + theta*sum(psi*X[t-k]) but it does not work and I get for X[t-k] every third element in my vector three times which looks as follows: X[t-k] [1] -0.25 -0.25 -0.25 0.50 0.50 0.

[R] Rhive queries failing

Hi folks, I am facing an error while firing Rhive queries. Its able to connect via hive server2, list tables do a simple select * operation. But, a count operation is failing or if I do select "field" from table_name, that fails too. Any query that involves processing (like in MR terms) is faili

[R] GSVD and qr questions

Hi All: I posted some of this earlier to the r-sig-mac because at that time I got a Mac specific error, but now I have a more general question. For a canonical variate analysis (CVA) used in some system identification routines, I need to calculate a generalized svd (GSVD). On the web there ar

Re: [R] Fwd: Question about paired plotting

Hi With R there are usually several ways Here is another one using lattice (I think it ships with R otherwise install) Reshape dat1 to long form and add a column for groups do give d2 d2 = structure(list(y = c(5L, 10L, 6L, 11L, 7L, 12L, 8L, 13L), x = c(1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L), gp = str

Re: [R] A problem with string handling to make a time duration

Regular expressions are the tool for this problem. This pattern matches your input data: t <- c("10H20M33S", "1H1M", "1M", "21M9S", "2H55S") patt <- "^(([0-9]+)H)?(([0-9]+)M)?(([0-9]+)S)?$" all(grepl(patt, t)) # TRUE We can use the pattern to extract hour/minute/second components hms <- lapply(c

Re: [R] Graphs for scientific publication ?

On May 4, 2015, at 3:13 AM, Jeremy Clark wrote: > Dear All, > > Many thanks for your very comprehensive replies. Here I provide some > coding which on my system has the following effects: > 1) The italic R is not rendered by CairoX11, but is rendered by quartz. > 2) Both geom_smooth and geom_abl

Re: [R] Graphs for scientific publication ?

Hi Jeremy, On Mon, May 4, 2015 at 6:13 AM, Jeremy Clark wrote: > Dear All, > > Many thanks for your very comprehensive replies. Here I provide some > coding which on my system has the following effects: > 1) The italic R is not rendered by CairoX11, but is rendered by quartz. I don't have a Mac

[R] [R-pkgs] Version 3.16-0 of Hmisc now on CRAN

Several updates have been made to the Hmisc package: Changes in version 3.16-0 (2015-04-25) * html.contents.data.frame: corrected html space character to add semicolon * ggplot.summaryP: added size of points according to denominators * colorFacet: new function * labelPlotmath: add

Re: [R] Need Help!

Dear Dawood, it is not a really good idea to open a new email with the same problem you already posted earlier. I see you accepted some suggestions (HTML, seq) but not others (pi is a defined constant in R, do not change it, type pi ENTER to see this.). I would like to add one more, c is a func

[R] A problem with string handling to make a time duration

I have a character string that represents a time duration. It has an hours minutes seconds structure(ish) but with letters denoting units (H,M or S) no leading zeros and no placeholder at all where one or other of the units are not required. It looks like this: t<-c("10H20M33S","1H1M","1M","21M9S

Re: [R] Mantel test

Assuming the 'matrix' format is a symmetrical distance 'matrix' stored as a data frame (which read.csv creates) rather a rectangular data 'matrix,' you can convert it to a dist object with as.dist(). ?dist - David L Carlson Department of Anthropology Texas A

Re: [R] Fwd: Question about paired plotting

Dear John and Tim! Thanks a lot for your replies. This is just what I needed!! All the best for you!! 2015-05-04 13:58 GMT-03:00 John Kane : > Hi Luis, > Ah yes, that paper. I was rather shocked at what it implied :( People > are still using dynamite plots, for heaven's sake! See > http://bio

[R] Block exogeneity in SVAR

Good day, My goal is to use a SVAR model to estimate monetary policy transmission. I have seven variables, five for a small country and two for a large country. I�d like to impose block-exogeneity restrictions so that the small county does not affect the large country (in the short run nor in the

[R] ggplot, Plot title at the bottom of plot

Hello All, I am looking for the answer to simple question: In ggplot, how to put plot title at the bottom of plot ? This code : qplot(rnorm(100)) + ggtitle("My Title") puts the title at middle and top of plot while I want it to be at middle and *Bottom* of the plot. Using +theme(plot.title =

Re: [R] GLM: What is a good way for dealing with new factor levels in the test set?

For anyone who is looking for an answer to this in the future... I went for "imputation". It's a way of filling in missing variables based off of what you see elsewhere in the data. Myself, I simply took a sample of the categorical from the rest of the test set. Some may argue that this is erro

[R] Mantel test

Dear R users, I'm having trouble getting my data into R in the correct format to run a Mantel test. I'm testing genome size differences by genetic distances of the 28S gene for ~30 species. I'm able to get my genome size data (as a single column of data) into matrix and dist formats in R but the

Re: [R] Request for functions to calculate correlated factors influencing an outcome.

Hi, >From my understanding of a model, you need to have one or more dependent parameter variables which you (y1, y2 etc) and the explanatory parameter variables (x1, x2 etc) which you fit along with certain coefficients to determine the dependent parameter with minimum error as possible. I don't t

Re: [R] Graphs for scientific publication ?

Dear All, Many thanks for your very comprehensive replies. Here I provide some coding which on my system has the following effects: 1) The italic R is not rendered by CairoX11, but is rendered by quartz. 2) Both geom_smooth and geom_abline here give stepped lines (I've realised the angle of the li

Re: [R] Fwd: Question about paired plotting

Hi Luis, Ah yes, that paper. I was rather shocked at what it implied :( People are still using dynamite plots, for heaven's sake! See http://biostat.mc.vanderbilt.edu/wiki/Main/DynamitePlots for some comments. I see Jim has given you one way to do the plots you want. Here is another way using

Re: [R] Downloading file.Rdata from internet

The answer to this question does not actually involve R. Directly accessing data files is not allowed by that web service. I suspect the reason is that they need to get paid somehow and they don't want you to bypass their ads. You may have better success if you use a paid web server to share you

Re: [R] Downloading file.Rdata from internet

Dear David and Robert, Thank you for your comments. In fact, I created the file. I would just like to load the file directly from the internet to R. Is there a command to do the R uncheck the "Use our download manager and get recommended downloads" and click the "DOWNLOAD" button? Thanks for you

Re: [R] gdata library 2.16.1

On 5/4/2015 9:01 AM, Marc Girondot wrote: Dear list-members, Since I update gdata library to 2.16.1 version this morning, I have an error on the two macs I use (details on system and R versions at the end). When I load the package, I have this error: > library("gdata", lib.loc="/Library/

Re: [R] Downloading file.Rdata from internet

On 5/3/2015 7:58 PM, Rafael Costa wrote: > Dear R users, > > To load the file into "http://www.datafilehost.com/d/c7f0d342";, I first > uncheck the "Use our download manager and get recommended downloads" option > and I click the "DOWNLOAD" button. How do I load and save the file directly > from

[R] gdata library 2.16.1

Dear list-members, Since I update gdata library to 2.16.1 version this morning, I have an error on the two macs I use (details on system and R versions at the end). When I load the package, I have this error: > library("gdata", lib.loc="/Library/Frameworks/R.framework/Versions/3.2/Resources/

Re: [R] Request for functions to calculate correlated factors influencing an outcome.

This would be better posted on a statistical list like stats.stackexchange.com, as it is largely about statistical methodology, not R code. Once you have determined what kinds of methods you want, you might then post back here -- or better yet, just search! -- for packages that implement those meth

Re: [R] [R-SIG-Mac] Unsubscribe from the list

Hi Silvia Solis, Re: > Hi, > > I wanted to ask you to scratch this email address from the mailing > list. This is a work address and I already receive your emails con my > personal account. My mistake this one time I replied from a different > email address. > > Thanks! > > Silvia > -- >

Re: [R] help - hoslem.test

Hi thanks for data. na.exclude does not excludes NA values. Actually it computes result without considering NAs but keeps those NA values in propper positions so that the result has the same length as input Instead of several ifelse you can use q131 <- as.numeric(cut(id3$q13,c(0,1.5,5)))-1 f

Re: [R] Problem in r help me

Hi Your question is rather cryptic. Data? Code? You can use internet search to get info about methods e.g R hierarchical clustering gives me many answer regarding clustering in R. Be more specific to get more specific answer Cheers Petr > -Original Message- > From: R-help [mailto:r

Re: [R] Request for functions to calculate correlated factors influencing an outcome.

Hi I used the MASS library library(MASS) (by reading about examples at http://www.statmethods.net/stats/regression.html ) fit <- lm(Mileage~Disp+HP+Weight+Reliability,data=newx) step <- stepAIC(fit, direction="both"

Re: [R] Fwd: Question about paired plotting

Hi Luis, Two plotting functions from the plotrix package might be of interest. One is "ladderplot" that is very similar to the upper plots in the panels, and the second is "ehplot" that is similar to the lower panels. Jim On Mon, May 4, 2015 at 5:52 PM, Luis Fernando García wrote: > Hello R exp

Re: [R] Downloading file.Rdata from internet

On May 3, 2015, at 5:58 PM, Rafael Costa wrote: > Dear R users, > > To load the file into "http://www.datafilehost.com/d/c7f0d342";, I first > uncheck the "Use our download manager and get recommended downloads" option > and I click the "DOWNLOAD" button. How do I load and save the file directly

[R] Fwd: Question about paired plotting

Hello R experts, I just found a new paper which shows the proper way (according to the authors) to show data, specially paired. I am very interested in presenting this kind of data, specially the scatter plott. I have found a way to present it using this link: http://journals.plos.org/plosbiolog