[R] applying random functions to multisets

2011-03-25 Thread zbynek.jano...@gmail.com
Hi, I am solving following problem: Suppose I have some multiset: > multiset <- c("a","a","c","d","d") and rules, which operate with it (for simplicity not writen in R functions) > rule1: "a" -> c("a","b") > rule2: "a" -> c("a","c") > rule3: "c" -> c("c","c") ... > ruleX: ... I want to apply rule

Re: [R] Incorrectness of mean()

2011-02-28 Thread zbynek.jano...@gmail.com
Oh, I see it now. I guess it confused me, when it did not give me warning and also the numbers were very much alike, so I expected wrong decimal places thanks zbynek -- View this message in context: http://r.789695.n4.nabble.com/Incorrectness-of-mean-tp3327701p3327714.html Sent from the R help m

Re: [R] Substituting inside expression

2011-02-28 Thread zbynek.jano...@gmail.com
Thanks, Your advice solved the problem for one substitution, but I needed to substitute twice: first s = (a+b+c)/2 and than c = sqrt(a^2 + b^2 -2*a*b*cos(gamma)) and I hoped I can do it simultaneously Luckily, I managed to go round the problem and this operation is not necessary anymore:) Thanks a

[R] Incorrectness of mean()

2011-02-28 Thread zbynek.jano...@gmail.com
I have found following problem: I have a vector: > a <- c(1.04,1.04,1.05,1.04,1.04) I want a mean of this vector: > mean(a) [1] 1.042 which is correct, but: > mean(1.04,1.04,1.05,1.04,1.04) [1] 1.04 gives an incorrect value. how is this possible? thanks, zbynek -- View this message in context: h

Re: [R] Substituting inside expression

2011-02-28 Thread zbynek.jano...@gmail.com
Thanks, but it is not exactly what i wanted to achieve. I need to calculate a derivation of this complex formula and I don´t think this would enable me to do so. But luckily, I have manages to built a "by-pass", so this is not necessary anymore zbynek -- View this message in context: http://r.7

[R] Substituting inside expression

2011-02-25 Thread zbynek.jano...@gmail.com
I am having following problem: I´m constructing model for calculation of area of triangle. I know sides a, b, and gamma angle. I wish to calculate the area using heron´s formula: S <- sqrt(s*(s-a)*(s-b)*(s-c)) where s <- (a+b+c)/2 and c is calculated using law of cosines: c <- sqrt(a^2 + b^2 -2*a*

Re: [R] How to optimize function parameters?

2010-12-20 Thread zbynek.jano...@gmail.com
Thanks for your suggestions, I will certainly look at that To answer your question... I can calculate root square error between empirical data a those predicted by model (I used this to optimize the parameters using gafit()). -- View this message in context: http://r.789695.n4.nabble.com/How-t

[R] How to optimize function parameters?

2010-12-20 Thread zbynek.jano...@gmail.com
Hi, I have a dataset and I want to fit a function to it. The function is variogram model (http://en.wikipedia.org/wiki/Variogram) The variogram model is defined by three parameters and I want them to be automatically optimized for real time data. I tried to use gafit {gafit} for this, but there a

Re: [R] How can I draw a line in empirical distribution function?

2010-12-16 Thread zbynek.jano...@gmail.com
I believe that functions yline {fields} xline {fields} are what you´re looking for try to type yline(c(.25,.75)) That shoud do that Best Regards Zbynek Janoska -- View this message in context: http://r.789695.n4.nabble.com/How-can-I-draw-a-line-in-empirical-distribution-function-tp3090675p309071

[R] gafit() - specifying range for results

2010-12-16 Thread zbynek.jano...@gmail.com
I have a problem regarding gafit() from gafit package. Is it possible to specify range of possible results of the function? I am using it for automated fitting if variogram model, which has 3 parameters. None of them can be negative, however some data samples lead to negative values. There are als

Re: [R] How to use lm() output for systemfit() 'Seemingly unrelated regression'

2010-09-06 Thread zbynek.jano...@gmail.com
Thanks, it works fine now. -- View this message in context: http://r.789695.n4.nabble.com/How-to-use-lm-output-for-systemfit-Seemingly-unrelated-regression-tp2525418p2527946.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-pro

[R] How to use lm() output for systemfit() 'Seemingly unrelated regression'

2010-09-03 Thread zbynek.jano...@gmail.com
I am having problem using output of lm() function for further analysing using systemfit package. Basicaly, the problem s following - I generate several formulas using lm() > fo1 <- lm(r98[,2] ~ f98[,1] + f98[,2] + ... + f98[,43]) > fo2 <- lm(r98[,1] ~ f98[,1] + f98[,2] + ... + f98[,43]) and than