search for: roslina

Dear r-users,   I have these data below:   I would like to compare each column with a certain value and count how many in each column less than that specified value.    odd column (1,3,5) will compare with 1.61 and even column (2,4,6) will compare with 75 and I would like to count how many for each column.      I tried these below but it give me just one set of data.   critical  <-

Hi r-users,   I have this data below and would like to plot a barplot overlap with a line. This is my data: > hist_50     pdf_obs pdf_gen.50 1  0.000000   0.000000 2  0.083156   0.125366 3  0.132196   0.158230 4  0.126866   0.149432 5  0.120469   0.127897 6  0.121535   0.104096 7  0.103412   0.082171 8  0.082090   0.063539 9  0.065032   0.048408 10 0.050107   0.036470 11 0.036247   0.027236 12

...0,0, 2.3,0,0,0, 2.1,1.4,0,0,0, 3,2,1,0,0,0... data_1: 1.5, 2.3               ( a single nonzero data between zeros data) data_2: 3.1, 2.5, 2.1,1.4   ( two nonzero data between zeros data) data_3: 3,1,2.5, 3,2,1       ( three nonzero data between zeros data) Thank you so much for any help given. Roslina [[alternative HTML version deleted]]

...: could not find function "tCopula" > norm.cop <- ellipCopula("normal", param = c(0.5, 0.6, 0.7), + dim = 3, dispstr = "un") Error: could not find function "ellipCopula" I'm not sure what is wrong.? Thank you so much for any help given. Regards, Roslina

Hi r-users, I would like to add a plot of vertical line segment with arrow from (77,.6) to (77,0) and also a horizontal  line segment with arrow from (0,0.6) to (77,.6) .  So far this is what I have: plot(sq, cdf, type="l", lwd=4,col="blue",xaxs="i",yaxs="i", xlab= "Rainfall (mm)", ylab= "Random no.", main="Random number  and

Dear r-users,   I have daily rainfall data from 1971 to 2000.  I would like to extract november and december data only.  I would also like to do column bind for november and december, therefore I would like to delete 31 December from december data so that the length of november and december are the same.  Hope somebody can help me.  I tried this below:   > kuantan.dt.1 <-

Hi, I have this code here and try to use sapply code.  But I got error message that I don't really understand to correct. bt   <- c(24.96874, 19.67861, 23.51001, 19.86868); round(bt,2) alp  <- c(2.724234, 3.914649, 3.229146, 3.120719); round(alp,2) bt_alp <- data.frame(bt,alp) sapply(bt_alp, function(bt,alp) ((bt_m/bt)^alp), bt_m = min(bt)) > sapply(bt_alp, function(bt,alp)

Hi, I need some opinion.  I would like to use graph that I generate from R code and save it into word document.  Which format is better? pdf, jpeg or tiff? Thank you. [[alternative HTML version deleted]]

Hi r-users, How do we evaluate the summation of (1/m!) from 0 to infinity (for example). Any help is very much appreciated. Thank you.

Hi R-users, I just change to the new version of R.? I just wonder why everytime?I run my function I will get this message: > source(.trPaths[4], echo=TRUE, max.deparse.length=10000) Thank you.

Hi r-users,   I would like to use Kolmogorov smirnov test but in my observed data(xobs) there are ties.  I got the warning message.  My question is can I do something about it?   ks.test(xobs, xsyn)           Two-sample Kolmogorov-Smirnov test data:  xobs and xsyn D = 0.0502, p-value = 0.924 alternative hypothesis: two-sided Warning message: In ks.test(xobs, xsyn) : cannot compute correct

Hi,   Does R has something similar  to vlookup function in excel?   Thank you for the info. [[alternative HTML version deleted]]

...1 934.9 3     72     1 107.2 4     73     1 236.4 5     74     1  10.5 6     75     1 744.6 7     76     1   9.2 8     77     1 108.7 9     78     1 251.5 10    79     1 197.3 11    80     1 144.1 12    81     1 104.5 13    82     1  17.7 14    83     1 151.8 ... Thank you so much for your help. Roslina [[alternative HTML version deleted]]

Hi, I have a data below and would like to search for positive pairs only and form a new data set. ???????????X1??? ??? ??? ?X2 31.0 9.0 11.0 1.0 1.0 0.0 0.0 0.0 8.0 0.0 0.0 0.0 2.0 2.0 18.0 3.0 0.0 0.0 0.0 0.0 0.0 0.0 10.0 0.0 6.0 0.0 ... The new data will be X1'?? X2' 31.0 9.0 11.0 1.0 2.0 2.0 18.0 3.0 I tried to write the function as: y1y2 <-?

  Hi r-users,   I would like to use gev and my data (annual rainfall ) is as follows:   > head(dat,20) A B C D E F G H I J 1 45.1 41.5 58.5 50.1 46.0 49.1 37.7 49.1 59.8 54.0 2 50.3 39.8 49.4 56.4 49.4 48.8 42.1 49.8 49.4 58.3 3 41.7 39.3 44.6 39.1 35.7 41.5 40.8 40.8 38.5 45.6 4 50.7 33.9 48.4 28.2 35.5 39.1 61.4 17.0 30.7 38.3 5 39.3 30.6 46.9 23.8 25.8

Hi, I'm Roslina, PhD student of University of South Australia, Australia from school Maths and Stats. I use S-Plus before and now has started using R-package. I used to analyse rainfall data using julian date. Is there any similar function that you can suggest to me to be used in R-package? Thank you so much...

...      39.4 at interval beginning 19760917 141200 19780629 080000                38.5 at interval beginning 19780508 195400   I hope somebody can help me with this.  I tried to save in ,csv.excel but quite difficult since it has numbers and words.   Thank you so much for any help given.   Regards,   Roslina [[alternative HTML version deleted]]

Hi R-users, Does R has a topic on newton's method? Thank you for the info.

Hi R-users,   I would like to search for the values of seq that match my rand values.  In excel I will use =VLOOKUP(G2,$E$2:$F$32,2).  For example, for rand=.262 it will give me approximately seq=120 and rand=0.964293344, seq=460 and etc.   E           F     G cdf         seq   rand 0.00E+00    0     0.262123478 1.56E-03    20    0.964293344 1.55E-02    40    0.494827113 5.30E-02    60   

Hi r-users,   I would like to know how to put all the data that is greater than certain value in certain cell for my histogram.  For example, since maximum value of p1 is 588 it doesn't fit in the breaks that we specified.  Can we write breaks >300?   max(p1[,2]) pre.hist <- hist(p1[,2],breaks = c(0,50,100,150,200,250,300),right=FALSE)   > max(p1) [1] 587.2761 > pre.hist <-