for example:
d = data.frame(gender=c(m,f,unkown), x=rnorm(300))
tapply(d$x,d$gender,mean)
f munkown
0.0787628409 0.0940534765 -0.0005323276
regards,
stefan
On Wed, Dec 09, 2009 at 04:50:37PM -0800, dolar wrote:
Hi there
I have a dataframe of a whole
On Wed, 09-Dec-2009 at 07:20PM -0600, Peng Yu wrote:
| There are a number of functions that are dispatched to from split().
|
| methods('split')
| [1] split.data.frame split.Date split.defaultsplit.POSIXct
|
| Is there a way to figure out which of these variants is actually
|
On Wed, 9 Dec 2009 19:20:47 -0600 Peng Yu pengyu...@gmail.com wrote:
Is there a way to figure out which of these variants is actually
dispatched to when I call split? I know that if the argument is of the
type data.frame, split.data.frame will be called? Is it the case that
if the argument is
Hello,
I came across a problem when building a randomForest model. Maybe someone can
help me.
I have a training- and a testdataset with a discrete response and ten
predictors (numeric and factor variables). The two datasets are similar in
terms of number of predictor, name of variables and
Hi all, I have a script running on a Windows XP 32 bit OS under R 2.9.2
I am experiencing a character string limitation of 64 chars. Is this a
general limitation of R or am I missing something?
Any help appreciated.
Regards,
Ian
--
View this message in context:
Hello
I am searching for a graph for finding the patterns in writing of Authors. I
have around 8 to 10 Authors for whom we have tabulated their writing times.
for e.g the data is as follows:
Author TimeofWriting
A 2009-09-16:57:35:45
A
Hi,
Doing some more reading, I think the problem is easier because the hull is
convex. Then an algorithm for testing points might be:
a) Define the convex hull as a set of planes (simplexes).
[as returned by convhulln!!]
b) Define one point, i, known to be interior
[e.g. mean of all
Hey R people,
I have just put a package on CRAN, mecdf 0.2.1.
It computes multivariate ECDFs.
i.e. Estimates (or perhaps I should say evaluates) a multivariate
cumulative distribution function, using data, without any assumptions
per se.
Plus contrary to my own advise, the vignette contains some
On 09/12/2009 11:29 PM, hadley wickham wrote:
On Wed, Dec 9, 2009 at 9:43 PM, Rolf Turner r.tur...@auckland.ac.nz wrote:
I am working with a somewhat complicated structure in which
I need to deal with a function that takes ``basic'' arguments
and also depends on a number of parameters which
On 10/12/2009 5:05 AM, Ian Platt wrote:
Hi all, I have a script running on a Windows XP 32 bit OS under R 2.9.2
I am experiencing a character string limitation of 64 chars. Is this a
general limitation of R or am I missing something?
You are missing something. There are line length limits on
Hi
r-help-boun...@r-project.org napsal dne 10.12.2009 11:05:15:
Hi all, I have a script running on a Windows XP 32 bit OS under R 2.9.2
I am experiencing a character string limitation of 64 chars. Is this a
general limitation of R or am I missing something?
Here is a snippet from
Hi,
Regarding your proposed algorithm (looks like it is indeed the correct
way to do it), there seem to be a somewhat similar Matlab
implementation,
http://www.mathworks.com/matlabcentral/fileexchange/10226-inhull
It should be possible to port this to R (you might want to check what
to do with
Hi everyone!
Im running the pooltest in plm package, like this
pooltest(cstfin12~lmaccdiscrz+lcobjur+lliq+lcollateral+ldimensao,
data = dados3,
model = within)
But I got the following error:
Error in FUN(X[[1L]], ...) : insufficient number of
observations
My data is an unbalanced panel
Hi,
Currently I am trying to solve a minimization problem using optim as method
Nelder-Mead. However, Neldel-Mead needs many iterations until it finally
converges. I have set $control.trace and $control.report such that I can see
the value of the function at each iteration. I do see that I
Hello Cecilia,
nice hearing from you again. I must restate a couple of my old hints,
though ;^)
1) please always put the authors c/c, as we are not guaranteed to browse
through the r-help every day
2) please provide reproducible examples.
As example(pooltest) keeps working fine, as do some
--
The University of Edinburgh is a charitable body, registered in
Scotland, with registration number SC005336.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
--
The University of Edinburgh is a charitable body, registered in
Scotland, with registration number SC005336.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
--
The University of Edinburgh is a charitable body, registered in
Scotland, with registration number SC005336.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Hi,
I am new to this list and my question may be trivial, I did however not
find an answer to it using the help.search or RSiteSearch features of R.
This is my question:
R seems to break lines of terminal output at a column which is
independent of how wide I have made my terminal. I would
On 10/12/2009 5:15 AM, Keith Jewell wrote:
Hi,
Doing some more reading, I think the problem is easier because the hull is
convex. Then an algorithm for testing points might be:
a) Define the convex hull as a set of planes (simplexes).
[as returned by convhulln!!]
b) Define one point, i,
Try
n=10
b=1
local(
for(i in 1:n) {
n=3
print(n)
b - b*i
}
)
print(n)
print(b)
or
n=10
b=1
local(
for(i in 1:n) {
n=3
print(n)
assign(b, b*i, env=.GlobalEnv)
On Thu, 10 Dec 2009, Michael Olberg wrote:
Hi,
I am new to this list and my question may be trivial,
But the answer depends on the 'at a minimum' information requested in
the posting guide. My comments below only apply to the standard
Unix-alike R and Rterm.exe. Many other front-ends do
Dear R-Community,
I am relatively new with R, so sorry for things which for you might be
obvious...
I am trying to automatically update lmer formulae.
the variables of the model are:
depM= my dependent measure
Sb2= a random factor
OS = a predictor
VR= another predictor
So, I am building
Hi all,
I'm having a problem when putting the four plots (Residuals vs Fitted, Normal
Q-Q, Scale - Location, and Residuals vs Leverage) from an Anova using
plot(aov()) into a pdf. A string aov(df[,i] ~ cat)) is added as the main
title, so when I use mtext() the two titles are overlapping.
Hello,
This may seem like a simple question but I am struggling to answer it.
I wish to discover if there is statistical significance between three
different habitats. I have used a variety of measurements which can be
analysed easily, i.e. GLM for height, distance and circumference of
Hi,everyone:
In MASS chapter 5, the codes that yield a ASH for the duration data is read
below:
breaks - seq(0, 5.9, 0.1)
counts - numeric(length(breaks))
for(i in (0:4)) counts[i+(1:55)] - counts[i+(1:55)] + rep(hist(duration,
breaks=0.1*i + seq(0, 5.5, 0.5), prob=TRUE,
Dear all,
I am trying to plot on a spatial map the punctual measurements of the data
located in the file (https://opengeo.metu.edu.tr/test.csv). I'd like to have
the dots in some way proportional to the magnitudes of the measurements.
I have difficulties with the code given below: The dot sizes
Hi all,
I'm trying to use heat.map to display some data.
The data is originally in csv format and when I read it in R using
read.table I get:
a b c d e
1 A 1 1 1 1
2 B 1 0 1 1
3 C 0 1 0 1
4 D 1 1 1 0
5 E 1 0 0 0
The problem is, that for heatmap I need a numerical matrix. When using
data.matrix
Dear R community,
I would be greatful if somebody could shed light on the following.
I have created a set of 6 points to check how centroid
agglomeration works in cluster analysis:
Y - data.frame(x=c(-1,1,1,-1,10,12),y=c(1,1,-1,-1,0,0))
It is quite intuitive to understand that the last
Hi,
I want to calculate the t- and p-values for a linear model using the Newey
West estimator.
I tried this Code and it usually worked just fine:
oberlm - lm(DYH ~ BIP + Infl + EOil, data=HU_H)
coeftest(oberlm, NeweyWest(oberlm, lag=2))
t test of coefficients:
Hi all,
I want to calculate hazard ratio within each covariate
Example, one covariate has 3 classes (1,2 and 3) and x2 has 2 classes
I want to compare the relative risk ratio within each class of the covariate.
How do I get this result ? .
The other question is that how do I interpret
On Dec 10, 2009, at 8:33 AM, Ashta wrote:
Hi all,
I want to calculate hazard ratio within each covariate
Example, one covariate has 3 classes (1,2 and 3) and x2 has 2 classes
I want to compare the relative risk ratio within each class of the
covariate.
How do I get this result ? .
Moreno,
I leave the discussion on the mixed models to others (you might consider
the SIG group on mixed models as well for this), but try a few hints to
make your code more accessible:
* The . in updating a formula is substituted with the respective old
formula (depending on the side), but is
Hello there Farrel,
I see no one is answering you yet, and I don't have answers to give you, but
I did see this thread:
http://tolstoy.newcastle.edu.au/R/e5/help/08/12/10204.html
http://tolstoy.newcastle.edu.au/R/e5/help/08/12/10204.htmlThat might be
of use to you.
Best,
Tal
David,
Thank you very much for your response.
I fitted the model as factor instead of numeric.
coxfit1 - coxph(Surv(sdat$time, sdat$cens)~factor(y1)+factor(x2)
coef exp(coef) se(coef) z
Pr(|z|)
factor(y1)2 0.036161 1.036822 0.083921 0.431
Hi,
I'm stuck with parsing data into R for heatmap representation.
The data looks like:
1 id1 x1 x2 x3 x20
2 id1 x1 x2 x3 x20
3 id1 x1 x2 x3 x20
4 id1 x1 x2 x3 x20
.
348 id2 x1 x2 x3 x20
349 id2 x1 x2 x3 x20
350 id2 x1 x2 x3 x20
351 id2 x1 x2
Yes, in this case it's comparing to the value '1'. In general with
factors it's comparing to whatever level is missing.
Sam
On Thu, Dec 10, 2009 at 10:40 AM, Ashta sewa...@gmail.com wrote:
David,
Thank you very much for your response.
I fitted the model as factor instead of numeric.
If I want to scale a histogram ie multiply by a constant how would I do
this? Thanks
[[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
Thank you Adrian. Your response was very informative.
?downloadDocument filled me with excitement untill I read, If you try to
download a spreadsheet with multiple worksheets into a 'csv' or 'tsv'
format, only the first worksheet will be downloaded.
So now there is a convenient fast way to read
On Dec 10, 2009, at 8:25 AM, Liat wrote:
Hi all,
I'm trying to use heat.map to display some data.
The data is originally in csv format and when I read it in R using
read.table I get:
a b c d e
1 A 1 1 1 1
2 B 1 0 1 1
3 C 0 1 0 1
4 D 1 1 1 0
5 E 1 0 0 0
The problem is, that for heatmap I need
Dear all,
I am facing problem with inserting the scheduled day of Observation
in the dataset. In the dataset I have only relative time (table 1) and not
scheduled day of observation (day 1, 4, 8, 15, 22, 29, 36, 43).
I would appreciate if any one could suggest me how to proceed.
Eg:
Table
Hi,
As a first exercise with R, I am misusing it purely as a plotting
tool. I have some data that look like this:
date group deptcputime
2008/01 A phys124
2008/01 B phys214
2008/01 C chem281
2008/02 A phys092
2008/02 B phys
Hi,
I want to calculate the t- and p-values for a linear model using the Newey West
estimator.
I tried this Code and it usually worked just fine:
oberlm - lm(DYH ~ BIP + Infl + EOil, data=HU_H)
coeftest(oberlm, NeweyWest(oberlm, lag=2))
t test of coefficients:
Estimate Std.
Farrel,
Please read the manuals. On the RGoogleData package page you can read:
The package provides R access to Google services through the Google
supported Java API.
[...]
A package with very similar functionality is maintained by Duncan Temple
Lang at
On 12/10/2009 04:01 PM, terry johnson wrote:
If I want to scale a histogram ie multiply by a constant how would I do
this? Thanks
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
I have looked at the help for data.matrix. But didn't find a solution.
data.matrix only gets two arguments the data and a logical rownames.force -
I have tried to change this logical to TRUE but that doesn't help.
Any advice would be greatly appreciated.
Thanks!!
Liat wrote:
Hi all,
I'm
Or like this:
d = data.frame(gender=c(m,f,unkown), x=rnorm(300))
ddply(d,gender,summarise,mean=mean(x))
Felipe D. Carrillo
Supervisory Fishery Biologist
Department of the Interior
US Fish Wildlife Service
California, USA
--- On Thu, 12/10/09, smu m...@z107.de wrote:
From: smu
Hello
My dataset is as follows:
jobno recruits
1100 18
1200 1
1850 5
2100 190
2789 25
3000 1
..
..
the dataset has 130 rows
I want to plot barplot from left side and cumulative curve from
On Dec 10, 2009, at 10:19 AM, Liat wrote:
I have looked at the help for data.matrix. But didn't find a solution.
data.matrix only gets two arguments the data and a logical
rownames.force -
I have tried to change this logical to TRUE but that doesn't help.
Any advice would be greatly
I will try to have something in place by Monday to allow you to download
a specific sheet not default to the first. I will let you know.
Adrian
On Thu, 10 Dec 2009, Farrel Buchinsky wrote:
Thank you Adrian. Your response was very informative.
?downloadDocument filled me with excitement
Hello
I am searching for a graph for finding the patterns in writing of Authors. I
have around 8 to 10 Authors for whom we have tabulated their writing times.
for e.g the data is as follows:
Author TimeofWriting
A 2009-09-16:57:35:45
A
On Wed, Dec 9, 2009 at 4:48 PM, David Reiss dre...@systemsbiology.org wrote:
Ideally I would like to be able to use the function f (in my example)
as-is, without having to designate the environment as an argument, or
to otherwise have to use e$x in the function body.
e - new.env()
e$x - 3
f -
Dear Helpful R Users,
I am graphing some data using the beanplot, but I am having trouble
getting the output I desire. I have five tanks (A-E) and 2 groups for
each tank grp1 or grp2, except tank C where there is only grp1. (I only
changed the grouprep to C grp1 for the example) When I plot
?merge
use all=L
On 10 Dec 2009, at 6:06AM, Venkatesh.P wrote:
Dear all,
I am facing problem with inserting the scheduled day of Observation
in the dataset. In the dataset I have only relative time (table 1)
and not
scheduled day of observation (day 1, 4, 8, 15, 22, 29, 36, 43).
I would
Dear R/Statistics-gurus!
I tried to find answer to my hypothetical question and in vain. Sorry, I
don't have a dataset that fits into this hypothetical question and pardon me
if my explanations/use of statistical terms are not accurate.
It does sound a weird question, but I want to rule out that
Y'all,
I would like to have most of the variables in my function to be global
instead of local. I have not found any references that tell me now to
do that. If I have missed a reference please excuse me and tell me
what I have missed.
Thanks,
Keith Jones
Dear R-users,
after several tries with lapply and searching the mailing list, i want
to ask, wheter and how it is possibly to avoid the for-loop in the
following piece of code?
set2-as.data.frame(matrix(rnorm(9),ncol=3))
set2[1,1] - NA
set2[3,2] - NA
set2[2,1] - NA
dimnames(set2)[1] -
I think you are looking for a macro facility. See defmacro in gtools
which is based on Thomas Lumley's function of the same name whose
article you can find in back issues of R News.
On Thu, Dec 10, 2009 at 11:00 AM, Keith Jones keit...@keithljelp.com wrote:
Y'all,
I would like to have most
Hi,
Is the following close enough?
apply(set2, 2, function(x) x[is.na(x)])
HTH,
baptiste
2009/12/10 Andreas Wittmann andreas_wittm...@gmx.de:
Dear R-users,
after several tries with lapply and searching the mailing list, i want to
ask, wheter and how it is possibly to avoid the for-loop in
I'm wondering where I can find the detailed descriptions on R memory
management. Understanding this could help me understand the runtime of
R program. For example, depending on how memory is allocated (either
allocate a chuck of memory that is more than necessary for the current
use, or allocate
Hi babtiste,
thank you very much for your fast answer. your solution is very good,
but i also need the dimnames as in the for-loop for further calculations.
best regards
Andreas
baptiste auguie wrote:
Hi,
Is the following close enough?
apply(set2, 2, function(x) x[is.na(x)])
HTH,
Try this:
lapply(1:ncol(set2), function(idx)set2[is.na(set2[,idx]),idx, drop = FALSE])
On Thu, Dec 10, 2009 at 2:33 PM, Andreas Wittmann
andreas_wittm...@gmx.de wrote:
Hi babtiste,
thank you very much for your fast answer. your solution is very good, but i
also need the dimnames as in the
On Dec 10, 2009, at 10:59 AM, Santosh wrote:
Dear R/Statistics-gurus!
I tried to find answer to my hypothetical question and in vain.
Sorry, I
don't have a dataset that fits into this hypothetical question and
pardon me
if my explanations/use of statistical terms are not accurate.
It
That would be fantastic. How about an option whereby all the worksheets are
downloaded and read into dataframes and appear as a a list of dataframes?
Farrel Buchinsky
Google Voice Tel: (412) 567-7870
On Thu, Dec 10, 2009 at 10:38, Adrian Dragulescu adria...@eskimo.comwrote:
I will try to
Hi all,
I was trying to test the assumption of proportional hazards
assumption, I used the cox.zph function
cox.zph(coxfit6)
Results are:
rhochisqp
x1 -0.03961.397 2.37e-01
x2 0.11079.715 1.83e-03
x3 -0.08857.743
VR's S PROGRAMMING has a lucid (to me) discussion of environments/lexical
scoping starting on p. 63. I am unable to judge whether it is sufficiently
complete to satisfy Rolf's (or others') needs, however.
Bert Gunter
Genentech Nonclinical Biostatistics
-Original Message-
From:
Related...
Rule of thumb:
Pre-allocate your object of the *correct* data type, if you know the
final dimensions.
/Henrik
On Thu, Dec 10, 2009 at 8:26 AM, Peng Yu pengyu...@gmail.com wrote:
I'm wondering where I can find the detailed descriptions on R memory
management. Understanding this
For the case below, you don't need to know anything about how R
manages memory, but you do need to understand basic concepts
algorithmic complexity. You might find The Algorithm Design Manual,
http://www.amazon.com/dp/1848000693, a good start.
Hadley
On Thu, Dec 10, 2009 at 10:26 AM, Peng Yu
I have a situation that I can not predict the final result's dimension.
In C++, I believe that the class valarray could preallocate some
memory than it is actually needed (maybe 2 times more). The runtime
for a C++ equivalent (using append) to the R code would still be C*n,
where C is a constant
On Dec 10, 2009, at 12:11 PM, Ashta wrote:
Hi all,
I was trying to test the assumption of proportional hazards
assumption, I used the cox.zph function
cox.zph(coxfit6)
Results are:
rhochisqp
x1 -0.03961.397 2.37e-01
x2 0.11079.715
HI,
I'm having trouble with a piece of Rscript which keeps outputting
incorrectly. it's something like this: the code reads in from a file which
contains (reformated) input
file-read.table(file=input_file.txt,sep=\t)[,c(1,3:5)]
file.rows- c(nrow(file)/288) # input_file.txt contains 288
Michael,
Thanks a lot for your reply, I have now understood how to fiddle
around with the formulae updates...my question (see my previous e-mail
where I was sketching this problem out) about LME models remains open...
whether:
depM ~ (1 |Sb2) + OS + (1 + OS | Sb2) + VR + (1 + VR | Sb2)
is
Hi All (especially Duncan and Baptiste),
Summary (of lengthy bits below):
I will have a convex hull in multiple (3 to 10) dimensions derived from many
(thousands) of points by geometry::convhulln.
I will need to categorise other 'test' points as inside/outside that convex
hull . e.g. given:
The symbols function may work better than plot for this situation.
--
Gregory (Greg) L. Snow Ph.D.
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
?aggregate
--
Gregory (Greg) L. Snow Ph.D.
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 dolar
Sent: Wednesday, December 09, 2009 5:51 PM
I notice that R basic installation does not take advantage of multi-thread or
multi-core processors. Is there an option for multi-threaded processors? I
currently run all of my analysis on thinkpad laptop and most processing is very
quick, however, when I run stepAIC with interaction terms it
... or perhaps simpler is ?ave
Bert Gunter
Genentech Nonclinical Biostatistics
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Greg Snow
Sent: Thursday, December 10, 2009 9:48 AM
To: dolar; r-help@r-project.org
Subject: Re: [R]
Does the survey package have a function similar to prop.test() Or is there a
way to use svyciprop() to perform a Chi-square test to see if the difference
in proportions is significant?
I'm comparing liberal congregations with conservative congregations in their
sponsorship of HIV/AIDS programs.
Scratch my suggestion. You have a data frame of variables which is what
aggregate() is for.
Bert Gunter
Genentech Nonclinical Biostatistics
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Greg Snow
Sent: Thursday, December 10,
Hello all,
Have some time series data stored in a data.frame, and am plotting it with
ggplot2 (which is totally awesome). I have explored the documentation and
mailing list archives, and I can't see any way to plot a 'smoother' that is
just the K-step moving average.
For example, imagine I had
On Thu, Dec 10, 2009 at 12:34 AM, nivas sbma...@gmail.com wrote:
Hi,
This is the time series data collected from 2001 to 2008 by every
month.so,there are 96 entries.I have done basic statistics.I need to find a
model fitted to forecast this data.This is the mixedpaper collection for
Hi Keith,
A more specific example of what you're looking for might be
helpful--ie do you want to read global variables or set them? You
probably want to look at environments and closures; ?- and ?assign
are good starting points (the latter has an example of Global
Assignment within a function)
Sorry for the previous error.
Dear Helpful R Users,
I am graphing some data using the beanplot, but I am having trouble
getting the output I desire. I have five tanks (A-E) and 2 groups for
each tank grp1 or grp2, except tank C where there is only grp1. (I only
changed the grouprep to C grp1
Renaming the .Rhistory file to .Rosxhistory didn't solve the problem
for me--the file is still empty. And I have confirmed that it is a
problem specific to R.app--running R from Terminal works fine, adding
records to the existing .Rhistory file.
Too bad there is no apparent workaround. Thanks,
On 11/12/2009, at 4:01 AM, terry johnson wrote:
If I want to scale a histogram ie multiply by a constant how would
I do
this? Thanks
This question is pretty vague. Please read the posting guide.
If you want to rescale the y-axis you could do something like:
set.seed(333) # Only half
On 11/12/2009, at 8:30 AM, Maria Gouskova wrote:
Renaming the .Rhistory file to .Rosxhistory didn't solve the problem
for me--the file is still empty. And I have confirmed that it is a
problem specific to R.app--running R from Terminal works fine, adding
records to the existing .Rhistory file.
Thanks to Gabor Grothendieck, Duncan Murdoch, Hadley Wickham, and
Bert Gunter
for their useful input. I'm beginning to get a glimmering of
understanding,
and I think I now have enough to make some progress.
cheers,
Rolf Turner
I have tow vectors one is the subset of another
x is a subset of X Both are vectors with n elements
X[X %in% x] would give me x again rite because it is a subset but i want all
those are not in x from X.
X[which(X != x)] should this do that
Thanks
Ramya
--
View this message in context:
On Thu, Dec 10, 2009 at 3:46 AM, Karl Ove Hufthammer k...@huftis.org wrote:
On Wed, 9 Dec 2009 19:20:47 -0600 Peng Yu pengyu...@gmail.com wrote:
Is there a way to figure out which of these variants is actually
dispatched to when I call split? I know that if the argument is of the
type
You probably want to use a time series package for this. There are
plotting facilities specifically aimed at time series in zoo, xts,
quantmod, timeSeries and latticeExtra. We illustrate with zoo which
has classic graphics and lattice graphics methods:
devAskNewPage(TRUE)
library(zoo)
Given X1,...,Xn ~ t_k(mu,sigma) student t distribution with k degrees
of freedom, mean mu and standard deviation sigma, I want to obtain the
MLEs of the three parameters (mu, sigma and k). When I try traditional
optimization techniques I don't find the MLEs. Usually I just get
k-infty. Does
In findInterval's help page, it says 'v[i[j]] = x[j] v[i[j] + 1]'.
I'm wondering if there is a variant of it such that 'v[i[j]] x[j] =
v[i[j] + 1]'.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
I have a non-linear regression with 8 parameters to solve however it
does not converge ... easily solves the excel ... including the initial
estimates used in the R were found in the excel ... Another question is how
to establish the increments of R by the parameters in the search ..
I am trying to create multiple dataset by group like the following using either
a loop or vectorization:
vocallGp1-subset(vocallsub, Group==1)
vocallGp2-subset(vocallsub, Group==2)
vocallGp3-subset(vocallsub, Group==3)
vocallGp4-subset(vocallsub, Group==4)
vocallGp5-subset(vocallsub, Group==5)
On Dec 10, 2009, at 2:55 PM, Ramya wrote:
I have tow vectors one is the subset of another
x is a subset of X Both are vectors with n elements
X[X %in% x] would give me x again rite because it is a subset but i
want all
those are not in x from X.
X[which(X != x)] should this do that
Happy New Year !
XLSolutions January 2010 R courses schedule is now available online at 9
USA cities for
*** R/Splus Fundamentals and Programming Techniques and R Advanced
Programming***courses in January 2010.
www.xlsolutions-corp.com/courselist.asp
* New York City **
Hi Ramya,
On Dec 10, 2009, at 3:29 PM, David Winsemius wrote:
On Dec 10, 2009, at 2:55 PM, Ramya wrote:
I have tow vectors one is the subset of another
x is a subset of X Both are vectors with n elements
X[X %in% x] would give me x again rite because it is a subset but i want all
Sender: microsoftexchange329e71ec88ae4615bbc36ab6ce411...@magnetar.com
On-Behalf-Of: mailinglist.honey...@gmail.com
Subject: Re: [R] vector help
Message-Id: 118a120f-a549-4bbe-86a4-8157e6529...@journal.report.generator
Recipient: 20...@boyersesf.com
---BeginMessage---
Sender:
Sender: r-help-boun...@r-project.org
On-Behalf-Of: mailinglist.honey...@gmail.com
Subject: Re: [R] vector help
Message-Id: 2154bd29-add9-4cfc-8447-ba5eea3a3...@gmail.com
Recipient: yingmei.la...@magnetar.com
---BeginMessage---
Hi Ramya,
On Dec 10, 2009, at 3:29 PM, David Winsemius wrote:
On
k - infinity gives the normal distribution. You probably don't care
much about the difference between k=1000 and k=10, so you might
try reparametrizing df on [1,infinity) to a parameter on [0,1]...
albyn
On Thu, Dec 10, 2009 at 02:14:26PM -0600, Barbara Gonzalez wrote:
Given X1,...,Xn ~
1 - 100 of 142 matches
Mail list logo