similar to: How make a x,y dataset from a formula based entry

Mathew Brown Institute of Bioclimatology University of G?ttingen B?sgenweg 2 37077 G?ttingen, Germany t: +49 551 39 9359 mathew.brown at forst.uni-goettingen.de On 9/24/2011 6:00 PM, r-help-request at r-project.org wrote: > Send R-help mailing list submissions to > r-help at r-project.org > > To subscribe or unsubscribe via the World Wide Web, visit >

Hello All, I have a question about glm in R. I would like to fit a model with glm function, I have a vector y (size n) which is my response variable and I have matrix X which is by size (n*f) where f is the number of features or columns. I have about 80 features, and when I fit a model using the following formula,? glmfit = glm(y ~ x1 + x2 + x3 + x4 + x5 + x6 + x7 + x8 + x9 + x10 + x11 + x12 + x13

I don't see the issue here. It would be helpful if people would report their sessionInfo() when reporting whether or not they see this issue. Mine is > sessionInfo() R version 3.4.3 (2017-11-30) Platform: x86_64-pc-linux-gnu (64-bit) Running under: Arch Linux Matrix products: default BLAS/LAPACK: /usr/lib/libopenblas_haswellp-r0.2.20.so locale: [1] LC_CTYPE=en_US.UTF-8

To my surprise, I can confirm on Windows 10 using R 3.4.3 . As tail is not recognized by Windows cmd, I replaced with: system('powershell -nologo "& "Get-Content -Path temp.csv -Tail 1') The last line shows only 7 digits after the decimal, whereas the first have 15 digits after the decimal. I agree with Dirk though, 1.6Gb csv files are not the best way to work with

Dear List, I've around 1000 *.txt files, I've generate with other software. I've now done the following code (below). My question is how can I automate this (with do.call () ?), so it could be done for all the *.txt files. Thanks in advance, Rog??rio names<- list.files() file <- "BLU_Var_%04d.txt" for(i in 1:1000){

Hi - a newbie question, if someone can please help.... I want to change X1, X2,,.....to X.1 X.2 etc in the names below. I am using the Principal Component Regression function (pcr) and it seems to want it this way > datap3.pcr <- pcr(water ~ X, 10, data = datap3, Validation ="cv") Error in model.frame.default(formula = water ~ X, data = datap3) : invalid type (list) for

Hello everyone, I am trying to fit the following model All X. variables are continuous, while the conditions are categoricals. model <- lm(X2

Hello, I have looked on https://www.r-project.org/bugs.html , but it seems that this is the only way to do it. The issue is that the precision used by write.csv is on consistant for big files. See the following code: First I create a large dataframe filled with random uniform values. Then I write it to .csv and print out the first and last lines. df = data.frame(replicate(100, runif(1000000,

Hallo, All, I have a question about changing the height or scale of the y axis. When I use following two R codes, I can get two plots. Please look at the y axes, the number of indices (x1, x2, ?) on the y axis in the first plot is smaller than that in the second plot, and hence the space between any two indices in the first plot is wider than that in the second plot. As the number of indices

Hi All, I have a query about time based sequences. I know such questions have been asked a lot on forums, but I couldnt find the exact thing that I was looking for. I want to create a time-based sequence which will mimic the trading window AND would span multiple days. Something like below: "2011-01-03 09:15:00 IST" "2011-01-03 09:15:01 IST" .... .... .... "2011-01-03

Dear All I would need to perform a MANOVA with both fixed (group, sex, group*sex) and random (brood) effects. I wonder if this is at all possible and if R does that. At the moment, I only know that I can run a classic MANOVA with R. Thank you David ______________________________________________ David Costantini, PhD http://www.davidcostantini.it NERC Postdoctoral research associate

Hi, I am using randomForests for a classification problem. The predict function in the randomForest library, when asked to return the probabilities, has precision of two digits after the decimal. I need at least four digits of precision for the predicted probabilities. How do I achieve this? Thank you, Nagu

Hi, I was wondering how it was possible to access the actual cluster exemplars from the APResult class. Currently it only spits it out onto the terminal if you type the object but there is no other way to see which one is the examplar. Would appreciate any help. Thanks, Sachin [[alternative HTML version deleted]]

If you want a character representation of a long formula (or a formula with long names), you can use: as.character(my.formula) However restriction on length of an as.character result returns only the beginning of a long formula, and without comment. In most cases, the following expression provides the complete result: paste(my.formula[[2]], " ~ ",

I am attempting to use randomForests to look for interesting genes in microarray data with 216genes, 2 classes and 52 samples. My data.frame is 52x217 with the last column, V217 being the class(1 or 2). When I try lung.rf <- randomForest(V217 ~ ., data=tlSA216cda, importance= TRUE, proximity = TRUE) the GUI crashes. I am running R-1.6.2 under windo$e98, and most

Hello everybody, I'm testing the randomForest package in order to do some simulations and I get some trouble with the prediction of new values. The random forest computation is fine but each time I try to predict values with the newly created object, I get an error message. I thought I was because NA values in the dataframe, but I cleaned them and still got the same error. What am I

Hi, I'm afraid that the function read.octave from package "foreign" has some problems with the ASCII data format exported by new versions of Octave (later than 3.2.X). It fails even for a simple case as: [Octave code:] octave:1> x=1; octave:2> save -ascii testdata.mat x [Now in R:] > octavedata <- read.octave('testdata.mat') Mensajes de aviso perdidos In

Hi All, I have data about seed predation (SP) in fruits of three differents colors (yellow, motted, dark) and in two fruiting seasons (2007, 2008). I performed a GLMM (lmer function, lme4 package) and the outcome showed that the interaction term (color:season) was significant, and some combinations of this interaction have significant Pr(>|z|), but I don't think they are the right

Hi, I am trying to run randomForests on a datasets of size 500000X650 and R pops up memory allocation error. Are there any better ways to deal with large datasets in R, for example, Splus had something like bigData library. Thank you, Nagu

Hello all, i'm a newbie of R trying to make some statistical work in R environment. Now i have to laptops, one is Thinkpad X40 with Debian Lenny and the other is Thinkpad T43 with Ubuntu 8.10. Recently i met such problem and am wondering if anybody can do some help. After upgrading my /etc/apt/sources.list , i install R by apt-get install command. It works fine in both laptops. Then i