similar to: number to volume weighted distribution

Hi Duncan I tried Ecdf and/or wtd.quantile from Hmisc and it is working (probably). Ecdf(x, q=.5) Ecdf(x, weights=xw,col=2, add=T, q=.5) wtd.quantile(x) 0% 25% 50% 75% 100% 10 10 10 100 300 wtd.quantile(x, weights=xw, type="i/n") 0% 25% 50% 75% 100% 10.0000 138.8667 192.5778 246.2889 300.0000 But could you please be more specific in this? >

Hi Petr, I think that Duncan suggests something like this: x<- c(rep(10,20), rep(300,5), rep(100, 10)) tx <- table(x) prop.x <- tx / sum(tx) vx <- as.integer(names(tx)) prop.wx <- tx * vx / sum(tx * vx) plot(ecdf(x)) plot(vx, cumsum(prop.x), ylim = 0:1) plot(vx, cumsum(prop.wx), ylim = 0:1) Best regards, Thierry ir. Thierry Onkelinx Statisticus / Statistician Vlaamse

On 24/11/2017 6:27 AM, PIKAL Petr wrote: > Dear all > > Strictly speaking it is not R question but as you are the most capable persons I know I give it a try. > > I am strugling with recalculation of number weighted to volume weighted distribution. > > Suppose I have objects (cubes) with size > > x<- c(rep(10,20), rep(100, 10), rep(300,5)) > I can get >

Hi Great. I did not think that such combination is posssible. Thanks. Petr > -----Original Message----- > From: peter dalgaard [mailto:pdalgd at gmail.com] > Sent: Tuesday, July 11, 2017 1:11 AM > To: PIKAL Petr <petr.pikal at precheza.cz> > Cc: r-help at r-project.org > Subject: Re: [R] fit lognorm to cdf data > > How about > > proc <- c(0.9, 0.84, 0.5,

Hi I copied your values to R, here it is > dput(temp) temp <- structure(list(par1 = structure(1:4, .Label = c("x1", "x2", "x3", "x4"), class = "factor"), y1 = c(-0.19, 0.45, -0.09, -0.16), y2 = c(0.4, -0.75, 0.14, -0.01), y3 = c(-0.06, -8.67, 1.42, 2.21), y4 = c(0.13, -0.46, 0.06, 0.06)), .Names = c("par1", "y1",

Hi > -----Original Message----- > From: Robert Baer [mailto:rbaer at atsu.edu] > Sent: Wednesday, August 16, 2017 3:04 PM > To: PIKAL Petr <petr.pikal at precheza.cz>; Duncan Murdoch > <murdoch.duncan at gmail.com> > Cc: r-devel at r-project.org > Subject: Re: [Rd] strange behaviour read.table and clipboard > > You said, "put a name in the cell".

> On 9 Nov 2017, at 14:58, PIKAL Petr <petr.pikal at precheza.cz> wrote: > > Hi > > Thanks for working example. > > you could use split/ lapply approach, however it is probably not much better than dplyr method. > > sapply(split(mydf, mydf$type), function(speed, n_vehicles) sum(mydf$speed*mydf$n_vehicles)/sum(mydf$n_vehicles)) > gives you averages > The

Hi Peter; The R code you provided works very well. Once again thanks so much for this. The number of variables in my data set that should appear on the y-axis is 733 and they are not numerical (for example the name of one variable is *palmitoyl-arachidonoyl-glycerol (16:0/20:4) [1]**. So, the plot looks very messy in one page. How can I make the plot to print out on multiple pages? Regards,

Hi Keep your messages coppied to R helplist, others could give you answers too. See in line From: Rosa Oliveira [mailto:rosita21 at gmail.com] Sent: Tuesday, August 1, 2017 4:38 PM To: PIKAL Petr <petr.pikal at precheza.cz> Subject: Re: [R] Superscript and subscrib R for legend x-axis and y-axis and colour different subjects in longitudinal data with different colours Hi Pikal, I looked

Hi Thanks for working example. you could use split/ lapply approach, however it is probably not much better than dplyr method. sapply(split(mydf, mydf$type), function(speed, n_vehicles) sum(mydf$speed*mydf$n_vehicles)/sum(mydf$n_vehicles)) gives you averages aggregate(mydf$n_vehicles, list(mydf$type), sum)$x gives you sums Cheers Petr > -----Original Message----- > From: R-help

Hi Maybe there are other ways but I would split data to several chunks e.g. in list and use for cycle to fill multipage pdf. With the toy data something like library(reshape2) library(ggplot2) temp <- melt(temp) temp.s<-split(temp, cut(1:nrow(temp), 2)) pdf("temp.pdf") for (i in 1: length(temp.s)) { p <- ggplot(temp.s[[i]], aes(x=par1, y=variable, size=abs(value),

Hi When you load external file to R, character variables are converted to factors by default and alphabetically sorted. I have limited connection to internet, so I cannot find the answer, you could try it yourself. Maybe you could try not to convert vector with names to factor, which, for plotting issue is not different from factor coding. See ?read.table for details However I am not sure if it

Hi Petr; Thanks so much. Exactly this is what I need. I will play to change color and so on but this backbound is perfect to me. I do appreciate your help and support. Regards, Greg On Thu, Feb 8, 2018 at 1:29 PM, PIKAL Petr <petr.pikal at precheza.cz> wrote: > Hi > > I copied your values to R, here it is > > > > > dput(temp) > > > > temp <-

Hi Petr; Thanks so much. This is great! Although last Sunday, alternatively, I have solved the problem using the following statement at the very end of the program. * ggsave('circle.pdf', p4, height = 70, width = 8, device=pdf, limitsize = F, dpi=300).* This works very well too. Asa my categorical variables are in my Y axis, my R program reorders the names on Y-axis. However, I would

Hi > -----Original Message----- > From: Martin Maechler [mailto:maechler at stat.math.ethz.ch] > Sent: Wednesday, August 16, 2017 12:01 PM > To: PIKAL Petr <petr.pikal at precheza.cz> > Cc: Martin Maechler <maechler at stat.math.ethz.ch>; Duncan Murdoch > <murdoch.duncan at gmail.com> > Subject: RE: [Rd] strange behaviour read.table and clipboard > >

I did the same transition, facing the same issues, so I settled for installation of RKWard. This pretty much reproduces my previous experience with Tinn-R and the Rgui. HTH Ruben ________________________________ Ruben H. Roa-Ureta, Ph. D. Senior Scientist, Center for Environment and Water, Marine Studies Section, King Fahd University of Petroleum and Minerals, KFUPM Box 1927, Dhahran 31261,

How about proc <- c(0.9, 0.84, 0.5, 0.16, 0.1) size <- c(0.144, 0.172, 0.272, 0.481, 0.583) plot(size, proc, xlim=c(0,1), ylim=c(0,1)) fit<-nls(proc~plnorm(size, log(xmid), sdlog, lower=FALSE), start=list(xmid=0.2, sdlog=.1)) summary(fit) lines(fitted(fit)~size) -pd > On 10 Jul 2017, at 16:27 , PIKAL Petr <petr.pikal at precheza.cz> wrote: > > Dear all > > I am

Thank you for the report, it is a bug in buffering in R (not specific to Windows) and will be fixed. Best Tomas On 08/17/2017 10:37 AM, PIKAL Petr wrote: > Hi > >> -----Original Message----- >> From: Robert Baer [mailto:rbaer at atsu.edu] >> Sent: Wednesday, August 16, 2017 3:04 PM >> To: PIKAL Petr <petr.pikal at precheza.cz>; Duncan Murdoch >>

Hi try str(varcompobject) to see structure of this object. You can extract parts by standard R means. Cheers Petr > -----Original Message----- > From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of Sharada > Ramadass > Sent: Tuesday, August 8, 2017 3:33 PM > To: r-help at r-project.org > Subject: [R] how to extract individual values from varcomp? > >

Q1. No. Q2. What do you do? Did you follow the instructions at [1]? If you go to the top level of the R help system you should be able to open the manuals. ?help Q3. Not that I know of. If you have difficulty with R on Linux then you should focus on learning how to use Linux in general... and learn what X-Windows is. But just as with Windows, googling your problem will usually lead to stepwise