similar to: Help on nested FOR loops

I am trying to simulate a series of ones and zeros (1 or 0) and I am using "rbinom" but realizing that the number of successes expected is not accurate. Any advice out there. This is the example: N<-500 status<-rbinom(N, 1, prob = 0.15) count<-sum(status) 15 percent of 500 should be 75 but what I obtain from the "count" variable is 77 that gives the probability of

Hi friends, I need suggestions/directions on how to producing a waterfall plot for present extend of change in tumour size for a set of respondents in a study.  Example of use of waterfall plot is in the following slides presented at ASCO 2007 by Axel Grothey. Link is

R has wonderful graphics but I am wondering whether there is anything in R to provide customization of tables like PROC Report does in SAS. A Smile costs Nothing      But Rewards Everything Happiness is not perfected until it is shared                                                               -Jane Porter   [[alternative HTML version deleted]]

Hi.. i have an expression of the form: model1<-nls(y~beta1*(x1+(k1*x2)+(k1*k1*x3)+(k2*x4)+(k2*k1*x5)+(k2*k2*x6)+(k3*x7)+(k3*k4*x8)+(k3*k2*x9)+(k3*k3*x10)+ (k4*x11)+(k4*k1*x12)+(k4*k2*x13)+(k4*k3*x14)+(k4*k4*x15)+(k5*x16)+(k5*k1*x17)+(k5*k2*x18)+(k5*k3*x19)+

I am trying to define groupings from levels of factor variables and this the warning message that R give "& not meaningful for factors". The nature of my task is this. I have a variable stage which has the levels (1B, 2A, 2B) - these are the AJCC TNM stages of cancer, and another variable diameter with factor levels ("=< 4", "4 - 6.5, > 6.5; limit values are

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

Dear R-helpers, I am struggling with an optimization problem at the moment and decided to write the list looking for some help. I will use a very small example to explain what I would like to. Thanks in advance for your help. We would like to distribute resources from 4 warehouses to 3 destinations. The costs associated are as follows: Destination >From 1 2 3 Total

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

Full_Name: Derek Elmerick Version: 2.4.0 OS: Windows XP Submission from: (NULL) (38.117.162.243) hello - i have some code that i run regularly using R version 2.3.x . the final step of the code is to build a multinomial logit model. the dataset is large; however, i have not had issues in the past. i just installed the 2.4.0 version of R and now have memory allocation issues. to verify, i ran

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,

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, This is not a problem with R, the calculated results are mathematically correct. This a matrix stability problem. Because of measuring errors, my matrix solution is a bit off. Here is what my equations look like: A11 x11+A12 x12 +A13 x13 = b1 A21 x21+A22 x21 +A23 x23 = b2 A31 x31+A32 x31 +A33 x33 = b3 A is a reading, X is a measured weight, and b is total. The 3 experiments give

Hello Folks, I have a database of 2000+ days with 35 observations each. I am trying to modeling a time series by day, but it seems a problem that I don<t have the time of the observation. I achieve something interesting by using the barplot function, but I`d rather working with ggplot2, since I have the book by Hadley Wickham in hands. I start by transforming my data into a dataframe, but

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

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]], " ~ ",

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]]

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

Dear all, I have encountered a problem when perform stepwise regression. The dataset have more 9 independent variables, but 7 observation. In R, before performing stepwise, a lm object should be given. fm <- lm(y ~ X1 + X2 + X3 + X11 + X22 + X33 + X12 + X13 + X23) However, summary(fm) will give: Residual standard error: NaN on 0 degrees of freedom Multiple R-Squared: 1, Adjusted