Re: [R] R cairo_pdf function does not respect plotting boundaries

do some more testing before committing that change. Paul On 1/03/19 8:13 AM, Lee Steven Kelvin wrote: Hello all, When producing a plot in R using the cairo_pdf device, the resultant plot does not respect the plotting boundaries. Lines and shaded regions will spill over the lower x-axis and the

[R] R cairo_pdf function does not respect plotting boundaries

: https://stackoverflow.com/questions/54892809/r-cairo-pdf-function-does-not-respect-plotting-boundaries Thank you in advance for any insights into this issue. Sincerely, Lee Kelvin -- Dr Lee Kelvin Department of Physics UC Davis One Shields Avenue Davis, CA 95616 USA Ph: +1 (530) 752-1500 Fax: +1 (530)

Re: [R] How to formulate quadratic function with interaction terms for the PLS fitting model?

a = gasTrain, validation = > "LOO") > >> where octane ~ NIR is the model that this example is fitting with. > >> NIR is a collective of variables, i.e. NIR spectra consists of 401 > diffuse > >> reflectance measurements from 900 to 1700 nm. > >> In

[R] Quadratic function with interaction terms for the PLS fitting model?

NIR[0,i]*NIR[1,i] + ... i.e. quadratic with interaction terms. But I don't know how to formulate this. May I have some help please? Thanks, Kelvin [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIB

[R] How to formulate quadratic function with interaction terms for the PLS fitting model?

i.e. quadratic with interaction terms. But I don't know how to formulate this. May I have some help please? Thanks, Kelvin [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see

[R] R2WinBUGS with Multivariate Logistic Regression

Dear R-User, I have written a simple code to analyze some data using Bayesian logistic regression via the R2WinBUGS package. The code when run in WinBUGS stops WinBUGS from running it and using the package returns no results also. I attach herewith, the code and a sample of the dataset. Any s

Re: [R] Clarification on Simulation and Iteration

Saturday, August 1, 2015 3:32 AM, David Winsemius wrote: On Jul 31, 2015, at 6:36 PM, Christopher Kelvin via R-help wrote: > Dear All, > I am performing some simulations for a new model. I run about 10,000 > iterations with a sample of 50 datasets and this returns one set of 50 &g

[R] Clarification on Simulation and Iteration

Dear All, I am performing some simulations for a new model. I run about 10,000 iterations with a sample of 50 datasets and this returns one set of 50 simulated data. Now, what I need to obtain is 10 sets of the 50 simulated data out of the 10,000 iterations and not just only 1 set. The model i

Re: [R] Problem with Newton_Raphson

Thank you very much for everything. Your suggestions were very helpful.  Chris - Original Message - From: Berend Hasselman To: Christopher Kelvin Cc: "r-help@r-project.org" Sent: Thursday, September 20, 2012 10:06 PM Subject: Re: [R] Problem with Newton_Raphson On 20-09-2

Re: [R] Problem with Newton_Raphson

s)*log(scale)+(shape-1)*sum(s)*log(1-((exp(-(scale)*sum(x) -(scale)*sum(x) + (shape)*log((exp(-(scale)*sum(x- (shape)*sum(s)*log((exp(-(scale)*sum(x   return(-log1)   }   start <- c(1,1)   zz<-optim(start,fn=z,data=q,hessian=T)   zz Thank you Chris Kelvin - Original Messa

[R] Problem with Newton_Raphson

)*sum(x) -(p[2])*sum(t) + (p[1])*log((exp(-(p[2])*sum(x- (p[1])*sum(s)*log((exp(-(p[2])*sum(x   return(-log1)   } }   start <- c(1,1)   zz<-optim(start,fn=z,data=q,hessian=T)   zz   m1<-zz$par[2]   p<-zz$par[1]  Thank you Chris Kelvin INSPEM. UPM ___

[R] Optim Problem

Hello, I want to estimate the exponential parameter by using optim with the following input, where t contains 40% of the data and q contains 60% of the data within an interval. In implementing the code command for optim i want it to contain both the t and q data so i can obtain the correct estim

[R] problem with ifelse

  s<-numeric(length(t ((       for(i in 1:length(t)){         for(j in 1:length(v-1))          { ifelse ((t[i]>v[j] & t< v[j+1] ),{z[i]<-v[j];s[i]<-v[j+1]},NA)}}        return(cbind(z,s))}  y(t,v) Chris Kelvin Institute for Mathematical Research UPM

[R] standard error

Dear all,  I want to determine the standard error or the mean squared error for the parameter estimate for beta and eta base on the real data.  Any help on how to obtain these estimated errors. library(survival) d <- data.frame(ob=c(149971, 70808, 133518, 145658, 175701, 50960, 126606, 82329),

[R] How to replace NA with zero (0)

Hello,  When i generate data with the code below there appear NA as part of the generated data, i prefer to have zero (0) instead of NA on my data. Is there a command i can issue to replace the NA with zero (0) even if it is after generating the data?  Thank you library(survival) p1<-0.8;b<-1.5;

[R] standard error

Hello, I have tried obtaining the value of standard error from the code below but i get different values when i compare it with the  standard error obtained from the hessian matrix. Can somebody help me out? Thank you n=100;rr=1000 p1=1.2;b=1.5 sq11=sq21=0 for (i in 1:rr){ t<-rweibull(n,shape=p1,

[R] interval censoring

Hello, May i know whether it is possible to generate data twice from Weibull distribution and use one as the start time and the  other as the end time, below is my code. Any suggestion on how to estimate the parameters of Weibull distribution with interval data will be highly appreciated. Thank y

[R] Standard error

Hello, I have tried obtaining the value of standard error from the code below but i get different values when i compare it with the  standard error obtained from the hessian matrix. Can somebody help me out? Thank you n=100;rr=1000 p1=1.2;b=1.5 sq11=sq21=0 for (i in 1:rr){ t<-rweibull(n,shape=p1,

[R] Interval censorin

Hello, May i know whether it is possible to generate data twice from Weibull distribution and use one as the start time and the  other as the end time, below is my code. Any suggestion on how to estimate the parameters of Weibull distribution with interval data will be highly appreciated. Thank y

[R] R: Help; error in optim

Hello, When i run the code below from Weibull distribution with 30% censoring by using optim i get an error form R, which states that Error in optim(start, fn = z, data = q, hessian = T) :    objective function in optim evaluates to length 25 not 1 can somebody help me remove this error. Is my

[R] R-Help: Censoring data

Hello,  I want to estimate weibull parameters with 30% censored data. I have below the code for the censoring  but how it must be put into the likelihood equation to obtain the desire estimate is where i have a problem with,  can some body help?  My likelihood equation is for a random type-I cens

[R] R-help; generating censored data

Hello,  can i implement this as 10% censored data where t gives me failure and x censored. Thank you p=2;b=120 n=50 set.seed(132); r<-sample(1:50,45) t<-rweibull(r,shape=p,scale=b) t set.seed(123);  cens <- sample(1:50, 5)  x<-runif(cens,shape=p,scale=b)  x Chris Guure Researcher, Institute for

[R] R-help; Censoring

Hello, I wish to censor 10% of my sample units of 50 from a Weibull distribution. Below is the code for it. I will need to know whether what i have done is correct and if not, can i have any suggestion to improve it? Thank you  p=2;b=120 n=50 r=45 t<-rweibull(r,shape=p,scale=b) meantrue<-gamma(

[R] simulation

Hello, i need to simulate 100 times, n=40 ,  the distribution has 90% from X~N(0,1) + 10% from X~N(20,10) Is my loop below correct? Thank you n=40 for(i in 1:100){ x<-rnorm(40,0,1)  # 90% of n z<-rnorm(40,20,10)  # 10% of n } x+z __ R-help@r-project.o

[R] matrix with Loop

Hello! I got something to ask..whether you can help me with the R program...i got this for example 5x4 matrix..and i want to find:  i) mean for each row of the matrix ii) median for each column of the matrix and i need to do this using a loop function...below is my program..u try to check it for

[R] R- Fisher Information

Dear All, Can you help me, with the code below how do I obtain the fisher information from it. Is my q<-replicate(1000,x) the right way to do simulation. thank you. x<-rweibull(100,0.8,1.5) q<-replicate(1000,x) z<-function(p){ beta<-p[1] eta<-p[2] log1<-(n*log(beta)-n*beta*log(eta)+(beta-1)*sum(l

[R] Fisher Imformation

Hello, i have used the code below to estimate the parameters of weibull distribution and i want to obtain the fisher information by providing the the next code but i receive errors anytime i try to, what do i do? by the way is my replication correct and is it placed at the right position for re

[R] r-help; fisher information

Hello, i have used the code below to estimate the parameters of weibull distribution and i want to obtain the fisher information by providing the the next code but i receive errors anytime i try to, what do i do? by the way is my replication correct and is it placed at the right position for rep

[R] r-help; parameter estimate

I need help, the codes below estimates the weibull parameters with complete failure, my question is how do i change the state to include some censoring (may be right, type-I or type-II) to generate and estimate the parameters. thank you x=rweibull(10,2,2) library(survival) d<-data.frame(ob=c(x),

[R] parameter estimate

I need help, the codes below estimates the weibull parameters with complete failure, my question is how do i change the state to include some censoring (may be right, type-I or type-II) to generate and estimate the parameters. thank you x=rweibull(10,2,2) library(survival) d<-data.frame(ob=c(x),

[R] r-help; weibull distribution

 Please, Help me, How do I generate data from the weibull distribution if the data contain both failure and interval censored, For example, I want to generate n=100, shape=2 and scale =4 with 30% interval censored.  What about right censoring Thank you  [[alternative HTML version deleted]

[R] r-help; weibull parameter estimate

Hello, If i write a function as below using log of weibull distribution i do not get the required results in estimating the parameters what do i do, please a/b * (t/b)^a-1 * exp(-t/b)^a n=500 x<-rweibull(n,2,2) z<-function(p) {(-n*log(p[1])+n*log(p[2])- (p[1]-1)*sum(log(x))+(p[1]-1)*log(p[2])

[R] R-help mailing list submissions

R-help mailing list submissions [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commen

[R] R-Help

Can somebody help me, How do I generate data from the weibull distribution if the data contain both failure and interval censored, For example, I want to generate n=100, shape=2 and scale =4 with 30% interval censored. Thank you [[alternative HTML version deleted]] ___

[R] R-HELP

please help; I want to know how to generate an interval-censored data of about 20% and a right censored data of about 30%  using the weibull distribution of say, x=rweibull(100,shape=1.2,scale=1.5) [[alternative HTML version deleted]] __ R-help

[R] Simulation numbers from a probability table

simulated from the probailties table above? A B C D ... 0 4 2 0 2 2 0 1 0 1 4 1 2 2 0 0 ... Thanks for help! Kelvin __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PL

[R] Calculate the percentages of the "numbers" in every column.

0.3 0.1 0.1 1 0.1 0.1 0.2 0.1 2 0.1 0.2 0.2 0.2 3 0.2 0.1 0.1 0 ... Thanks your help! Kelvin __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

[R] Using WCSLIB in R

if anyone in the community has used this library before, and if so how? Thankyou in advance, Lee Kelvin -- School of Physics & Astronomy University of St Andrews North Haugh St Andrews KY16 9SS United Kingdom Phone (+44) [0]1334461668 Email l...@st-andrews.ac.uk

Re: [R] Problems with loading 'wordnet' in a standalone pc

Thanks Uwe. The wordnet is working now. My previous configuration that causes the error was: R - version 2.9.1 wordnet - version 0.1-4 rJava 0.6-2 After updating rJava to version 0.7-0 everything works now. Uwe Ligges-3 wrote: > > Looks like nobody answered so far: > >

[R] Problems with loading 'wordnet' in a standalone pc

ckage/namespace load failed for 'wordnet' My guess is I miss something from the namespace but unfortunately I can't get to the internet. Does anyone know what I'm missing (from the namespace?). Thank you very much! Kelvin -- View this message in context: http://www.nabble.com/

Re: [R] extracting text from wordnet using synonym()

quot; > [6] "assistance" "help" > > On Sat, Aug 15, 2009 at 5:55 PM, Henrique Dallazuanna > wrote: >> Try this: >> >> Str <- c("c(\"aid\", \"assist\", \"assistance\", \"help\")", >>        

[R] extracting text from wordnet using synonym()

] "c(\"aid\", \"assist\", \"assistance\", \"help\")" [2] "c(\"aid\", \"assistance\", \"help\")" How can I get a one character string at the end that looks like this: [1] &qu

Re: [R] re placePatterns() for multiple words replacement in tm

r bad name for your dataset, as it is the > name of a commonly use R function. Ergo, possibility of confusion. > > Bert Gunter > Genentech Nonclinical Biostatisics > > > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]

[R] re placePatterns() for multiple words replacement in tm

Hi all, I wonder how you can replace all words that need to be changed using replacePatterns(). The following is my code. I want to replace both "abc" and "def" to " Yes ". However, I can only replace the first occurrence in sample[[1]]. > sample[[1]] [1] abc def ghi > change <- c("abc","def"

Re: [R] How to Import Excel file into R 2.9.0 version

Use odbcConnectExcel() in the RODBC package. rajclinasia wrote: > > Hi Every one, > I have a problem with Reading Excel file into R 2.9.0 version. In older > versions it is working with "xlsReadWrite" package. But in 2.9.0 version > there is no package like that. so help me out in this aspect

[R] Sorting text docs based on document meta values in tm()

Hi all, I wonder if there's any way to reshuffle the text collection by the document meta values. For instance, if I have 5 documents that correspond to the following meta data: MetaID Sex Age 0 M38 0 M46 0 F 24 0 F 49 0 F 33 Can I re

Re: [R] Entire Organization Switching from SAS to R - Any experience?

We've decided not to proceed with the change but the senior management has been alerted with the uprising SAS license cost. Thank you so much for all the comments!!! Frank E Harrell Jr wrote: > > Jeffrey J. Hallman wrote: >> Hmmm, I sense a story in the offing. Was that an accidental emergen

Re: [R] Entire Organization Switching from SAS to R - Any experience?

I should elaborate the situation a bit more. We store our data in UNIX and have been using UNIX SAS for our work. My Biostat dept has 40 SAS users from which at most 10 also use R. The Epi/Grad Students/Investigators combine for another 30-40 not-so-frequent SAS users let alone R. So we are ta

Re: [R] (newbie) sum for certain number of rows

Thanks. I tried it and it worked wonderful. Wishing for the "DAY' to come. Life needs to be 'reset'. --- On Wed, 7/15/09, Dimitris Rizopoulos wrote: > From: Dimitris Rizopoulos > Subject: Re: [R] (newbie) sum for certain number of rows > To: "kelvin lau&q

[R] (newbie) sum for certain number of rows

on the Internet but so far haven't managed to find it. I am extremely grateful if someone can point me how to go about it. Thanks. Kelvin __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posti

[R] Fractional Polynomials (mfp) for Weibull Model

27;s getting nowhere. I wonder if any of you can give me some tips on how to modify it or any experience modelling other survival models besides Coxph. Thank you very much! Regards Kelvin Lam, MSc. Analyst, Programming & Biostatistics Institute for Clinical Evaluative Sciences (ICES) G1

[R] Levelplot of percentages always using 0 to 100 in the color scheme

I am trying to create levelplot's of cpu usage for systems. print(levelplot(util.mean ~ x.hour * x.day, colorkey=T, cut=20, scales=list(x=list(at=seq(0,96,length=25), labels=ifelse(seq(0,24) %% 4 == 0, seq(0,24), ''))), # add tick marks at the hour main="CPU Utilization During No

[R] Bivariate Survival Analysis

Hi group, Can someone kindly inform which package deals with Bivariate Survival Analysis. I did a search but nothing came up. Thanks! Regards, Kelvin ___ This email may contain confidential and/or privileged

[R] mtry in ctree_control()

Dear Group, What is the actual usage of "mtry" in ctree(), or specifically, ctree_control() since it's a single tree? Thanks in advance. Regards, Kelvin Lam, MSc. Analyst, Programming & Biostatistics Institute for Clinical Evaluative Sciences (ICES) 2075 B