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

2017-11-15 Thread tonja . krueger
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

2016-01-16 Thread tonja . krueger
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 <-

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

2014-02-28 Thread Tonja Krueger
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

[R] Fit Gumbel Distribution using Method of Moments

2012-11-20 Thread Tonja Krueger
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

Re: [R] own function: computing time

2012-10-11 Thread Tonja Krueger
That's perfect, thanks a lot! Tonja Gesendet: Mittwoch, 10. Oktober 2012 um 21:37 Uhr Von: William Dunlap wdun...@tibco.com An: tonja.krue...@web.de tonja.krue...@web.de, r-help@r-project.org r-help@r-project.org Betreff: RE: [R] own function: computing time Your original

[R] own function: computing time

2012-10-10 Thread tonja . krueger
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

[R] Remove similar rows from matrix

2012-08-22 Thread Tonja Krueger
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.

[R] scatterplot3d(); customise axes

2012-04-13 Thread Tonja Krueger
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

[R] local minima/ maxima

2011-12-04 Thread Tonja Krueger
__ 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

2011-07-27 Thread Tonja Krueger
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

[R] Accelerating the calculation of the moving average

2011-03-22 Thread Tonja Krueger
$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 ___ Handy Internet-Flat ¿ gratis ¿ mit

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

2011-02-14 Thread Tonja Krueger
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

[R] intercept point coordinates

2011-01-17 Thread Tonja Krueger
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

2010-09-15 Thread Tonja Krueger
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

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

2010-09-14 Thread Tonja Krueger
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()

2010-08-16 Thread Tonja Krueger
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 =

[R] QQ-plot – Axes

2010-08-02 Thread Tonja Krueger
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] Confidence intervals (Weibull, LogNormal, Gumbel)

2010-08-02 Thread Tonja Krueger
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] lmomRFA package: error bounds/confidence intervals

2010-07-21 Thread Tonja Krueger
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 =

Re: [R] Peak Over Threshold values

2010-05-31 Thread Tonja Krueger
you in advance, Tonja -Ursprüngliche Nachricht- Von: William Dunlap wdun...@tibco.com Gesendet: 27.05.2010 22:13:21 An: Hutchinson,David [PYR] david.hutchin...@ec.gc.ca,Tonja Krueger tonja.krue...@web.de Betreff: RE: [R] Peak Over Threshold values -Original Message- From: r

[R] Peak Over Threshold values

2010-05-27 Thread Tonja Krueger
)) plot(oseg,add=TRUE,col=2) oseg$psi best, vito Tonja Krueger ha scritto: 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

[R] Peak Over Threshold values

2010-05-26 Thread Tonja Krueger
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

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

2010-05-25 Thread Tonja Krueger
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