Re: [R] Should help of estimate in t.test be corrected?

Hi Thanks for your feedback. I didn't think about that. Still, the mean difference is computed for paired, not because there are two samples. IMHO, the help should be updated. Best, Samuel Le 2023-04-03 à 12:10, PIKAL Petr a écrit : Hi You need to use paired option t.test(x=0:4, y

[R] Should help of estimate in t.test be corrected?

estimate mean of x     2 > t.test(0:4, 5:9)$estimate mean of x mean of y     2     7 Best, Samuel __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the po

Re: [R] Why R >= 4.1 warns about installing Rtools on Windows?

R package. From your answers, I may contact the maintainer of devtools. Thanks for reading. Bonne journée, Samuel __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the pos

[R] Why R >= 4.1 warns about installing Rtools on Windows?

Hi, I don't understand why Rtoosl seems to be required even if a package is pure R. For sure, this is just a warning, but it frightens new comers. So, why encouraging people add stuff that will be probably unneeded? Have a nice, Samuel __ R-help@r

[R] [R-pkgs] simplePHENOTYPES v1.3.0

ired mean (intercept) value. Please check it our here: https://cran.r-project.org/package=simplePHENOTYPES Vignettes for the most common scenarios one would want to simulate may be found here: https://cran.r-project.org/web/packages/simplePHENOTYPES/vignettes/simplePHENOTYPES.html Best regar

[R] [R-pkgs] simplePHENOTYPES: SIMulation of Pleiotropic, Linked and Epistatic PHENOTYPES

PHENOTYPES.html Best regards, Samuel Fernandes [[alternative HTML version deleted]] ___ R-packages mailing list r-packa...@r-project.org https://stat.ethz.ch/mailman/listinfo/r-packages __ R-help@r-proj

[R] stats:: spline's method could not be monoH.FC

Hi, I have just noticed that the argument method of the spline function of the stats package does not allow to specify monoH.FC although the documentation tells it should be possible. I know how to program a workaround. This post intends to alert the maintainers. Stay safe, Samuel

Re: [R] Identify does sort the locations

Thanks Paul. I have started using locator(). I also found the gatepoints package that sounds interestering. https://CRAN.R-project.org/package=gatepoints Happy R, Samuel On 15-02-2018 21:50, Paul Murrell wrote: Hi Sorry, I think this has always been the behaviour (and the documentation has

Re: [R] Identify does sort the locations

Hi Paul, Thanks for your answer. I am wondering if in the previous versions, let's say 2.1x, the data were in the selection order... Let me know if there a R fonction taht I can code. Do you think I should use locator function and match points on my own? Best, Samuel

[R] Identify does sort the locations

, 5, 7 as output. set.seed(0); x = rnorm(10); y = rnorm(10); plot(x,y); text(x, y, seq(x)); sel = identify(x,y); sel [1] 1 5 7 Any help appreciated, Samuel R version 3.4.3 (2017-11-30) Platform: x86_64-w64-mingw32/x64 (64-bit) Running under: Windows 7 x64 (build 7601) Service Pack 1 Matrix

[R] Suppress horizontal mean line in beanplot()

;median" . I have tried overallline = F, overallline="n", and overallline="", but without success. Strangely, I could also not find any suggestions on the internet. Does anybody know how to do this? Thanks, Samuel ## Example code using the dataset InsectSprays from d

[R] Unbalanced data in split-plot analysis with aov()

nd, why the structure of the ANOVA changes suddenly. Is there any argument, I could supply to change this behaviour? When I do the same with lm() and subsequent anova(), and calculate F-value for Variety by hand, the estimates are still quite robust. Best regards, Sam -- Samuel Knapp Leh

[R] png: cannot load any more object with static TLS

Hi all, I have an issue running png after loading various packages. After running library(drc) library(RCurl) library(XLConnect) library(ROracle) png("test.png") dev.off() The file fails to write, and I get Warning messages: 1: In png("test.png") : unable to load shared object

[R] Issue writing correct code for Coursera assignment--performing matrix inversion

Hello All, I'm new to this mailing list, so please let me know if I've committed any posting faux-pas. I'm working on an assignment for my Coursera course; please see my code below in which I have tried to write two functions--to perform the task of matrix inversion and then caching that data.

[R] Which is the final model for a Boosted Regression Trees (GBM)?

-way interactions: print(plot(gbm1, i.var=c(1,2), n.trees=1000, return.grid=TRUE)) print(plot(gbm1, 2, return.grid=TRUE)) Hope that helps for others searching for the same thing! Samuel [[alternative HTML version deleted]] __ R-help@r

[R] Bug in rep() function

Did I make any mistake? Or where else should I address this problem? Thanks and best regards, Samuel __ 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

Re: [R] Bug in rep() function

] 1.999555911 [1] 1 [1] /// Why are the first two yielding an integer after multiplying, and the last two don't? Apparently, c(0.8,0.6,0.4,0.2) can't be represented exactly. What would be your approach? Always round numbers first, before giving them to rep() ? Thanks, Samuel

Re: [R] R / JMP interface

https://communities.sas.com/message/199936 see link. Looks like R 3.1 may not yet be supported with JMP 11. _ From: Robert Douglas Kinley Sent: Wednesday, May 07, 2014 12:18 PM To: r-help@r-project.org Cc: Samuel J Gardner Subject: R / JMP interface

[R] lasso cross-validation

Dear all, I am using cv.glmnet in r and I have the following question: The default is 10-fold cross-validation, but it is not clear to me how many times are repeated? Is it 50 repeats? I am sore if me question will be very easy for some people! Many thanks in advance, Samuel

[R] How to catch both warnings and errors?

I tried to use this solution (from over two years ago, but it remains an official demo), but found that it only captured the last warning rather than all of them. Instead, the code below collects all warnings. tryCatch.W.E - function(expr) { W - list() w.handler - function(w){ # warning

[R] R for commercial use

Grüßen i. A. Samuel Zastko ALTE LEIPZIGER Bauspar AG Controlling Tel. (06171) 66-4103, Fax (06171) 66-7500-4103 mailto:samuel.zas...@alte-leipziger.de www.alte-leipziger.dehttp://www.alte-leipziger.de/ ALTE LEIPZIGER Bauspar AG Die Bank für private Baufinanzierung im ALTE LEIPZIGER – HALLESCHE

[R] predict.gam, a question about offset variable and spatial resolution

Dear all, I'm using predict.gam (mgcv package) to predict count data (y) from line transect to a regular grid. My model have this form: y=offset(log(x1*0.6))+s(x2)+s( x3)+s(x4), family=quasipoisson,... the offset is the area covered by a portion of a transect line (length(x1)*observation

[R] Undocumented behavior around daylight savings time?

. -- Samuel Brown [[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 commented

Re: [R] Repeated cross-validation for a lm object

help me once again. Best regards, - Bc.Sc.Agri. Alessandro Samuel-Rosa Postgraduate Program in Soil Science Federal University of Santa Maria Av. Roraima, nº 1000, Bairro Camobi, CEP 97105-970 Santa Maria, Rio Grande do Sul, Brazil -- View this message in context: http://r.789695.n4

Re: [R] Repeated cross-validation for a lm object

Dear Max Thank you for your attention. The train function in the caret package realy does what I need. Best regards, - Bc.Sc.Agri. Alessandro Samuel-Rosa Postgraduate Program in Soil Science Federal University of Santa Maria Av. Roraima, nº 1000, Bairro Camobi, CEP 97105-970 Santa Maria

[R] Repeated cross-validation for a lm object

regards, - Bc.Sc.Agri. Alessandro Samuel-Rosa Postgraduate Program in Soil Science Federal University of Santa Maria Av. Roraima, nº 1000, Bairro Camobi, CEP 97105-970 Santa Maria, Rio Grande do Sul, Brazil -- View this message in context: http://r.789695.n4.nabble.com/Repeated-cross-validation

[R] 3D confidence ellipsoid with ellipse projections onto 2D plane

I have a 3xN matrix of parameters obtained from N regressions where the 3 parameters are jointly statistically significant. I would like to reproduce a 3D confidence ellipsoid projecting 2D ellipses onto the XY plane as in Figure 5.2 in this

[R] R Source Code Request Office For National Statistics UK

are very happy to discuss them. Could you please let me know if ONS would be allowed to view your software's source code? I would appreciate if you could reply at your earliest convenience. Regards, David Samuel ONS Segensworth Road Titchfield PO15 5RR UK (+44) 01329 444663 david.sam

Re: [R] Populate a matrix

mat - matrix(ncol = length(x), nrow = length(y)) for(i in 1:length(x)) { mat[,i] = y} HTH, Samuel -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of fernando.cabr...@nordea.com Sent: 05 October 2011 17:11 To: r-help@r-project.org

Re: [R] check availability of a file in R

Does file.exists answer to your question? file.exists(.RData) If you are not sure of the exact name of the file but know it contains .RData, you can try: List.files(directory,.RData) -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf

[R] glm

? Many thanks, Samuel [[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 commented

Re: [R] Conditional plots in the soiltexture package

Dear Julien It is exactly what I need. Thanks a lot. - Bc.Sc.Agri. Alessandro Samuel-Rosa Postgraduate Program in Soil Science Federal University of Santa Maria Av. Roraima, nº 1000, Bairro Camobi, CEP 97105-970 Santa Maria, Rio Grande do Sul, Brazil -- View this message in context: http

Re: [R] Avoiding Biplot Labels Overllaping

)*. See more detailed description at http://www.statmethods.net/advgraphs/axes.html. Alessandro Samuel-Rosa Post-Graduation Program in Soil Science Federal University of Santa Maria Av. Roraima, nº 1000, CEP 97105-970 Santa Maria, RS, Brazil -- View this message in context: http://r.789695.n4

[R] Conditional plots in the soiltexture package

Dear R users I'd like to know if it is possible to create conditional plots in the soiltexture package. I have data from the particle size distribution of soils derived from two parent materials and I want to use the parent material as conditioner. Thanks a lot. Alessandro Samuel-Rosa Post

[R] convert an xml object into a list on R 2.13

, Samuel [[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 commented, minimal, self

Re: [R] execute r-code stored in a string variable

Yes, you can use: eval(parse(text=c)) On the other hand I would not recommend to use c as a variable name as it is the name of a very important function in the R language to aggregate data. HTH, Samuel -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r

[R] formula used by R to compute the t-values in a linear regression

if someone could provide some clarifications. Samuel [[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

Re: [R] formula used by R to compute the t-values in a linear regression

Exactly. My formula holds only for k=1, this is how I generated it. Do you have any references concerning the rather more careful algorithms? Thanks, Samuel -Original Message- From: peter dalgaard [mailto:pda...@gmail.com] Sent: 01 August 2011 14:45 To: Samuel Le Cc: r-h

Re: [R] formula used by R to compute the t-values in a linear regression

Yes, that's what I was looking for. Many thanks, Samuel -Original Message- From: S Ellison [mailto:s.elli...@lgcgroup.com] Sent: 01 August 2011 15:16 To: Samuel Le; r-h...@stat.math.ethz.ch Subject: RE: formula used by R to compute the t-values in a linear regression -Original

[R] proportion explained by each term in a GAM

Dear list, I have read several posts on this topic. I would use the same methodology as proposed by Simon Wood in this post: http://r.789695.n4.nabble.com/variance-explained-by-each-term-in-a-GAM-td836513.html My first question is: Does anyone know a scientific source (paper, book,...) that

Re: [R] set.seed and for loop

What about: set.seed(1001) total - 0 data - vector(list, 30) for(i in 1:30) { data[[i]] - runif(50) } set.seed(1001) data[[23]] - runif(50) HTH Samuel -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Soyeon Kim Sent: 09 June 2011

[R] residual checking for GAM (mgcv)

Dear list, i'm checking the residuals plots of a gam model after a processus of model selection. I found the best model, all my terms are significant, the r-square and the deviance explained are good, but I have strange residuals plots: http://dl.dropbox.com/u/1169100/gam.check.png

[R] converting a matrix or data.frame into an html table

Dear all, Is there a function in R to convert a matrix or a data.frame into an html table? Many thanks, Samuel [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

[R] help needed with LU decomposition: problem with sparse

matrix? Thanks, Samuel [[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 commented

[R] Model checking for gam (mgcv) result

Dear list, i'm checking the residuals plots of a gam model after a processus of model selection. I found the best model, all my terms are significant, the r-square and the deviance explained are good, but I have strange residuals plots: http://dl.dropbox.com/u/1169100/gam.check.png

[R] Question about GAM (mgcv)

Dear list, i'm using the GAM function from mgcv package. I'm using this syntax: model=gam(y~offset(x)+s(log1p(x1))+s(log1p(x2))+s(x3)+s(x4)+s(5),family=quasipoisson,data=data) and I'm sequentially dropping the single term with the highest non-significant p-value from the model and re-fitting

Re: [R] converting call objects into character

Hi, David and Douglas, thanks for the effort in helping me. It seems that deparse(match.call()) is doing the trick. I learned that the class call is not easy to handle in R. Samuel -Original Message- From: David Winsemius [mailto:dwinsem...@comcast.net] Sent: 03 April 2011 18:56

[R] converting call objects into character

don't manage to convert it into a character to write it into a log file with other informations. Can anyone help? Many thanks, Samuel [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

Re: [R] Graph many points without hiding some

it can do. Thankyou Dennis in particular for translating my code into ggplot, this will be a great help as I get started. Samuel On 1 April 2011 05:07, Greg Snow greg.s...@imail.org wrote: Just a note, Base graphics does support transparency as long as the device plotting to supports

[R] controlling the labels width of a barplot

name,second name,third name,fourth name,fifth name,sixth name,seventh name,eigth name,nineth name)) Does someone know a way to control the size of the font in the barplot function, or to give them an inclination angle? Many thanks, Samuel [[alternative HTML version deleted

[R] Graph many points without hiding some

data as above, using transparent symbols, or something else? Is there some other way of solving this issue that I haven't thought of? Thankyou, Samuel Dennis [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https

[R] passing arguments via ...

I would like to do something like the following: Fancyhist-function(x,...) { # first, process x into xprocess somehow, then ... if (is.null(breaks)) { # yes, I know this is wrong # define the histogram breaks somehow, then call hist: hist(xprocess,breaks=breaks,...) } else { # use

[R] generalized mixed linear models, glmmPQL and GLMER give very different results that both do not fit the data well...

Hi, I have the following type of data: 86 subjects in three independent groups (high power vs low power vs control). Each subject solves 8 reasoning problems of two kinds: conflict problems and noconflict problems. I measure accuracy in solving the reasoning problems. To summarize: binary

[R] Date arithmetic coerces POSIXlt to POSIXct?

Hi. This feels like a bug to me, or at least an undocumented feature, but I thought I'd see what people here thought of it. Consider a POSIXlt object like this one: a - as.POSIXlt (2011-01-23 12:45:45) class (a) [1] POSIXlt POSIXt Fine. Now, if I do some arithmetic on that object, the

Re: [R] text file problem-Help needed

Hello, Make sure you have your file in your workspace directory (you can get it on your R console with the command getwd()). HTH, Samuel -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of George Pantop Sent: 08 February 2011 13:19

Re: [R] Matrix' with Functions

Hello, Here is a quick suggestion: F-function(i,j,a,b,c,d) { res-eval(parse(text=paste(f,i,j,(a,b,c,d),sep=))) return(res) } HTH, Samuel -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Alaios Sent: 03 February

Re: [R] Constrained Regression

should just do that for you. HTH, Samuel -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Jackie Chen Sent: 20 January 2011 16:31 To: R-help@r-project.org Subject: [R] Constrained Regression Hi everyone, I'm trying to perform

[R] UniCox in R

Hello, I am interested in Figure 2 in http://www-stat.stanford.edu/~tibs/ftp/cus.pdf Can anyone tell please how to create this plot? Many thanks Samuel [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] UniCox in R

holtman jholt...@gmail.com wrote: From: jim holtman jholt...@gmail.com Subject: Re: [R] UniCox in R To: Samuel Okoye samu...@yahoo.com Cc: r-help@r-project.org Date: Wednesday, 8 December, 2010, 5:25 Do a little reading on how to use the graphics commands.  I would look at plot lines segments

Re: [R] Integral of PDF

your problem: integrate(function(x) dnorm(x, 500,50), 500-1000, 500+1000) 1 with absolute error 0.74 I centred the interval of integration on 500 instead of on 0 if you integrate between -Inf and Inf. HTH, Samuel -Original Message- From: r-help-boun...@r-project.org [mailto:r-help

[R] t-stat for the coefficients of an ARIMA model

be greatly appreciated. Samuel [[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

[R] t-stat for the coefficients of an ARIMA model

be greatly appreciated. Samuel [[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

[R] adding blocks to arrays

I have a four-dimensional array. I have a three-dimensional array I want to add to it. My particular problem is conceptually very simple because the new array is just going to go into the end of the 4D array; i.e. will be added into the final index. I think I have seen an R command that can do

Re: [R] help cannot put multiple chart Stacked Bar (from PerformanceAnalysis library) in a single plot

I am unfamiliar with the chart.StackedBar function, but since barplot works as expected I thought I'd check that you do know you can obtain stacked bar charts using barplot also? (beside = FALSE). On 10 November 2010 09:41, patrick nguyen patrickqngu...@yahoo.com wrote: Hi I'm having

Re: [R] how do i plot this hist?

To superimpose anything onto a barplot you just need to know the x axis locations of the bars, then everything behaves as expected. You can find this by writing the barplot to an object, ie: bp - barplot(c(1:5),beside=TRUE) bp is now a matrix containing the actual x axis locations of the bars

Re: [R] approxfun-problems (yleft and yright ignored)

. Thanks a million!! Best, Sam On 11 September 2010 15:53, Martin Maechler maech...@stat.math.ethz.ch wrote: MM == Martin Maechler maech...@stat.math.ethz.ch     on Sat, 11 Sep 2010 16:04:37 +0200 writes: SW == Samuel Wuest wue...@tcd.ie     on Thu, 26 Aug 2010 14:34:26 +0100 writes:    SW

[R] adding variable to netCDF file

I would like to open an existing netCDF file and add a variable to it. I am using the ncdf package. This test code gives the idea of what I am trying to do: library(ncdf) print('here we go') print('first, construct netCDF file') t-dim.def.ncdf('t','',1:1,unlim=T,create_dimvar=F)

Re: [R] approxfun-problems (yleft and yright ignored)

. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- project.org] On Behalf Of Samuel Wuest Sent: Wednesday, August 25, 2010 8:20 AM To: r-help@r-project.org Subject: [R

[R] approxfun-problems (yleft and yright ignored)

a namespace (and not attached): [1] affyio_1.16.0 preprocessCore_1.10.0 -- - Samuel Wuest Smurfit Institute of Genetics Trinity College Dublin Dublin 2, Ireland Phone: +353-1-896 2444 Web: http://www.tcd.ie/Genetics/wellmer-2/index.html

[R] AUC

Hello, Is there is any R function computes the AUC for paired data? Many thanks, Samuel [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] contours: lines and colors

Is there a way to get both contour lines and filled colors on the same contour plot? My group here doesn't want me to use the lattice package (long story). I tried over-plotting with contour() on top of filled.contour(). This will work in principle, but the required manipulation of the plotting

[R] What is wrong with this code?

Here it is: read.inputs-function(infile) { for (counter in 1:length(input.list)) { seek(infile,where=0,origin='start') newline-readLines(infile,n=1) while(length(newline)0) { if(!is.na(grep(as.character(input.list[[counter]][1]),newline)[1])) { break }

Re: [R] What is wrong with this code?

? Anyway, I still wouldn't mind some advice on character matching. Thanks. Cable, Samuel B Civ USAF AFMC AFRL/RVBXI wrote: Here it is: read.inputs-function(infile) { for (counter in 1:length(input.list)) { seek(infile,where=0,origin='start') newline-readLines(infile,n=1

Re: [R] What is wrong with this code?

Anyway, I still wouldn't mind some advice on character matching. Thanks. If so we need a reproducible example of what you are doing. OK, let's say I have three strings. Str1=abc. Str2=abcd. Str3=efgh. I want to compare Str1 and Str2 in such a way that R detects that Str2 indeed

Re: [R] What is wrong with this code?

Okay, here is a reproducible example of a possible solution (actually it probably does not do what you want but it's a start for discussion. Note it is actual code and data which is what is needed. === str1 - abc str2 - abcd str3 - efgh

Re: [R] glm

Thank you ver much. Is there is a function in R which is doing penalized cubic regression, say spl.plr(), that if I have weeks = 1:9 I can use somthing like pp - spl.plr(weeks,c(1,3,5,7)) and for 8 and 9 will be linear? Is rcs() library(Design)  doing this? Many thanks, Samuel --- On Tue

[R] glm

in predicting the count variale for the weeks 10, 11 and 12 with treat 10mg and 15mg. Thanks in advance. Regards, Samuel [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r

[R] installing RExcel package

Was wondering if anyone has any experience installing the RExcel package by hand. I think I have all the files needed, but our firewall here prevents RExcelInstaller from going through the internet to get them like it wants to do, and it just gives up. Any ideas? Thanks. --Sam

[R] Power calculation

Hello, Is there any R function which does power calculation for unbalanced groups (n1 neq n2)? Since power.t.test has n Number of observations (per group). Many thanks, Samuel [[alternative HTML version deleted]] __ R-help@r

[R] Problem with library(SSPA)

Hello, I have the fellowing problem and I am thankful for any advice! Regards, Samuel    source(http://bioconductor.org/biocLite.R;) BioC_mirror = http://www.bioconductor.org Change using chooseBioCmirror(). biocLite(SSPA

Re: [R] Problem with library(SSPA)

I can't find his email and I have asked the same question to bioconduc...@stat.math.ethz.ch Regards, Samuel --- On Wed, 9/6/10, Uwe Ligges lig...@statistik.tu-dortmund.de wrote: From: Uwe Ligges lig...@statistik.tu-dortmund.de Subject: Re: [R] Problem with library(SSPA) To: Samuel Okoye samu

Re: [R] Problem with library(SSPA)

Dear Maarten, I have problem to use your package and I would be very thankful if you could help me to solve this. Regards, Samuel --- On Wed, 9/6/10, Uwe Ligges lig...@statistik.tu-dortmund.de wrote: From: Uwe Ligges lig...@statistik.tu-dortmund.de Subject: Re: [R] Problem with library(SSPA

[R] LDA with unstandardized discriminant functions

-rnorm(10,3,2) X-cbind(x1,x2,x3) group-c(a,a,a,b,b,b,b,c,c,c) (xlda1-lda(X,group,scale=TRUE)) (xlda2-lda(X,group,scale=FALSE)) # It gives the same results Thanks for your help, Samuel Royer Tardif Ph.D student Soil ecology Université de Sherbrooke, Qc, Canada

[R] AM/PM strptime %p failing 2.11.0 WinXP

it was a problem with that particular installation of R and Windows. What could I do to get this function working on my Windows XP machine? Thankyou, Samuel Dennis sjdenn...@gmail.com [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

[R] Regression using R

searching on google for such functions but to no avail. Thank you, -- Samuel Bravo Department of Computer Science The University of Texas at El Paso El Paso, TX 79968 [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

[R] FW: fairly simple file I/O

OK, my apologies. I am sure this is a question that has been answered before. But I have looked all over the web and can't find an answer for it. I promise, wasting your time and bandwidth is my last resort. So here goes: I have an ASCII file formatted like so: Label 1.1 Time 1

[R] simulation of censoring data

] - min(tt[i],cen[i])  }  moda - lm(tt~x)  modb - lm(y~x)   Is this simulation correct? Thank you in advance! Samuel [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] arrow plots

Thanks, all, for the help. Much obliged. I realize now that I should have said that I am using lattice graphics. The par() command has not been helpful in convincing lattice to plot outside of the default window. Any other advice is appreciated. Thanks again. [[alternative HTML

[R] arrow plots

Am doing some vector plots with the arrows() function. Works well. But what I need to do is supply an arrow for scaling for the reader. I need to plot an arrow of some known magnitude somewhere on the page (preferably outside the bounds of the plot, so that it can be seen clearly) with some

[R] variation in one variable

Hello, Could you please tell me wether there is any function in R that tell me how many subgroup in one variable I have? So for example if my data are x - c(rnorm(50,50,3),rgamma(50,2,1),runif(50,0,1)) I want to know how many group I have? Many thank in advance, Samuel --- On Thu, 9/17/09

[R] Optimise parameter of a generic function

I would like to know how to determine the best value of a particular parameter in a generic function. My function is: nbin - function(k,Dt) { R - Dt / (k + Dt) q - (k + Dt) / k Pt - c(0:7) for (r in c(0:7)) { Pt[r+1] - (factorial(k+r-1)/(factorial(r)*factorial(k-1)))*((R^r)/(q^k)) } Pt } This

Re: [R] Optimise parameter of a generic function

Typical, spotted a very stupid typo in my code after posting, and when corrected I have solved my own problem. Sorry for bothering everyone. Line below should have read: optimise(f = nbinopt, interval = c(0,100), maximum=FALSE, tol=0.01) Samuel Dennis wrote: optimise(f = nbin063, interval

Re: [R] SVM

TRUE samples in group A and FALSE samples in group B? Best wishes, Samuel PS: I will not give up to understand statistics! --- On Thu, 9/17/09, Steve Lianoglou mailinglist.honey...@gmail.com wrote: From: Steve Lianoglou mailinglist.honey...@gmail.com Subject: Re: [R] SVM To: Samuel Okoye samu

[R] SVM

Hello, I have 12 sample each sample has got 1000 observation, i.e I have a matrix X with 1000 rows and 12 columns! m - svm(t(X)) p - predict (m) Can anyone tell me how to use svmtrain() in R! Many Yhanks, Samuel [[alternative HTML version deleted

Re: [R] SVM

Thank you for your reply! Yes, I am using the svm and I do not have new data (how do you create new data?), all I have these 12 samples which I want to classify (predict) these into two group so that I do have six sample in each group? Many thanks, Samuel --- On Thu, 9/17/09, Steve Lianoglou

[R] setting levels in contourplot()

Does the contourplot() routine have an argument analogous to the levels argument in the contour() routine? More generally, is there a way for the user to fix the contour levels in contourplot()? Thanks. --Sam Cable __ R-help@r-project.org mailing

[R] weird error with contourplot

Sorry if I am cluttering up this list with too many rank beginner questions. But I *am* a rank beginner, and I am at my wit's end with this one: I have a data frame df, with components x,y,z. This command to make a contour plot with lattice graphics works just fine:

[R] weird error with contourplot

a difference, I would be interested to know. Thanks again! --Sam -Original Message- From: Cable, Samuel B Civ USAF AFMC AFRL/RVBXI Sent: Thursday, April 16, 2009 9:55 AM To: 'r-help@r-project.org' Subject: weird error with contourplot Sorry if I am cluttering up this list with too many

[R] (hopefully) simple array op

I have a multidimensional array a, for example, a , , 1 [,1] [,2] [1,]13 [2,]24 , , 2 [,1] [,2] [1,]57 [2,]68 So a is 2x2x2. I have another array b, for example, b [,1] [,2] [1,]9 11 [2,] 10 12 So b is 2x2. I want to tack b onto a

[R] (hopefully) simple array op

Regarding my last mail, I typed too soon. Figure I can do a-array(c(a,b),dim=c(2,2,3)) Seems to work fine. But I would like to hear other ideas if you have any. In particular, I am going to have to do this operation over and over. That last index 3 will have to increase by one every time.

  1   2   >