similar to: List-manipulation

Hi, How do I retrieve the var.test() function code? I had a similar problem once before with another function but getMethods() solved the problem then. Now I tried several combinations for var.test() without success. Regards benjamin -- Benjamin Otto Universitaetsklinikum Eppendorf Hamburg Institut fuer Klinische Chemie Martinistrasse 52 20246 Hamburg [[alternative HTML

Hi, Can I rotate a plot (e.g. a distribution plot) by 90 degrees? The barplot function provides the "horiz" command but that's not availeable for the base package functions. I found an old advice from Paul Murrell on a similar problem suggesting to use viewports (grid package). Yet I couldn't reproduce his examples successfully. And going through the examples in the current grid

Hi, Is there a printed version of S-poetry which can purchased or would I have to print out and bind it myself? Regards, Benjamin -- Benjamin Otto Universitaetsklinikum Eppendorf Hamburg Institut fuer Klinische Chemie Martinistrasse 52 20246 Hamburg [[alternative HTML version deleted]]

Hi, I have the following function: > kurtosis <-function(x) (mean((x-mean(x))^4))/(sd(x)^4) #x is a vector and data > print(mydata) V1 V2 V3 V4 V5 1 1007_s_at DDR1 2865.1 2901.3 1978.3 2 1053_at RFC2 103.6 81.6 108.0 3

Hi, Suppose I have a set of values x and I want to calculate the distribution of the data. Ususally I would use the "density" command. Now, can I use the resulting "density-object" model to generate a number of new values which have the same distribution? Or do I have to use some different function? Regards, Benjamin -- Benjamin Otto Universitaetsklinikum Eppendorf Hamburg

Dear all, Trying to test for heteroskedasticity I tried several test from the car package respectively lmtest. Now that they produce rather different results i am somewhat clueless how to deal with it. Here is what I did: 1. I plotted fitted.values vs residuals and somewhat intuitively believe, it isn't really increasing... 2. further I ran the following tests bptest (studentized

Hi Everyone, I am very new to R. When i run the code yesterday, it was working fine. I was able to find gene annotation. somehow, today, when i try to run it again, it is giving me errors message that i don't understand. It says that it cannot open file. what file is it looking for? ###################################################### ##make sure to install and load annotationtool in R

I want to apply Nyblom-Hansen test with the strucchange package, but I don't know how is the correct way and what is the difference between the following two approaches (leeding to different results): data("longley") # 1. Approach: sctest(Employed ~ Year + GNP.deflator + GNP + Armed.Forces, data = longley, type = "Nyblom-Hansen") #results in: # Score-based CUSUM

Hi, There is a point which is irritating me currently quite a bit and that is an aspect of different behaviour between the %in% function and the smaller/bigger than signs (<>). Here is are two examples to demonstrate what I mean: Example1: > c(1,1,2,2,3,4,4,6,7) %in% c(1,2,3) [1] TRUE TRUE TRUE TRUE TRUE FALSE FALSE FALSE FALSE Right, that is what I expect. Example2: > ps

Hey everybody, i am using the rugarch-package and its great! I have a pretty easy problem, but i just dont get it, so thanks if you can help me. Normally i use: / data(DATANAME) spec = ugarchspec() fit = ugarchfit(data = x[,1], spec = spec) fit slotNames(fit) names(fit at fit) coef(fit) infocriteria(fit) likelihood(fit) nyblom(fit) signbias(fit) head(as.data.frame(fit)) head(sigma(fit))

Hi all, I want to plot the grouped means of some variables. The dependent variables and the grouping factor are stored in different columns. I want to draw a simple line-plot of means, in which the x-axis represents the variables and y-axis represents the means. The means of the groups should be connected by lines. So far, the only function that I could find comes closest to what I'm

Hello, I am interested in fitting a generalized nonlinear regression (gnlr) model with negative binomial errors. I have found Jim Lindsay's package that will do gnlr, but I have having trouble with the particular model I am interested in fitting. It is a threshhold model, where below a certain value of one of the parameters being fitted, the model changes. Here is a sample: Cones:

Hi all, I have Xen 3.0.3 installed on Debian Etch 4.0. The problem I have with it is that vncservers are not run for guest domains (also Etch), nomatter I specify vfb = [ ''type=vnc'' ], or not. Does Xen 3.0.3 support vnc frame buffers, or do I have to upgrade Xen? I know that Xen 3.1 creates vncservers with no problems. Thanks _______________________________________________

hi all i started recently using R and i found myself stuck when i try to analyze microarray data. i use the "affy" package to obtain the intensities of the probes, i have two CTRs and two treated. HG.U133A.Experiment1.CEL HG.U133A.Experiment2.CEL HG.U133A_Control1.CEL HG.U133A_Control2.CEL 1007_s_at 2156.23115 467.75615 364.60615 362.11865

hell all: I have a vector as follows: > head(res) 1007_s_at.value 1053_at.value 117_at.value 121_at.value 1255_g_at.value 0.225801033 0.009747404 0.709517954 0.008825740 0.496859178 1294_at.value 0.005091231 after I convert the res into a data frame I got the following: resx<- data.frame(res) > head(resx) res

Hello! I have compiled Asterisk 1.8.1.1 on my SheevaPlug. It works all right for me, except for one problem that I have encountered: I can only register a SIP client (X-Lite in my case) if the secret field of the extension is left blank. Otherwise it throws a 'bad auth' error. Does anybody have any clue? -- Arik Goldfeld -------------- next part -------------- An HTML attachment was

Hello, I have created a few guest domains, which use virtual network (i.e. their NIC''s sit on 192.168.122.0). It works ok, but I have to figure out how to access these machines with VNC from outside. I guess there should be some way to map dom0''s IP to the internal IP/port. Please advise. Thanks, Arik _______________________________________________ Xen-users mailing list

I'm new to R and some what new to the world of stats. I got frustrated with excel and found R. Enough of that already. I'm trying to test and correct for Heteroskedasticity I have data in a csv file that I load and store in a dataframe. > ds <- read.csv("book2.csv") > df <- data.frame(ds) I then preform a OLS regression: > lmfit <- lm(df$y~df$x) To

Dear mailing list, I have a C function that gives me a wrong result when I run it under Windows 7. This is the code under Linux (RHEL5): > library(phenoTest) > data(epheno) > sign <- sample(featureNames(epheno))[1:20] > score <- getFc(epheno)[,1] > head(score) 1007_s_at 1053_at 117_at 121_at 1255_g_at 1294_at -1.183019 1.113544 1.186186 -1.034779 -1.044456

I am trying to compute the white heteroskedasticity-robust standard errors (also called the Huber standard errors) in a linear model, but I can't seem to find a function to do it. I know that the design library in S+ has something like this (robcov?), but I have not yet seen this library ported to R. Anyone know if there is already a function built into R to do this relatively simple job?