[R] ks.test() with 2 samples vs. 1 sample an distr. function

Dear all, I have a question concerning the ks.test() function. I tryed to calculate the example given on the German wikipedia page. xi <- c(9.41,9.92,11.55,11.6,11.73,12,12.06,13.3) I get the right results when I calculate: ks.test(xi,pnorm,11,1) Now the question: shouldn't I obtain the same or a

[R] use gcheckbox in gWidgets to switch on/off gframe

Dear All, I’m trying to create a GUI using gWidgets. I would like to have a checkbox to “switch on/off” different frames within the GUI. Ideally if one frame is switched on, all other frames would be switched off. Unfortunatly I only came as far as this:  library(gWidgets)  Population <- c("A"

[R] do calculation within expression(), text in plot

Dear Helpers, I would like to add some text to my plot, containing a variable and the calculated value for the variable. As I would like to produce many plots, I hope this can be done automatically. Unfortunately I can't get R to do the calculation for cor(nou,dak) when I

[R] Fit Gumbel Distribution using Method of Moments

Hi all! I'm sure this is a stupid question but I can't find an answer. How can I fit the Gumbel distribution to my data using The Method of Moments in R? Thank you for helping me, Tonja __ R-help@r-project.org mailing list https://stat.eth

Re: [R] own function: computing time

That's perfect, thanks a lot! Tonja Gesendet: Mittwoch, 10. Oktober 2012 um 21:37 Uhr Von: "William Dunlap" An: "tonja.krue...@web.de" , "r-help@r-project.org" Betreff: RE: [R] own function: computing time Your original method would be the following function f <- funct

[R] own function: computing time

Hi all, I wrote a function that actually does what I want it to do, but it tends to be very slow for large amount of data. On my computer it takes 5.37 seconds for 16000 data points and 21.95 seconds for 32000 data points. As my real data consists of 1800 data points it would take ages to u

[R] Remove similar rows from matrix

Hi everybody, I have a matrix (mat) from which I want to remove all rows that differ from other rows in that matrix only by having one ore two NA’s instead of a numbers. I would like to remove rows with more NA’s preferably, so in the end the matrix would look like mat2. Ha

[R] scatterplot3d(); customise axes

Hi all! I’m using scatterplot3d() to show the distribution of data for different locations. As I wound like to show distances between the locations and also label the locations, I was wondering whether there is a function similar to axis() for a 2D plot that works with scatterplot3d

[R] local minima/ maxima

__ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

[R] fitting sine wave

Dear R-helpers   I have 7 data points that I want to fit a continuous curve to, that should look similar to a sine wave My data points would mark the local minima and maxima respectively. This is what I’ve got so far. And I would keep doing so, but sadly nls() then says that it has reached the ma

[R] Accelerating the calculation of the moving average

(y in min(dat$index):max(dat$index)){ m<- mean(dat[(dat$index>=y)&(dat$index<=y+qs+1),3]) erg<-c(erg,m) } It does works, but it takes ages. Is there a faster way to compute the moving average? Thank you, Tonja Krueger __

[R] fit two-parameter lognormal distribution, l-moments

Dear R helpers, I would like to fit a two-parameter lognormal distribution to my data using l-moments. Is there a package that provides this feature? I used the “lmom”-package to fit the three-parameter lognormal distribution to my data as shown beneath. I would like something similar for the

[R] intercept point coordinates

Hi List, Can someone help me to calculate the coordinates of the red and green points? In this example I found their approximate location by trying, but as I have to analyse many similar curves, I’d rather calculate the exact location. data<- c(0.008248005, 0.061242387, 0.099095516, 0.189943027,

[R] Contour line coordinates

Hi all, I used contour() to add contour lines to a plot. Now I’m wondering if there is a way to get an output of the calculated x- and y- coordinates of the contour lines? Tonja __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo

[R] solve integrate(,..) varying limits of integration

Dear List, Is there a way to solve integrate(func.1,x, Inf) $value =0.05 by varying the lower limit of integration (x in the example above)? So far I got: r<- 0.730163 s<--2 func.1<- function(t) {1/(2*pi*sqrt(1-r^2))*exp(-1/(2*(1-r^2))*(s^2-2*r*s*t+t^2))} I can change the lower limit manually,

[R] lmomRFA-package: regsimq()

Hi List! I’m using regsimq() from the “lmomRFA”-package to calculate error bounds for diverse distributions. For example: regsimq(gumfit$qfunc, nrec = lmom.data$n, f = lcdfgum, boundprob = c(0.025, 0.975)) Several times I got this error massage: Fehler in quantile.default(ou, probs = boundpro

[R] Confidence intervals (Weibull, LogNormal, Gumbel)

Is there a way to calculate confidence intervals for other distributions than the GEV distribution (like Weibull, LogNormal, Gumbel). I used the maximum likelihood method to estimate the parameters. (For the GEV distribution I used the extReme package)

[R] QQ-plot – Axes

I would like to change the position of the major tick marks in my qq-plot? Right now the ticks are set at 5.5, 6.0, 6.5 and 7.0. I would like them to be at 5.4, 5.6, 5.8, 6.0, 6.2, 6.4, 6.8 and 7.0. So I would have to remove some of the present ticks. So far I can only add ticks

[R] "lmomRFA" package: error bounds/confidence intervals

Dear List I’m using the ”lmomRFA” package to fit different distributions to my data sample. To calculate the error bounds I used: regsimq(…) and sitequantbounds(…) So my questions are: Are error bounds and confidence intervals the same thing? And: Does regsimq(… boundprob = c(0.05,

Re: [R] Peak Over Threshold values

you in advance, Tonja -Ursprüngliche Nachricht- Von: William Dunlap Gesendet: 27.05.2010 22:13:21 An: "Hutchinson,David [PYR]" ,Tonja Krueger Betreff: RE: [R] Peak Over Threshold values >> -Original Message- >> From: r-help-boun...@r-project.org >&

[R] Peak Over Threshold values

values for psi > oseg<-segmented(olm, seg.Z=~day, > psi=NA,control=seg.control(stop.if.error=FALSE,K=30)) > > plot(oseg,add=TRUE,col=2) > oseg$psi > > > best, > vito > > > Tonja Krueger ha scritto: >> Dear List &

[R] Peak Over Threshold values

Dear List I hope you can help me: I’ve got a dataframe (df) within which I am looking for Peak Over Threshold values as well as the length of the events. An event starts when walevel equals 5.8 and it should end when walevel equals the lower threshold value (5.35). I tried “clu

[R] Peak Over Threshold values/ length of the events

Dear List I hope you can help me: I’ve got a dataframe (df) within which I am looking for Peak Over Threshold values as well as the length of the events. An event starts when walevel equals 5.8 and it should end when walevel equals the lower threshold value (5.35). I tried “clu