search for: zroslina

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

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   

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

Hello! I just uploaded a photo on zroslina's DailyFlog page that I want you to see. Please come and see: http://www.mydailyflog.com/go/invite_register/zroslina/56282492&stc=18 Thanks! Roslina Zakaria ======================================== Got a digital camera? MyDailyFlog is a personal photo-blogging space where you can eas...

Hi R-expert, I try to eliminate 29 Feb but I got an error message below: feb_data1 <- Pooraka_data[Pooraka_data$Month=="2",] feb_28days <- feb_data1 [feb_data1$Day=="28",] feb_29days <- feb_data1 [feb_data1$Day=="29",] ## delete 29 Feb feb_no_29 <- feb_data1 [-(feb_29days),] feb_no_29 <- feb_data1 [-(feb_29days),] Error in xj[i] : invalid subscript

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 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 R-users, Can Maple function be exported to R? I have a jacobian matrix (4X4)?from maple in algebraic form which involve modified Bessel function of the first kind. I just wonder whether we can use algebraic form into R before the value of the parameters can be estimated. Thank you so much for your attention and help.

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

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, I try to use sn package but it give me the following message: > install.packages(repos=NULL,pkgs="c:\\Tinn-R\\sn_0.4-12.zip") Warning: package 'sn' is in use and will not be installed updating HTML package descriptions I did tried a few time to save the .zip file but it give me the same error message. Thank you so much for any help given.

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

Hi r-users,   I try to draw histogram (in terms of probabilty) and superimpose with the gamma pdf. Using this code below I can get the plots BUT the y scale for the density is out of scale.  How do I change the y-axis scale to max 1?  Another thing, how I do draw smooth line rather that points?   Nota that my observed data is hume_pos and the fitted data is rgam1.   hist(hume_pos,prob=TRUE)

Hi R-users, I installed new version of R and Tinn-R and I just wonder why all my latex document change to Tinn-R symbol? Thank you. [[alternative HTML version deleted]]

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

Hi r-users,   Is there any quick way to write this label?   c("0","50","100","150","200","250","300","350","400","450", "500","550","600","650","700","750","800","850","900")   Thank you.

Hi r-users,   I have this code here, but I just wonder how do I use 'sapply' to make it more efficient lamda_cor <- eigen(winter_cor)$values   > lamda_cor [1] 1.3459066 1.0368399 0.8958128 0.7214407   lamda_cxn <- function(dt) { n       <- length(dt)   term    <- vector(length=n, mode="numeric")     for (i in 1:n)   { term[i] <- (dt[i]/n)*log(dt[i]/n) }    

Hi, I have the binned data (observed and generated from model)  that I would like to test using the anderson-darling goodness of fit test.  But I'm not sure which package in R to use. I tried ad.test(...) but it does not recognise the test by Vito Ricci in FITTING DISTRIBUTIONS WITH R   > ad.test(hist_hume_beec[,1],hist_hume_beec[,2]) Error: could not find function "ad.test"

Hi, I would like to overlap the cdf curve for observed and generated data  Here is my code: plot(cdf,main ="CDF of the sum for winter season-Hume",cex.axis=1.2,xlab="Rainfall (mm)", xaxs="i",yaxs="i",col=c("black","red"), lty=c(1,1),ylab="Cumulative probability", xlim=c(0,800),lwd=1) lines(ecdf(datobs))