On 11-02-12 05:08, beky wrote:
There is a code from SPSS Syntax
do if sub(ATVK,2,2)=01.
comp strata=1.
else if sub(ATVK,2,2)=05 and sub(ATVK,2,2)=27.
comp strata=3.
else if sub(ATVK,4,2)20 or sub(ATVK,6,2)20.
comp strata=4.
Note that your first statement is overwritten by the one below.
sammyny wrote:
I did a fourier transform on a function in time domain to get the following
functions in frequency domain (in latex):
$Y_1[\omega] = \frac{1}{1-\phi_1 e^{-jw}}$
$Y_2[\omega] = \frac{1}{1-(\phi_1 + \phi_2)e^{-jw} +\phi_1\phi_2e^{-2jw}}$
How do I find the spectrum of this
Dear Felipe,
Thank you very much, this is exactly what I need (did find it in the
ggplot2 book).
Best,
Shige
On Sat, Feb 12, 2011 at 10:19 PM, Felipe Carrillo
mazatlanmex...@yahoo.com wrote:
you can just probably add
+ labs(x=Time Elapsed,y=Predicted Probability)
Felipe D. Carrillo
From: greg.s...@imail.org
To: ment...@gmx.net; r-help@r-project.org
Date: Sat, 12 Feb 2011 18:04:34 -0700
Subject: Re: [R] Test for equivalence
Does it make sense for you to combine the 2 data sets and do a 2-way anova
with treatment vs. control as one factor and experiment number as the
Hello,
I have some data file, say, mydata
1,2,3,4,5,6,7
3,3,4,4,w,w,1
w,3,6,5,7,8,9
4,4,w,5,3,3,0
i want to replace some percentages of mydata file in to NAs for those values
that are NOT w's. I know how to apply the percentage thing here but don't know
how to select those values that are not
Daniel,
I think you might benefit from (re)reading the Intro to R materials, since you
seem to misunderstand both ifelse() and possibly how subsetting works.
I can't provide an actual example since I have no idea what sample or var
or sm are (and it's a VERY good idea to not use names for
Dear All,
I am trying to recode a variable using the functions provided by
memisc package. Actually I am following the examples on page 9-10 of
the vignette:
--
d.fig - within(d.fig,{
sev - recode(sev,
1 - 0.9,
Daniel:
1. Please read up on how to create reproducible examples to post in R,
e.g. via dput. What you inluded below is useless to reproduce your
example.
2. If I understand you correctly, the trick here is to get the
indices, e.g. with ?which:
# First create an example:
mydata -
Hi there, I have a numeric vector let say:
Vect - c(12.234, 234.5675, 1.5)
Now I want a string vector like:
changedVec - c(012.234, 234.568, 001.500)
Would be grateful if somebody help me how can I do that.
Thanks and regards,
[[alternative HTML version deleted]]
On Feb 13, 2011, at 10:43 AM, Bogaso Christofer wrote:
Hi there, I have a numeric vector let say:
Vect - c(12.234, 234.5675, 1.5)
Now I want a string vector like:
changedVec - c(012.234, 234.568, 001.500)
?sprintf
David Winsemius, MD
West Hartford, CT
Bogaso wrote:
Vect - c(12.234, 234.5675, 1.5)
Now I want a string vector like:
changedVec - c(012.234, 234.568, 001.500)
sprintf(%06.3f, c(12.234, 234.5675, 1.5))
Dieter
--
View this message in context:
Is this what you want:
Vect - c(12.234, 234.5675, 1.5)
sprintf(%07.3f, Vect)
[1] 012.234 234.567 001.500
On Sun, Feb 13, 2011 at 10:43 AM, Bogaso Christofer
bogaso.christo...@gmail.com wrote:
Hi there, I have a numeric vector let say:
Vect - c(12.234, 234.5675, 1.5)
Now I want a
On 2011-02-13 07:43, Bogaso Christofer wrote:
Hi there, I have a numeric vector let say:
Vect- c(12.234, 234.5675, 1.5)
Now I want a string vector like:
changedVec- c(012.234, 234.568, 001.500)
Just for completeness, let's add formatC to your options:
formatC( Vect, digits=3,
Shige Song wrote:
Dear All,
I am trying to recode a variable using the functions provided by
memisc package. Actually I am following the examples on page 9-10 of
the vignette:
--
d.fig - within(d.fig,{
sev -
wanfairos wrote:
Dear members;
I have one problem to display my simulated data. Basically, I've generated
a set of rnorm (std.norm.mat)and rnbinom (gen.data) data to test on the
performance of FE Poisson model estimates
... Long self-contained example shorted
wanfairos wrote:
Dear members;
I have one problem to display my simulated data. Basically, I've generated
a set of rnorm (std.norm.mat)and rnbinom (gen.data) data to test on the
performance of FE Poisson model estimates
... Long self-contained example shorted
wanfairos wrote:
I have one problem to display my simulated data. Basically, I've generated
a set of rnorm (std.norm.mat)and rnbinom (gen.data) data to test on the
performance of FE Poisson model estimates
... Long self-contained example shorted
my.data-rep(,sim.size)
for(i in
wanfairos wrote:
I have one problem to display my simulated data. Basically, I've generated
a set of rnorm (std.norm.mat)and rnbinom (gen.data) data to test on the
performance of FE Poisson model estimates
... Long self-contained example shorted
my.data-rep(,sim.size)
for(i in
On 2011-02-13 07:05, Shige Song wrote:
Dear All,
I am trying to recode a variable using the functions provided by
memisc package. Actually I am following the examples on page 9-10 of
the vignette:
--
d.fig- within(d.fig,{
On 11-02-13 10:48, David Winsemius wrote:
On Feb 13, 2011, at 10:43 AM, Bogaso Christofer wrote:
Hi there, I have a numeric vector let say:
Vect - c(12.234, 234.5675, 1.5)
Now I want a string vector like:
changedVec - c(012.234, 234.568, 001.500)
?sprintf
Just for my own
just using as.character does not add the leading/trailing zeros:
Vect - c(12.234, 234.5675, 1.5)
as.character(Vect)
[1] 12.234 234.5675 1.5
On Sun, Feb 13, 2011 at 12:04 PM, Sascha Vieweg saschav...@gmail.com wrote:
On 11-02-13 10:48, David Winsemius wrote:
On Feb 13, 2011, at 10:43 AM,
Hello,
Could somebody please tell me what the following is about?
I try to run BMS function for quite a small dataset (400 rows and 50
columns) to get model suggestions, but I end up getting following
error message:
bms-bms(dat)
Error in chol.default(XtX.start) :
the leading minor of
Hello
I have the following model
y ~ a + b + c + d + e + f
I would like to look at all possible 3 way interaction, is it possible to do
this without typing out each interaction manually? Is there some way to run a
saturated model permitting only 2 and 3 way interactions (As opposed to
How to interprete the results of panel data models of R? I estimate a adapted
form of Koenker's (2004) suggestion for a quantile regression approach with
panel data, for my data:
rq.fit.panel - function(X,Y,s,w,taus,lambda)
{
require(SparseM)
require(quantreg)
K - length(w)
if(K !=
Quantile Regression for Longitudinal Data
--
View this message in context:
http://r.789695.n4.nabble.com/Question-of-Quantile-Regression-for-Longitudinal-Data-tp883458p3303675.html
Sent from the R help mailing list archive at Nabble.com.
__
Please consider following string:
MyString - ABCFR34564IJVEOJC3434
Here you see that, there are 4 groups in above string. 1st and 3rd groups
are for english letters and 2nd and 4th for numeric. Given a string, how can
I separate out those 4 groups?
Thanks for your time
[[alternative
Dear all,
I found a bug in the which() function.
When trying to remove elements with the which function,
if the criteria is not matched, numeric(0) is returned instead of the
array itself.
This is very weird.
a = c(1,2,3,4,5)
a[!a==6]
[1] 1 2 3 4 5
a[-which(a==6)]
numeric(0)
I think your rewrite is overdue, because returning the array from the which
function seems counterintuitive.
---
Jeff Newmiller The . . Go Live...
DCN:jdnew...@dcn.davis.ca.us Basics: ##.#. ##.#. Live Go...
Live:
If by bug you mean function behaving exactly as documented.
which() returns only the matches, the TRUE values. If there are no
matches, it doesn't return anything.
If I understand what you are trying to do, and I may not,
a[which(a != 5)]
is really what you want, and it is precisely to preserve
On Sun, Feb 13, 2011 at 10:27 AM, Megh Dal megh700...@gmail.com wrote:
Please consider following string:
MyString - ABCFR34564IJVEOJC3434
Here you see that, there are 4 groups in above string. 1st and 3rd groups
are for english letters and 2nd and 4th for numeric. Given a string, how can
I
On 2011-02-13 06:54, 7...@queensu.ca wrote:
Hello
I have the following model
y ~ a + b + c + d + e + f
I would like to look at all possible 3 way interaction, is it possible to do
this without typing out each interaction manually? Is there some way to run a
saturated model permitting only 2
If you have an indeterminate number of the patterns in the string, try
the following:
MyString - ABCFR34564IJVEOJC3434
# translate to the pattern sequences
x - chartr('ABCDEFGHIJKLMNOPQRSTUVWXYZ0123456789'
+ , '0011'
+ , MyString
+
Hi, all,
I have a question on the acf function in time sereis analysis. After I run
the function
acf(data);
I got two blue lines above and below the x axis. I guess they are confidence
inverval borders, right?
Could someone tell me what the confidence level is for the blue lines? Is
there some
Thanks Sarah,
Yes, the function behaves Exactly as documented:
check this out:
a = c(1,2,3,4,5)
a[which(a!=6)]
[1] 1 2 3 4 5
a[!which(a==6)]
numeric(0)
a[-which(a==6)]
numeric(0)
a[!a==6]
[1] 1 2 3 4 5
I guess this is just a gotcha, since
I often use !which and -which to remove elements,
Sybil:
Have a look at
?map.grid
for a start.
You will need to fiddle if you really don't want the grid lines.
HTH,
Ray Brownrigg
On Sun, 13 Feb 2011, sybil wrote:
Hi,
I am new to R and I want to display longitude and latitude for a projected
map.
map.axes won't do it since it only works
Hi
Can some one please point out where i am wrong.
I am trying to position set of nodes column-wise in cytoscape using
RCytoscape
AD
BE
CF
---
g - new ('graphNEL', edgemode='undirected')
cw - CytoscapeWindow ('smallExample', graph=RCytoscape::makeSimpleGraph())
Hi all,
I'm using
xyplot(closingDataXTS)
To get a page with any number of seperate plots on it (as many plots as
columns in closingDataXTS).
Each plot is named according to colnames(closingDataXTS).
I would like to control the size of the text each plot name appears in.
I've seen a number of
typing 'markov switching' into the R site search box gives 127 results,
certainly a few of which on the first page are possibly of interest to you.
among them packages MSBVAR and bayesGARCH
hth, best, Ingmar
On Sat, Feb 12, 2011 at 12:53 PM, Nanda Mendez nand...@hotmail.com wrote:
I'm looking
Hello ,
I have a general question:
What does it means that a package was removed from the CRAN repository?
There were problems in the package?
I want to use 'mlice' package, but the package doesn't exist in the CRAN now
and I worry about...
I will be very thankful for the answer!
Pavel
Yeah, it works!
I use lty=0 to suppress the grids.
If I can get the longitude to display on the bottom instead of the top, i
would be happier.
But I am happy with it now.
Thanks,
--
View this message in context:
On Feb 13, 2011, at 2:38 PM, Pavel Goldstein wrote:
Hello ,
I have a general question:
What does it means that a package was removed from the CRAN
repository?
There were problems in the package?
The CRAN Package Checks page may have the answer:
I want to use 'mlice' package, but the
On Feb 13, 2011, at 2:38 PM, Pavel Goldstein wrote:
Hello ,
I have a general question:
What does it means that a package was removed from the CRAN
repository?
There were problems in the package?
It appears it may never have been on CRAN. There are other sources of
packages: BioC,
Dear all, I admit this is not anything to do R and even with
Statistics perhaps. Strictly speaking this is a math related question.
However I have some reasonable feeling that experts here would come up
with some elegant suggestion to my question.
Here my question is: What is sum of all Integers?
-Original Message-
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On Behalf Of Mike Lawrence
Sent: Friday, February 11, 2011 9:28 AM
To: r-h...@lists.r-project.org
Subject: [R] Help optimizing EMD::extrema()
Hi folks,
I'm attempting to use the EMD
I did not try to emulate the ndata nad ndatam1 arguments
to extrema(), as I didn't see what they were for. The
help file simply says they are the length of the first
argument and that minus 1 and that is what their default
values are. If they do not have their default values
then extrema()
On Sun, Feb 13, 2011 at 4:42 PM, Megh Dal megh700...@gmail.com wrote:
Hi Gabor, thanks (and Jim as well) for your suggestion. However this is not
working properly for following string:
MyString - ABCFR34564IJVEOJC3434.36453
strapply(MyString, (\\D+)(\\d+)(\\D+)(\\d+), c)[[1]]
[1] ABCFR
On Feb 13, 2011, at 4:47 PM, Maithula Chandrashekhar wrote:
Dear all, I admit this is not anything to do R and even with
Statistics perhaps. Strictly speaking this is a math related question.
However I have some reasonable feeling that experts here would come up
with some elegant suggestion to
Hello,
I had the same problem with R 2.11. After upgrading to R 2.12 it worked.
Hope it helps for you too.
--
View this message in context:
http://r.789695.n4.nabble.com/Can-not-invoke-maxent-of-library-dismo-in-GNU-linux-tp3242023p3304211.html
Sent from the R help mailing list archive at
Just add '.' to the pattern specifier:
MyString - ABCFR34564IJVEOJC3434.16ABC123.456KJHLKJH23452345AAA
# translate to the pattern sequences
x - chartr('ABCDEFGHIJKLMNOPQRSTUVWXYZ0123456789.'
+ , '00111'
+ , MyString
+ )
x.rle -
Hi all,
I have some problem with the tinnR editor. At
each end of the line, there is a strange symbol
and there is a dot between any two words.
Does anyone know how to correct this?
Thank you very much!
Hannah
[[alternative HTML version deleted]]
Dear all,
I have a data set which has latitude, longitude (not in increasing order)
and observed precipitation corresponding to those locations. I am using
image.plot command for that. But I want to overlay a map of US on that plot.
I know we can use map('usa') to draw it, but I dont know how can
Hi Gabor, thanks (and Jim as well) for your suggestion. However this is not
working properly for following string:
MyString - ABCFR34564IJVEOJC3434.36453
strapply(MyString, (\\D+)(\\d+)(\\D+)(\\d file://d+)(//d+)(//D+)(//d+),
c)[[1]]
[1] ABCFR 34564 IJVEOJC 3434
Therefore there is decimal
Hi,
I am trying to determine how many points fall ouside the confidence interval
range.
This is the code I have so far but it does not work. Any help would be
appreciated.
Count - vector ()
for (i in 1: nrow (dataname)){
if (dataname[i] l.ci.post[1]//
dataname[i] u.ci.post[i]){
count[i] - 1
Hi,
The logical operators are actually vectorized, so I do not think you
need a loop. Does this do what you want?
## Some data
set.seed(10)
dat - matrix(rnorm(500, sd = 3), nrow = 80)
## Hypothetical confidence interval
ci - c(-5, 5)
## Find the number of points outside interval
sum(dat
On Saturday, February 12, 2011 05:08:06 am beky wrote:
There is a code from SPSS Syntax
do if sub(ATVK,2,2)=01.
comp strata=1.
else if sub(ATVK,2,2)=05 and sub(ATVK,2,2)=27.
comp strata=3.
else if sub(ATVK,4,2)20 or sub(ATVK,6,2)20.
comp strata=4.
else if sub(ATVK,4,2)00.
comp strata=2.
Wow, great work! This makes a huge difference (i.e. 5 mins to process
an electrode instead of 5 hours).
I believe that the ndata and ndatam1 are there to let EMD::emd() and
EMD::extractimf() (the functions that call EMD::extrema()) implement
various boundary corrections. I personally don't intend
There are lots of modifications to make it faster.
E.g., replace the call to 'sequence' by a call to
'Sequence' and define Sequence by
Sequence - function(nvec) {
seq_len(sum(nvec)) - rep(cumsum(c(0L,nvec[-length(nvec)])), nvec)
}
That seemed to make the time for
testing the null hypothesis of no interaction is not the same as a
test of equivalence for the two differences. There is a literature on
tests of equivalence. First you must develop a definition of
equivalence, for example the difference is in the interval (-a,a).
Then, for example,
This is a seemingly simple problem, but I'm unable to make any progress. I
have a large dataset, within it there is a column of location names, some of
which repeat; the adjacent column lists values for each location name. I
need to sum the values which correspond to the same location name,
So do not use '$' or '[..]' terms in model formulae - this is going to cause
problems when it comes to predict, because your formula will not associate
with the names it has in its formula in the new data frame. When you think
about it, this is obvious.
I really don't understand why this is
Hi,
I understand that dichotimization of the predicted probabilities after
logistic regression is philosophically questionable, throwing out
information, etc.
But I want to do it anyway. I'd like to include as a measure of fit %
of observations correctly classified because it's measured in
On Feb 14, 2011, at 12:31 AM, Daniel Weitzenfeld wrote:
Hi,
I understand that dichotimization of the predicted probabilities after
logistic regression is philosophically questionable, throwing out
information, etc.
But I want to do it anyway. I'd like to include as a measure of fit %
of
Hi:
One (of many) ways is to use the aggregate() function in base R:
df - data.frame(g = rep(LETTERS[1:5], each = 5), y = rnorm(25))
aggregate(y ~ g, data = df, FUN = sum)
g y
1 A -1.3165949
2 B 1.9850015
3 C -0.6460323
4 D -1.6459293
5 E 6.7001807
HTH,
Dennis
On Sun, Feb 13,
Hey,
I am trying to graph a series of XY plots from a set of .csv files. Getting the
csv files read in and the original plot set up is not a problem, but I cannot
figure out a way to call get(iterator) within a separate dataset. A chunk of my
code is below.
data - read.csv(file= filename)
I believe it's actually 2/sqrt(n) where n is the sample size...
Interesting..
--
View this message in context:
http://r.789695.n4.nabble.com/autocorrelaltion-function-plots-tp3303966p3304490.html
Sent from the R help mailing list archive at Nabble.com.
Hi,
I have a huge table with 134,000 entries. I want to choose only rows 60,000
to 90,000.
Please help me on how to do this.
V1 V2 V3 V4
1 31.10267 168.3204 41 0.4
2 31.19941 168.1488 41 0.4
3 31.29580 167.9764 41 0.4
4 31.39183 167.8034 41 0.4
5 31.48750 167.6297 41 0.4
6
Hi,
Well, how did you make that little 10 row sample for us? The same
subsetting principles hold in R regardless of hugeness. I am kind of
assuming you have a matrix or data frame class object, in which case,
here is a little example for you:
## Big matrix
dat - matrix(1L, nrow = 134000, ncol
On Fri, Feb 11, 2011 at 9:26 PM, John Smith zmr...@gmail.com wrote:
Hello All,
I am using the following code to draw a figure. But the legend given buy
auto.key has mismatched colors. Could any one help me?
See ?simpleTheme
-Deepayan
I am using R2.12.1 and most current lattice on windows
Hi all,
Is there a way to save the currently displayed plot to an image format just
after we view it?
I think this would be more intuitive as a user if I wish to save it just
after I visualize the plot.
I am aware that we need to do some thing like this
jpeg('somefilename.jpg')
... plot...
look at
?savePlot
It will save the current plot you are viewing.
--
View this message in context:
http://r.789695.n4.nabble.com/saving-plots-tp3304540p3304548.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org
Thanks!
thats just what I was looking for!
regards,
KM
On Mon, Feb 14, 2011 at 1:16 PM, Joel joda2...@student.uu.se wrote:
look at
?savePlot
It will save the current plot you are viewing.
--
View this message in context:
http://r.789695.n4.nabble.com/saving-plots-tp3304540p3304548.html
71 matches
Mail list logo
