search for: kyong

...from R is based on number of runs (up or down) not based on a longest length of runs of same events. To be more specific, for example, from a series, HHTHTTTTHHH, the number of runs are 5, and the longest length of runs of the same events is 4. I'll check for the website you mentioned below. Kyong -----Original Message----- From: bogdan romocea [mailto:br44114@gmail.com] Sent: Thursday, December 27, 2007 12:08 PM To: kyong.ho.park@us.army.mil Cc: r-help Subject: RE: [R] A function for random test based on longest run (UNCLASSIFIED) > require(tseries) > ?runs.test Also, take a lo...

...tried read.csv to import an excel file after saving an excel file as csv. But it added alternating rows of fictitious NA values after row number 16. When I applied read.delim, there were trailing several commas at the end of each row after row number 16 instead of NA values. Appreciate your help. Kyong [[alternative HTML version deleted]]

Classification: UNCLASSIFIED Caveats: NONE Dear R users, How can I retrieve an evaluated gradient value from "deriv" function provided below? I want to retrieve b0 value. Appreciate your help. Kyong junk1<-deriv(~(1/b1)*k-b0/b1,"b0",c("k","b0","b1"),formal=T) junk1(k=0,b0=-14.0236,b1=2.44031) [1] 5.746647 attr(, "gradient"): b0 [1,] -0.409784 Classification: UNCLASSIFIED Caveats: NONE [[alternative HTML version deleted]]

...ing “the integral is probably divergent.” integrate(function(y){ + sapply(y,function(y){ + integrate(function(x)5*y*(1/sqrt(2*pi))*exp(-y^2/2)*(1-(1/sqrt(2*pi))*exp(-x^2/2))^4,-Inf,y)$value}) + },-Inf,Inf) I’m not sure whether the syntax is correct or not. Appreciate your help. Kyong [[alternative HTML version deleted]]

Classification: UNCLASSIFIED Caveats: NONE Hello, R users, Has anybody written a function for random test based on the length of longest run of same events. I really appreciate your help. Kyong Park Classification: UNCLASSIFIED Caveats: NONE [[alternative HTML version deleted]]

...question. Sandra suggested using lm instead of lme for a model without a random effect, and Simon suggested anova.lme instead of anova.lm for comparing two models with and without a random effect. For the anova test, both anova.lme and anova work. Thanks again for their quick helpful responses. Kyong -----Original Message----- From: Park, Kyong H Mr ECBC Sent: Friday, March 14, 2008 3:26 PM To: 'r-help@r-project.org' Subject: Lme does not work without a random effect (UNCLASSIFIED) Classification: UNCLASSIFIED Caveats: NONE Dear R users, I'm interested in finding a random eff...

...effect. I'm not sure how to group the data without continuous value which is shown in the error message at the bottom line. If I use 'aov' with Error(Block), is there a test method comparing between with and without the Block random effect. I'm using R 2.4.1. Appreciate your help. Kyong LCU ST1 SURF Block 1 6.71 A N 1 2 6.97 A Y 1 3 6.77 B N 1 4 6.90 B Y 1 5 6.63 C N 1 6 6.94 C Y 1 7 6.79 D N 1 8 6.93 D Y 1 9 6.23 A N 2 10 6.83 A Y 2 11 6.61 B N 2 12 6.86 B Y...

Dear R users, I'm using R 2.1.0 with Windows 2000. In the middle of a R session, an error message pops up saying Rgui.exe generated errors and shuts down the R session. How can I correct this problem? Appreciate your help. Kyong [[alternative HTML version deleted]]

Hello, R users, I applied segmented regression method contributed by Muggeo and got different slope estimates depending on the initial break points. The results are listed below and I'd like to know what is a reasonable approach handling this kinds of problem. I think applying various initial break points is certainly not a efficient approach. Is there any other methods to deal with segmented

Hi R bayesians, I need an advise how to resolve the two different estimates applying a traditional glm (TG) and a bayes glm (BG), and different results depending on the data formats of response data and the prior specs using bayesglm in R. I'm not familiar with bayes estimate and my colleague asked me to look into this because the EPA from France reported a quite different estimates for

...Ramasamy, Vito Ricci. I tried Marc's simple suggestion which was to copy only the data area and paste to a new excel sheet, and that solved my problem. I'll try other suggestions while I'm learning more about R. Unedited answers are provided below followed by my question. Thanks again. Kyong ---------------------------------------------------------------------- From Richard Müller: I import my OpenOffice calc files as follows (OOo or Excel won't make any difference, the csv-format is the same): inp <- (scan(file, sep=";", dec=",", list(0,0), skip = 13, nl...

I appreciate to Roger Koenker, Achim Zeileis and Vito Muggeo for their informative answers. Listed below is unedited replies I got followed by the question I posted. Kyong 1. Roger Koenker: You might try rqss() in the quantreg package. It gives piecewise linear fits for a nonparametric form of median regression using total variation of the derivative of the fitted function as a penalty term. A tuning parameter (lambda) controls the number of distinct segments....