similar to: fit lognorm to cdf data

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

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,

* fitdistr? * it seems unusual (to me) to fit directly to the data with lognormal... fitting a normal to the log of the data seems more in keeping with the assumptions associated with that distribution. -- Sent from my phone. Please excuse my brevity. On July 10, 2017 7:27:47 AM PDT, PIKAL Petr <petr.pikal at precheza.cz> wrote: >Dear all > >I am struggling to fit data which form

Thank you Martin. If I understand correctly, OP could do wheat.list <- nlsList(Prop ~ SSfpl(end, A, B, xmid, scal), data=wlg) or add some small value to all zeroes wlg$prop < -wlg$Prop+1e-7 wheat.list <- nlsList(prop ~ SSlogis(end,Asym, xmid, scal), data=wlg) which gives fairly reasonable results. plot(augPred(wheat.list)) Am I correct? Cheers Petr > -----Original Message-----

>>>>> PIKAL Petr <petr.pikal at precheza.cz> >>>>> on Fri, 20 Oct 2017 06:33:36 +0000 writes: > Hi > Keep your messages in the list, you increase your chance to get some answer. > I changed your data to groupedData object (see below), but I did not find any problem in it. > plot(wlg) > gives reasonable picture and I am

Hi Duncan The simples spreadsheet is: Put a name in the cell, let say "a1" Put number e.g. 1 below "a1" Copy the number to enough rows Select this column and press ctrl-c result is > temp<- read.delim("clipboard") > str(temp) 'data.frame': 1513 obs. of 1 variable: $ a1: Factor w/ 2 levels "1","a1": 1 1 1 1 1 1 1 1 1 1 ...

Dear all I have 3 questions. Due to some reason I switched from Vista to Ubuntu on home PC. I was used to start with Rgui.exe. However I am not able to find it under Ubuntu and R starts as terminal (probably Rterm). Question 1. Is Rgui.exe available on linux? In Windows doc folder I can find manuals, however I did not find doc folder in Ubuntu. I found somewhere that manuals need to be

Hi Petr; Thanks for your reply. It is much appreciated. A small example is given below for 4 independent and 4 dependent variables only. The values given are regression coefficients.I have looked ggplot documents before writing to you. Unfortunately, I could not figure out as my experience in ggplot is ignorable Regards. Greg y1 y2 y3 y4 x1 -0.19 0.40 -0.06 0.13 x2 0.45 -0.75 -8.67 -0.46 x3

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".

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 Keep your messages in the list, you increase your chance to get some answer. I changed your data to groupedData object (see below), but I did not find any problem in it. plot(wlg) gives reasonable picture and I am not such expert to see any problem with data. Seems to me, that something has to be wrong with nlsList function. > wheat.list <- nlsList(Prop ~ SSlogis(end,Asym, xmid,

Dear all I used to transfer data from excel to R by simple ctrl-c and read.delim("clipboard") construction. I know it is a bad practice but it is easy and for quick exploratory work it is OK. However after changing to new R devel few days ago I encountered weird behaviour. I tried one or two columns. In case of 2 columns, header is repeated after 526 items >

Hi Petr; I would like to get a plot with names as they are in the original file. They are chemical names and I have 733 in the my file. For example, let me give to chemical names "*2-hydroxybutyrate/2-hydroxyisobutyrate*" and "*palmitoyl-arachidonoyl-glycerol (16:0/20:4) [1]**" .So, what should I put [c(2,3,1)] part in the command: iris$MySpecies<-factor(iris$Species,

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,

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,

Hello, I am trying to use varcomp to decompose the variance across multiple nested levels on a lme object. I am able to successfully do this and when I view the varcomp object I can see the individual values / estimates for the variance at different levels. However, I want to be able to extract each of them separately, as I need to build a confidence interval using bootstrapping on the sample

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

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 all; I have 733 discrete categories that will go on y-axis in ggplot2. I used the following command to put the name of x-axis. scale_x_discrete (limits = c("SI", "HOMAIR", "AIR","HOMAB","SG","DI","FI","FG")) Since there are only 8 categories on x it was easy to do. Is there any way to do the same for 733 discrete

Dear friends - I have a data.frame with "," instead of "." and found the discussion http://r.789695.n4.nabble.com/How-to-replace-all-commas-with-semicolon-in-a-string-tt4721187.html#a4721192 so copying the code of Ulrik(I hope:-)) I tried (making some data) AX <-