Re: [R] Sort data and symbol change in Jitter plot (ggplot2)

Thanks all of you guys! Your answers were very useful! Thanks for the answer John, but I will try to keep the dotplot, but it is very useful to know both techniques anyway :) Petr, many thanks for your help, it was just what I needed, btw, the arrangment based on the median was a headache for me

Re: [R] Not able to submit work assignment4 with submitscript3.R

Hi there, this is not a Coursera Forum. Please ask your question there or search the Forums, these are quite regular problems there. PS: do you really have a user name and a password in your e-mail? Feladó: R-help [r-help-boun...@r-project.org] ; meghatalm

[R] nlminb supplying NaN parameters to objective function

Dear list, I am doing some maximum likelihood estimation using nlminb() with box-constraints to ensure that all parameters are positive. However, nlminb() is behaving strangely and seems to supply NaN as parameters to my objective function (confirmed using browser()) and output the following: $pa

[R] Not able to submit work assignment4 with submitscript3.R

Hi, I tried to submit the work assignment4 with submitscript3.R and i got an error message: Error in assign(".CourseraLogin", r, globalenv()) : unused arguments (r, globalenv()) I tryied to do it manually, but the files that I got to upload was empty. Please, let me know, how can i do to send th

[R] How to finding a given length of runs in a series of data?

I'm trying to study times in which flow was operating at a given level or greater. To do so I have created a way to see how long the series has operated at a high level. But for some reason the data is calculating the runs one hour to long. Any ideas on why? Code: Date<-format(seq(as.POSIXct(

[R] Convert csv to xpt file in R?

Dear Rxperts.. Was wondering if there is a way in R to read a csv file and generate an XPT file? For some reason the function write.xport() does not seem to work for me i get the following error... "error in label.data.frame(df,default=""): length of default same as x A sample dataframe is g

[R] package implementing continuous binomial?

Hi I'm wondering if anyone is aware of an R package implementing (i.e., providing a pdf, cdf, and/or quantile function) for the continuous binomial distribution? Specifically the one characterized here: http://www2.math.uni-paderborn.de/fileadmin/Mathematik/AG-Indlekofer/Workshop/Satellite_meetin

Re: [R] lm model exported from R to excel

At the expense of extra bandwidth. I also like latex + booktabs; hate word I format the table header and footer in latex with extra vertical space between lines and type; output looks more like latex than being cramped from xtable. Thanks anyway to xtable creators -- would be lost without it Dun

Re: [R] Drawing the regression line and the 95% confidence intervals

Jim, Perfect Many thanks once more for your support... Best regards SV Envoyé de mon iPhone Le 7 mai 2015 à 00:01, Jim Lemon a écrit : > Hi Sacha, > As Brian noted, you will have to expand the vertical axis. Here is one > way to do it. > > plot(GDP.per.head, fitted(LinearModel.1),ylim=c(-60,1

Re: [R] Drawing the regression line and the 95% confidence intervals

Hi Sacha, As Brian noted, you will have to expand the vertical axis. Here is one way to do it. plot(GDP.per.head, fitted(LinearModel.1),ylim=c(-60,1200)) devlm1<-lm(fitted(LinearModel.1)~GDP.per.head) abline(devlm1) conflm1<-confint(devlm1) abline(coef=conflm1[,1],lty=2) abline(coef=conflm1[,2],lt

Re: [R] Drawing the regression line and the 95% confidence intervals

Brian, thanks for the precisions. David, many thanks for your code that perfectly works.  Best, Sacha De : David L Carlson À : "Cade, Brian" ; varin sacha Cc : R-help Mailing List Envoyé le : Mercredi 6 mai 2015 21h52 Objet : RE: [R] Drawing the regression line and the 95% confidence

Re: [R] order by which the eigenvalues are presented

On 06/05/2015 12:58 PM, Luis Borda de Agua wrote: > Thank you, Bill, for your reply. However, I'm afraid I didn't explain myself > properly. > > > > Imagine you have a 2x2 matrix > > Then the eigenvalues lambda_1 and lambda_2 are analytically calculated from > > > > lambda_1 = (-b+sqrt(de

Re: [R] Drawing the regression line and the 95% confidence intervals

Something like this? # Compute the prediction limits and get their range to set ylim= plim <- predict(LinearModel.1, interval = "prediction") rnge <- c(min(plim[ , 2]), max(plim[ , 3])) plot(GDP.per.head, fitted(LinearModel.1),ylim=rnge) # As before devlm1<-lm(fitted(LinearModel.1)~GDP.per.head)

Re: [R] Drawing the regression line and the 95% confidence intervals

The prediction intervals are likely to be much wider than the confidence intervals so you will need to be sure you scale the yaxis limits large enough to see them. Brian Brian S. Cade, PhD U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818

Re: [R] Package Build Recommendations

On 05/05/2015 8:16 PM, Glenn Schultz wrote: Hi All, I have my R package built and it passes the CRAN tests. Now, I have a question. The file structure is not standard. There are addition data files as follows outlined below. Each, if you will, represents separation of concerns with respec

Re: [R] Drawing the regression line and the 95% confidence intervals

Dear Jim, I really thank you lots, it perfectly works ! The reproducible example is below. Last thing : I can easily get the predictions intervals but I don't get to draw them. If I want to draw on the same graph (the one I already have the confidence intervals) the prediction intervals. I hav

Re: [R] Special character is graph label

Also have a look at http://astrostatistics.psu.edu/su07/R/html/grDevices/html/plotmath.html for a list of terms. Oh and it's called Xbar or X-bar. John Kane Kingston ON Canada > -Original Message- > From: dcarl...@tamu.edu > Sent: Wed, 6 May 2015 16:33:01 + > To: tacsun...@yahoo

Re: [R] Package Build Recommendations

Hi Glenn, Generally data files are stored in the 'data' directory. If you visit Hadley's R packages site on the data page (http://r-pkgs.had.co.nz/data.html) this is described quite clearly. You can use the devtools package functions like `use_data` to have data files properly stored in your pac

[R] order by which the eigenvalues are presented

Thank you, Bill, for your reply. However, I'm afraid I didn't explain myself properly. Imagine you have a 2x2 matrix Then the eigenvalues lambda_1 and lambda_2 are analytically calculated from lambda_1 = (-b+sqrt(delta))/2a lambda_2 = (-b-sqrt(delta))/2a where delta = b^2-4ac If

Re: [R] does segfault mean (always) a bug?

It looks like a problem in the Matrix package. I made the file KE.rda containing the Matrix objects K and E constructed in calc.diffusion.kernel by adding a save() call just before where R dies in the original example: K = lam$values[1] * I - M E = I - matrix(1, ncol = ncol(I), nr

Re: [R] Special character is graph label

You do this with plotmath (see the manual page, ?plotmath): > plot(rnorm(10), rnorm(10), xlab=expression(bar(X))) - David L Carlson Department of Anthropology Texas A&M University College Station, TX 77840-4352 -Original Message- From: R-help [mailto:

Re: [R] lm model exported from R to excel

If you are trying to save the html version of the table, Duncan's example can be simplified to > print(xtable(mod), type="html", file="lm.html") Open this file in Excel or Word. Alternatively open it in your browser and paste it into Word or Excel. If you are using Windows and Excel, you can s

Re: [R] lm model exported from R to excel

Hi Duncan Thank you so much - it worked :) Best Livia Fra: Duncan Mackay [dulca...@bigpond.com] Sendt: 6. maj 2015 14:26 Til: R; Livia Maria Vestergaard Emne: RE: [R] lm model exported from R to excel Hi Livia There are several html packages that ?coul

[R] Special character is graph label

Dear R users, I am having issues finding a special character (and how to insert it) in the lab of a graph axis. Let us say that the label of my axis is "X", i would like the X to have a "line" over it, indicating that it is the "mean of X values" (i don't even know how to properly state that i

Re: [R] lm model exported from R to excel

And for those of us who know close to nothing about HTML I found just now that under a basic print.xtable commmand we get those horrible HMTL borders that in Apache OpenOffice seemed impossible to remove safely. No idea about Word--I have not used it in years. I did find that adding html.tabl

Re: [R] Looking up the code for a function

Here's a discussion on stack overflow with more hints ... http://stackoverflow.com/questions/19226816/how-can-i-view-the-source-code-for-a-function ... and a link to an R-News article by Uwe Ligges on the topic http://cran.r-project.org/doc/Rnews/Rnews_2006-4.pdf (p. 43 f) Cheers, Boris On Ma

Re: [R] Looking up the code for a function

Probably these are non-exported functions. Try: getAnywhere() Or if you know which package a function comes from: ::: Best, Wolfgang -- Wolfgang Viechtbauer, Ph.D., Statistician Department of Psychiatry and Neuropsychology School for Mental Health and Neuroscience Faculty of He

[R] Looking up the code for a function

Hi all, I am trying to find out how a certain functionality is implemented in R respectively what a certain found does exactly. Specifically I am interested in multivariate kernel density estimation. I found the "ks" package and its "kde" function. Usually, my preferred way to "look under the hoo

Re: [R] Mantel test

Please keep the discussion on the list. It is hard to answer your question without context. You give us a warning message from a function that is not in base R. Perhaps it is in package ade? You also don’t include anything about the data or the commands that produced the warning message. A dist

Re: [R] lm model exported from R to excel

If you know some basic html language you can jazz up the table headings to your liking by writing that before the xtable statement. It save having to muck around in Microsoft to fix it. If you are going to do a lot of it - a little study of html basics can go far. I was changing the headings to wh

Re: [R] MOnth over Month Variance in %

I think we need to know a bit about your level of knowledge of R, what the data currently looks like (In R format preferably) and maybe see what you have tried so far. Please have a look at the following link, and pay particular attention to the use of dput() as a way to present data on R-help.

Re: [R] MOnth over Month Variance in %

It seems like the ?diff function might be useful here, but without a reproducible example (at the very least some data provided as dput output and results of your Excel calculations for that data) we would have to guess quite a bit and likely be wrong. Take the recommendations at [1] to heart

[R] Decomposing irregular time series

​Dear fellow members of the r-help list, ​May someone be so kind as to point out​ how an irregular time series should be decomposed into season and trend? My code, the related error and my data are enclosed herein. ​Best regards and thanks in advance! ​-- Jue​ > x<-read.table(file="Data3.txt"

Re: [R] lm model exported from R to excel

Hi Livia There are several html packages that ?could also do it Heres a way with xtable library(xtable) y = rnorm(100) x= rnorm(100)+rnorm(100) mod <- lm(y ~x) # latex example easy view xtable(mod) # html file.create("lm.htm") ff <- file("lm.htm", "a+") fchars <- print(xtable(mod),type = "ht

Re: [R] lm model exported from R to excel

Hi AFAIK result of lm model is a list so you cannot expect easy direct output to Excel. List of 12 $ coefficients : Named num 61.2 ..- attr(*, "names")= chr "cas" $ residuals: Named num [1:2] 0 0 ..- attr(*, "names")= chr [1:2] "1" "2" $ effects : Named num [1:2] -61.2 0 ..- at

Re: [R] lm model exported from R to excel

Hi Livia, One way is to use the "delim.table" function in the prettyR package. See the examples, in particular the final one. The resulting TAB delimited file will usually import directly into Excel. Jim On Wed, May 6, 2015 at 7:36 PM, Livia Maria Vestergaard wrote: > Hi all > I all. I am wonde

[R] lm model exported from R to excel

Hi all I all. I am wondering whether anybody know how to export an output of an lm model from R to excel in order to have excel recognize the table that comes and divide the numbers in the table into columns and rows? I really hope it is possible? :) Best Livia ___

[R] MOnth over Month Variance in %

Hi All, I have data based on truck load for various states. The data points range from Oct'14 To Mar'15. Now I need to know what was the difference in load in Nov as compared to Oct in both real numbers as well as in %. Similarly for all the month in comparison to the previous month. I am able to

Re: [R] Plot by FIPS Code using Shapefiles

Dear Anthony, Thanks again for your reply. The following worked for merge: merged_data <- merge(shapes_fips,max_change,by="FIPS",all.x=T, all.y=F) However, I think I am doing something wrong - as I have 3109 FIPS code in my original data but when I merge with the shapes file SpatialPolygonsDat

Re: [R] Plot by FIPS Code using Shapefiles

Corey Sparks utsa.edu> writes: > > Joining data the way you're doing it is dangerous, Roger Bivand and others > describes a standard way to do this process here: > http://r-sig-geo.2731867.n2.nabble.com/Merging-shapefiles-and-csv-td7586839.html Quite right - the chunks Corey is referring to ar

Re: [R] predict function in regression analysis

Maybe this is what you wanted : Data <- structure(list(y = c(4.5, 4.5, 4.7, 6.7, 6, 4.4, 4.1, 5.3, 4, 4.2, 4.1, 6.4, 5.5, 3.5, 4.6, 4.1, 4.6, 5, 6.2, 5.9, 3.9, 5.3, 6.9, 5.7), lot = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L), duration = c(0

Re: [R] Rgraphviz and NA indices error

Hello, I ran into the same error. From what I understood, if the code : plot(dtm, terms=findFreqTerms(dtm, lowfreq=100) [1:30], corThreshold=0.75) gives you this message it's because there aren't 30 terms which happens at least a hundred times in your corpus. Therefore you mi

[R] R write to JSON format error

Dear all There are a stardand JSON format file from R code. It have a error on this file. How to modify this source code for a stardand JSON format file? Test1: Soure code for R studio: Issue_Name<-c("Tag over 30% in minutes") URL_Link<-c("http://61.2

Re: [R] does segfault mean (always) a bug?

> lejeczek > on Wed, 6 May 2015 08:20:46 +0100 writes: > On 05/05/15 20:36, Duncan Murdoch wrote: >> On 05/05/2015 2:54 PM, lejeczek wrote: >>> hi eveybody >>> >>> I'm trying something simple (Biocunductor packages), so >>> simple I believe it's example from

Re: [R] Need Help!

Hi, type res # (now it contains your results as a vector) or head(res) #so you only see the first couple of values or plot(res) #to see a basic plot (you could set axis labels etc, but I suggets you to go ahead and read a tutorial on R! Best, have a great time with R, kd

Re: [R] Transformation of the circular variable into linear variable

> David Winsemius > on Tue, 5 May 2015 08:41:08 -0700 writes: > On May 5, 2015, at 3:59 AM, Nasr Al-Dhurafi wrote: >> Dear All >> >> I am working with a circular data ( wind direction ) and i have problem >> in transforming the circular variable into linear vari

[R] Package Build Recommendations

Hi All, I have my R package built and it passes the CRAN tests.  Now, I have a question.  The file structure is not standard.  There are addition data files as follows outlined below.  Each, if you will, represents separation of concerns with respect to structured securities like MBS and REMIC

Re: [R] Plot by FIPS Code using Shapefiles

Joining data the way you're doing it is dangerous, Roger Bivand and others describes a standard way to do this process here: http://r-sig-geo.2731867.n2.nabble.com/Merging-shapefiles-and-csv-td7586839.html And I do an example using US Census data here, using merge(): http://spatialdemography.org/w

Re: [R] does segfault mean (always) a bug?

On 06/05/15 09:32, Martin Maechler wrote: Duncan Murdoch on Tue, 5 May 2015 15:36:59 -0400 writes: > On 05/05/2015 2:54 PM, lejeczek wrote: >> hi eveybody >> >> I'm trying something simple (Biocunductor packages), so >> simple I believe it's example from docs but I

Re: [R] does segfault mean (always) a bug?

> Duncan Murdoch > on Tue, 5 May 2015 15:36:59 -0400 writes: > On 05/05/2015 2:54 PM, lejeczek wrote: >> hi eveybody >> >> I'm trying something simple (Biocunductor packages), so >> simple I believe it's example from docs but I get segfault. >> I don't suppos

Re: [R] does segfault mean (always) a bug?

On 05/05/15 20:36, Duncan Murdoch wrote: On 05/05/2015 2:54 PM, lejeczek wrote: hi eveybody I'm trying something simple (Biocunductor packages), so simple I believe it's example from docs but I get segfault. I don't suppose incorrect scripting can cause segfault, right? In R, a segfault always