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
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
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
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
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
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
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
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
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
, 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
;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
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
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
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.
-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
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
] 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
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
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
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
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
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
.
--
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
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
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
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
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
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
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
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
?
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
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
)*. 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
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
,
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
. 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
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)
.
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
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
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
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
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
}
?
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
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
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
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
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
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
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
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
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
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
-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
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
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
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
] - 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
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
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
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
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
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
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
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
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
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
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:
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
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
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 - 100 of 123 matches
Mail list logo