similar to: activate console

Hello, I'd like to understand 'what' is predicting the response for library(mgcv) gam? For example: library(mgcv) fit <- gam(y~s(x),data=as.data.frame(l_yx),family=binomial) xx <- seq(min(l_yx[,2]),max(l_yx[,2]),len=101) plot(xx,predict(fit,data.frame(x=xx),type="response"),type="l") I want to see the generalized function(s) used to predict the response

I've a data frame with this structure: 'data.frame': 1987 obs. of 11 variables: $ site : Factor w/ 24 levels "B(W)","BC-1",..: 1 1 2 2 2 1 1 1 ... $ sampdate : Date, format: "2000-07-18" "2000-07-18" ... $ tclass : Factor w/ 8 levels "Annelida","Arachnida",..: 1 5 5 5 5 ... $ torder : Factor

Ein eingebundener Text mit undefiniertem Zeichensatz wurde abgetrennt. Name: nicht verf?gbar URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20130513/fe848ce7/attachment.pl>

Dear All, would you please provide your thoughts on the following: let us say I have: a <-c(1,5,8,15,32,69) b <-c(8.5,33) and I would like to extract from "a" the two values that are closest to the values in "b", where the length of this vectors may change but b will allways be shorter than "a". So at the end based on this example I should have the result

Hello, I am brand new in this site as well as in R. Please accept my early apology if it's not the right place to ask such a question. I have the following equation: b=[2(1-a)(1+2a)^0.5]/(1+3a) How can I solve the above equation for 'a' using R? Thanks in advance. KD [[alternative HTML version deleted]]

Dear all, I am a R beginner, and I am looking for a way to do the same thing for all levels of a column in a table. Basically, I have a bunch of protein sequences composed of different amino acid residues, and each residue is represented by an uppercase letter. I want to calculate the ratio of different amino acid residues at each position of the proteins. Here is an example table: Proteins

hi, I use R - 2.15(32bit), and want to make a code to clear a console. Actually, I used to run following code to do that but after update the version of R from 2.14 to 2.15, it doesn't work. cls <- function (t) { require(RDCOMClient) wsh <- COMCreate("Wscript.Shell") wsh$SendKeys("\f") invisible(wsh) } cls() or cls <-

I have a data in the following form : CIN TRN_TYP 9079954 1 9079954 2 9079954 3 9079954 4 9079954 5 9079954 4 9079954 5 9079954 6 9079954 7 9079954 8 9079954 9 9079954 9 . . . . . . there are 100 types of CIN (9079954,12441087,15246633,...) and respective TRN_TYP first of all, I want this data to be

suppose I have two sets of vectors: x1,x2,...,xN and y1,y2,...,yN. I want N correlations: cor(x1,y1), cor(x2,y2), ..., cor(xN,yN). my sets of vectors are arranged as data frames x & y (vector=column): x <- data.frame(a=rnorm(10),b=rnorm(10),c=rnorm(10)) y <- data.frame(d=rnorm(10),e=rnorm(10),f=rnorm(10)) cor(x,y) returns a _matrix_ of all pairwise correlations: cor(x,y)

Dear all, knitr writes a comment sign in front of each line of console output. This extremely useful, especially for beginners who are starting to write scripts. One could just compose their script in the console (getting immediate output) and then copy the whole chunk into their .r-file. I have been able to configure the options such that the prompt and the continuing line labels are replaced

Hello all, I cannot seem to figure out this seemingly simple procedure. I want to sort a data frame by a specified character vector. So for : df.. <- data.frame(Season=rep(c("Summer","Fall","Winter","Spring"),4),Obs= runif(length(rep(c("Summer","Fall","Winter","Spring"),4)))) I want to sort the data frame

I have a stupid problem that is currently driving me crazy... Let us suppose that I want to draw a big red square in the middle of my current device (say X11) I tried the following code : pushViewport(viewport(xscale=c(0,1), yscale=c(0,1)), just=c("center", "center")) vp1 <- viewport(x=unit(0.5, "native"), y=unit(0.5, "native"), width=unit(0.4,

Hello, Does the (core) designation in the CRAN Task View indicate packages that provide basic/recommended functionality in that domain? I don't see it documented anywhere on the site. Thanks for your help.

Following the beginning tutorial, I typed help.start() and was given a choice of browsers. I picked chrome but got back that chrome is not found. I cannot seem to change or get rid of it now. I looked at /etc/R/Renviron but don't see anywhere to set browser (R_BROWSER = ${R_BROWSER ...etc) - nothing about chrome. I tried starting as "R_BROWSER= R" - same thing. How do I change the

Hi all, I am looking for a command for doing 2D filtering (rectangular or Gaussian) in R... I have looked at ksmooth, filter and convolve but they seem to be 1D... Any thoughts? Thanks a lot! [[alternative HTML version deleted]]

Hi, is there a way to replicate a data.frame like you can replicate the entries of a vector (with the repeat-function)? I want to do this: x <- data.frame(x, x) (where x is a data.frame). but n times. And it should be as cpu / memory efficient as possible, since n is pretty big in my case. thanks for any suggestions!

Hello :) I'm just starting out with R and would appreciate your help with a couple of problems I am running into. I have used Sys.glob to get a list of all filenames having a particular file extension (in my case, *.txt) I would now like to use this list in the following manner: I would like to use each filename from the list and open that file into a tab separated matrix and proceed. How can

Hi everybody, I would be very grateful if you could give me your thoughts on the following issue. I need to perform Box-Cox (bcPower€) transformation on my data. To do this, I calculated lambda using the function '€powerTransform'€. powerTransform(data) However, I got an error message when performing this function: Convergence failure: return code = 52 I was told by John Fox

1. Minimum Mean Square Estimation (for Multi-Dimensional Linear Regression) 2. Two-Sample Kolmogorov-Smirnov Tests 3. Empirical ANOVA Methods All this to be done on empirical data. I am new to the R world, any help would be appreciated. Thank You, Bill Gilliland Lead Sensor Systems Engineer E303, Tactical Communications, Transport, and Sensors Group: Tactical

Hi! Yesterday I accidentally discovered this: > a <- LETTERS[1:5] > a [1] "A" "B" "C" "D" "E" > > a[1] <- factor(a[1]) > a [1] "1" "B" "C" "D" "E" BUT: > b <- factor(LETTERS[1:5]) > b [1] A B C D E Levels: A B C D E > b[1] <- factor(b[1]) > b [1] A B C D E