similar to: converting "by" to a data.frame?

Dear list members, I am trying to estimate parameters of the AR(1)-GARCH(1,1) model. I have one additional dummy variable for the AR(1) part. First I wanted to do it using garchFit function (everything would be then estimated in one step) however in the fGarch library I didn't find a way to include an additional variable. That would be the formula but, as said, I think it is impossible to add

Dear list members, is there any maximum amount of lines of code that can be copy pasted into the R console? Or some maximum amount of levels of nested parantheses? Whenever I copy paste some big chunk of code into the R console I get errors of missing parantheses etc. even though the parantheses are balanced (i.e. the code runs when run via source or inside some function). Cheers Jannis

Dear list, is it just me or is everybody experiencing this? Since the upgrade to R 2.13.0 on my Mac (OS X 10.6.7, R GUI 1.40 (5751)) the highlighting of matching parantheses doesn't work the way it used to. It doesn't highlight the matching parantheses while typing, but it still does while scrolling over it with the cursor with the arrow keys. This is quite inconvenient since I have to go

Dears, do you know whether it is possible to include any square parantheses (brackets) in an expression to use it as an axis label? e.g. I would like to have a label like (with the sub and superscripts correct): "speed [m s^-1]" I know how to combine an expression with text via paste, but as I run soimething like: a='m*s^{-1}' plot(1:10,main=parse(text=a)) I found now way

Dear all, I just stumbled over the following instruction in the LLVM IR of a C program compiled with clang: %26 = call i32 (...)* bitcast (i32 (i32, i32, i32, i32, i32)* @KeWaitForSingleObject to i32 (...)*)(i32 %23, i32 %24, i32 %25, i32 0, i32 0) Since our LLVM Parser choked on this instruction, I tried to check the documentation, but did not find anything about such nested bitcasts within

I have tried to install the RODBC package with the next command: > install.packages(("ORDBC"), (lib="usr/lib/R/library")) But it is failed. Can I get some information about the installation of this package anywhere? Best regards, Nader -------------------------------------------------------------- Zie ook/see also: http://www.knmi.nl/maildisclaimer.html

Hello to the whole group. I am a newbie to R, but I got my way through and think it is a lot easier to handle than other software packages (far less clicks necessary). However, I have a problem with respect to the summary of regression results. The summary function gives sth like: Residuals: Min 1Q Median 3Q Max -0.46743 -0.09772 0.01810 0.11175 0.42252

Hi, I am trying to convert the following NLMIXED code to NLME, but am running into problems concerning 'Singularity in backsolve'. As I am new to R/S-Plus, I thought I may be missing something in the NLME code. NLMIXED *********** proc nlmixed data=kidney.kidney; parms delta=0.03 gamma=1.1 b1=-0.003 b2=-1.2 b3=0.09 b4=0.35 b5=-1.43 varu=0.5; eta=b1*age+b2*sex+b3*gn+b4*an+b5*pkn+u;

Not paying close attention to detail, I entered the equivalent of pstr<-c("b1=200", "b2=50", "b3=0.3") when what I wanted was pnum<-c(b1=200, b2=50, b3=0.3) There was a list thread in 2010 that shows how to deal with un-named vectors, but the same lapply solution doesn't seem to work here i.e., pnum<-lapply(pstr, as.numeric) or similar vapply

Hey hey, I`m analyzing a data set containing the element contentrations of various samples... I wanted to construct notched boxplots and got quite ugly results for some of the boxplots. The notches are often larger then the hinges which resulted in weird looking edges (even though I`m using a log-boxplot). To avoid this problem I thought about using "normal" log-boxplots and adding some

ruby 1.8.6 (2007-09-24 patchlevel 111) str = ''\&123'' puts "abc".gsub("b", str) puts "abc".gsub("b", "#{str}") puts "abc".gsub("b", str.to_s) puts "abc".gsub("b", ''\&123'') puts "abc".gsub("b", "\&123") >ruby test.rb ab123c

Previously, I've posted queries about this, and thanks to postings and messages in response have recently had some success, to the extent that there is now a package called nlmrt on the R-forge project https://r-forge.r-project.org/R/?group_id=395 for solving nonlinear least squares problems that include small or zero residual problems via a Marquardt method using a call that mirrors the nls()

I seem to recall discussion of an language definition file for S for use with the lgrind utility but I can't find any trace of it in an R Site Search. The lgrind utility takes a file of code in a particular programming language and prepares it for "pretty printing" in LaTeX. In my version the available language definitions are $ lgrind -s When specifying a language case is

Hello I want to make some variables with the ifelse-function, but i don't know how to do it. I want to make these five variables; b2$PRRSvac <- ifelse(b2$status=='A' | b2$status=='Aa',1,0) b2$PRRSdk <- ifelse(b2$status=='B' | b2$status=='Bb',1,0) b2$sanVac <- ifelse(b2$status=='C' | b2$status=='sanAa',1,0) b2$sanDk <-

Dear all, is there anyone who can help me extracting p-values from an Anova object from the car library? I can't seem to locate the p-values using str(result) or str(summary(result)) in the example below > A <- factor( rep(1:2,each=3) ) > B <- factor( rep(1:3,times=2) ) > idata <- data.frame(A,B) > fit <- lm( cbind(a1_b1,a1_b2,a1_b3,a2_b1,a2_b2,a2_b3) ? sex,

Dear All, My doubt about how to integrate a simple kernel density estimation goes on. I have seen the recent posts on integrate density estimation, which seem similar to my question. However, I haven't found a solution. I have made two simple kernel density estimation by: kde.1 <-density(x, bw=sd(x), kernel="gaussian")$y # x<- c(2,3,5,12) kde.2 <-density(y,

Hi I searched in the list and only found questions without answers e.g. http://finzi.psych.upenn.edu/R/Rhelp02a/archive/19955.html : Is there a way to get the same results with lme as with aov with Error()? Can anybody reproduce the following results with lme: id<-c(1,1,1,2,2,2,3,3,3,4,4,4,5,5,5,1,1,1,2,2,2,3,3,3,4,4,4,5,5,5,1,1,1,2,2,2,3,3,3,4,4,4,5,5,5)

Hi folks, Suppose I create the character vector charvec by > charvec<-c("a1.b1","a2.b2") > charvec [1] "a1.b1" "a2.b2" and then I use strsplit on charvec as follows: > splitlist<-strsplit(charvec,split=".",fixed=TRUE) > splitlist [[1]] [1] "a1" "b1" [[2]] [1] "a2" "b2" I was wondering

Dear Wiz[R]ds, I am deeply grateful for the help from Duncan Murdoch, Gray Calhoun, and others. We are almost there. For whatever reason, I can't change the symbol from a circle to a triangle in the upright posture plots. Any ideas? I have included the problem in full. # tritiated (3H)-Norepinephrine(NE) disappearance from plasma # concentrations supine and upright # supine datasu <-

Hello I have a dataset named "b2" with 1521 rows, in that dataset i have 64 rows containing specific information. the rownumbers with specific info are: + i [1] 22 53 104 127 151 196 235 238 249 250 263 335 344 353 362 370 389 422 458 459 473 492 502 530 561 624 647 651 666 671 [31] 715 784 791 807 813 823 830 841 862 865 1036 1051 1062 1068