Hi
r-help-boun...@r-project.org napsal dne 19.01.2011 00:03:10:
On Jan 18, 2011, at 5:25 PM, kurt_h...@nps.gov wrote:
Greetings
I have a bunch of NAs in a column of categorical variables
designating
the size classes (e.g., smallest to largest: 1,2,3,4) of cave
crickets.
Hi:
Here's a variation on Jim Holtman's solution - same packages.
df - data.frame(V1 = rep(c('X1', 'X2'), each = 3),
V2 = rep(c('Y1', 'Y2', 'Y3'), 2),
Y1 = rep(1, 6), Y2 = rep(2, 6), Y3 = rep(3, 6))
dd - melt(df, id = c('V1', 'V2'))
dcast(dd, V2 + variable ~
tcc.new2 - droplevels(tcc.new)
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Ben Harrison
Sent: 19 January 2011 12:24
To: r-help@r-project.org
Subject: [R] Using subset to filter data table
I am having difficulty understanding
I get some results back from running an iterative analysis in the form of a
list of matrices. What I would like to do with this list is combine it such
that all the similar components get combined into a multidimensional array. If
possible I'd like to put results[[1]]$resultmean and
Hello list.
I feel like an idiot.
There exists a method called expand.grid which, from the documentation,
appears to do just what I want, but then it doesn't, and I can't get it to
behave.
Given a dataframe
dfr-data.frame(c1=c(a, b, NA, a, a), c2=c(d, NA, d, e, e),
c3=c(g, h, i, j,
Hi,
If you don't find droplevels(), then you should install the latest
version of R (2.12.1). It's always a good thing, even more if you're
starting.
Ivan
Le 1/19/2011 01:27, Ben Harrison a écrit :
Thanks for the reply Peter.
On 18 January 2011 22:52, Peter Ehlersehl...@ucalgary.ca
Well, my solution with the loop might be slower (even though I don't see
any difference with my system, at least with up to 100 lines and 3
strings to separate), but it works whatever the number of strings.
But I should have renamed the columns outside of the loop:
names(df)[2:3] - paste(a,
Hi:
Well, you fit a saturated model. How many degrees of freedom do you have
left for error? The fact that the standard errors are so huge relative to
the estimates is a clue.
Taking a look at your data, it's pretty clear that nation 3 is an
outstanding outlier on its own. It is clearly - nay,
On 2011-01-18 16:27, Ben Harrison wrote:
[...snip...]
plot(Ozone ~ Temp, data = aq)
This highlights one of the very confusing aspects of R language for me;
is plot(x, y) the same as plot (y ~ x)? Seems to be, but maybe I'm
missing some nuance.
plot has many methods; see with
Thank you very much Hugo for your answer.
Yesterday I was out of my office and I couldn't test the advise you
gave me. Today I'll do it.
I have never used C, so I have to ask.
For creating an executable file with your code, I copied and pasted
the text into a text file, and then tried from the
Hello everyone,
I have a data set like this:
head( fish_transect)
ID_TRANSECT ID_PROJECT DE_ZONE DE_LOCALITY DE_SECTOR MES
1 42 MB TarragonaCreixell Control I 9
2 42 MB TarragonaCreixell Control I 9
3 42 MB Tarragona
We are pleased to announce the launch of a new series of books on R.
Chapman Hall/CRC: The R Series
Aims and Scope
This book series reflects the recent rapid growth in the development and
application of R, the programming language and software environment for
statistical computing and
We are happy to inform you that abstract submission and registration for
useR! 2011 is now available online, see
http://www.R-project.org/useR-2011
This meeting of the R user community will take place at the University
of Warwick, Coventry, UK, August 16-18, 2011.
The conference schedule
Here is a more clear example:
# Create the ftable object
oftable-ftable(Titanic, row.vars = 4,col.vars=c(2,3))
require(gridExtra)
#Draw the grid, but the rows and heads are missing
grid.draw(tableGrob(oftable,
show.csep=TRUE, show.rsep=TRUE, show.box=TRUE,
Hi, all
I wanted to fetch data from Bloomberg for govt bonds, and analyse it
further.
I am having trouble in getting data as when I use field=PX_LAST, it is
giving the prices but when I use field=CPN, or ISSUE_DT, it is not giving
the results and just bouncing back NA for that.
This
Hi
r-help-boun...@r-project.org napsal dne 19.01.2011 09:51:43:
Hello everyone,
I have a data set like this:
head( fish_transect)
ID_TRANSECT ID_PROJECT DE_ZONE DE_LOCALITY DE_SECTOR MES
1 42 MB TarragonaCreixell Control I 9
2 42 MB
Hi:
Try this:
lapply(results, function(x) array(unname(unlist(x)), c(3, 3, 2)))
HTH,
Dennis
On Wed, Jan 19, 2011 at 12:31 AM, Maas James Dr (MED) j.m...@uea.ac.ukwrote:
I get some results back from running an iterative analysis in the form of a
list of matrices. What I would like to do
On 01/19/2011 04:04 AM, Andy Aldersley wrote:
Hello all,
Hoping that there is a fairly simple solution to my query...
I'm trying to overlay a line plot of some data onto a barplot of different
data. The y-axes are different for each set of data but the x-axes are the same
(1:12,
G'day Nick,
On Wed, 19 Jan 2011 09:43:56 +0100
Nick Sabbe nick.sa...@ugent.be wrote:
Given a dataframe
dfr-data.frame(c1=c(a, b, NA, a, a), c2=c(d, NA, d, e,
e), c3=c(g, h, i, j, k))
I would like to have a dataframe with all (unique) combinations of
all the factors present.
Easy:
R
slaps self in forehead/
I appear to have misinterpreted the help: considering that it explicitly
makes note of factors, I wrongly assumed that it would use the levels of a
factor automatically. My bad.
For completeness' sake, my final solution:
getLevels-function(vec, includeNA=FALSE,
Running ‘table’ on a factor with levels containing non-ASCII characters
seems to result in *extremely* bad performance on Windows. Here’s a simple
example with benchmark results (I’ve reduced the number of replications to
make the function finish within reasonable time):
library(rbenchmark)
Having upgraded to R version 2.12.1 I still have the same problem:
The combination of facet_grid and geom_hline produce (for me) 4 panels
of which two are empty of any data or lines (labelled 1 and 2).
Removing either the facet_grid or the geom_hline gives me the result I
would then expect.
I
Hi,
why does predict.lme() give missings in the predict.Subject column for
Subjects which are not in the original training set?
library(nlme)
example(predict.lme)
## c.f. Subject==F30 has missing values.
How can I predict/impute values for such new Subjects based on the fitted
model?
Thanks
Hello,
I have a file with very big values. I want to display the real values with
classic ways (without exposant)
summary(a[,1]);
Min. 1st Qu.Median Mean 3rd Qu. Max.
1.198e+09 1.199e+09 1.200e+09 1.200e+09 1.201e+09 1.202e+09
Can you help please?
Regards
Adel
--
Try this:
options(scipen = 20)
summary(a[,1])
On Wed, Jan 19, 2011 at 9:36 AM, Adel ESSAFI adeless...@gmail.com wrote:
Hello,
I have a file with very big values. I want to display the real values with
classic ways (without exposant)
summary(a[,1]);
Min. 1st Qu.Median
Hey everybody!
I have to use R's tm package to do some text analysis, first thing would be to
create a term frequency matrix.
Digging in tm's source code it seems like it uses some logic like this to
create term frequencies:
data(crude)
(txt - Content(crude[[1]]))
(tokTxt -
I got the following warning:
memorytest.c: In function ‘main’:
memorytest.c:5: warning: return type of ‘main’ is not ‘int’
Is this important ?
Hello Mauricio,
No, your gcc version is unduly puristic here. The traditional return
type of the main function in C should be int, and if that
Thanks a lot Petr!! It worked!
--
View this message in context:
http://r.789695.n4.nabble.com/count-of-factors-tp3224791p3224924.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
Hello R-help,
I recently spent a lot of time debugging a program, and the problem
turned out to be my own misuse of the det() function. A simple
recreation of the problem would be
X - diag (c(2,2,2))
det (X,log=TRUE)
## returns 8 as opposed to log 8
I had simply assumed that (along with many
Dear Hugo,
I tried your memory test program (without further modifications) just
after the gc() command:
print(gc())
print(gcinfo(TRUE))
system(/mypath/memorytest.out)
and the result that I got was:
Number of simulations read from ' Particles.txt ' : 9000
Dear R-users:
I am R-usser begginer and at the same time a beginner running
discriminant analysis;
I wanted to perform a DA using just the 80ß% of the original data but I
have some problems with simmilarity in variables,
Here my Skript
set.seed(123)
data80 - data[sample(472, 378), ]
Folks,
1. can you pelase cite the messages you are repsonding to?
2. Can you please include the original poster of the question who is not
necessarily subscribed to the list?
Thanks,
Uwe Ligges
On 18.01.2011 09:52, Joel wrote:
You could try useing '127.0.0.1' instead of 'localhost' and
Date: Wed, 19 Jan 2011 01:20:06 -0800
From: djmu...@gmail.com
To: laomen...@gmail.com
CC: r-help@r-project.org
Subject: Re: [R] question about result of loglinear analysis
Hi:
Well, you fit a saturated model. How many degrees of freedom do you have
left for error? The fact that the
X959X967X968
X9591 -0.04 0.45
X967-0.04 1 -0.09
X9680.45-0.09 1
X968_2 0.76-0.16 0.82
__
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 nueral networks and R. i want run the NN using nnet.
Following are my Questions.
1 To run regression using nnet, why should i specify linout=T always?. if
not it is not prediciting properly
2In the help, i have read that linout, entropy, softmax and censored are
mutually
Hi,
I am using nnet function to run regression and classification algorihtms.
following are my questions
1. test-nnet(Strength~.,cement,skip=T,size=0,linout=T), why linout=T always
for Regression algorithm
2. In the help, linout, entropy, softmax and censored are mutually
exclusive. but i can
Dear Mauricio,
what I do not understand at all is the message:
Error: cannot allocate vector of size 476.2 Mb
Have you tried to allocate a big matrix in between, whith
the R - statement not being shown in the output? Probably
not. If not, your local R version is buggy for sure. The test shows
Thanks again Hugo,
2011/1/19 Hugo Mildenberger hugo.mildenber...@web.de:
Dear Mauricio,
what I do not understand at all is the message:
Error: cannot allocate vector of size 476.2 Mb
Have you tried to allocate a big matrix in between, whith
the R - statement not being shown in the output?
Hello everyone,
I am trying to understand how covariance work. So I created a vector called
sr-c(2,5,7,5,2)
so according to wikipedia Cov(X,X)=E[XX]-E[x]*E[x] which in R is
mean(sr*sr)-mean(sr)*mean(sr)
[1] 3.76
but also
cov(sr,sr)
[1] 4.7
why is this difference between these two
Dear
R-users,
I was
wondering if someone could give me some advices on the following problem.
I tried to
apply moranâs test to a small dataset and couldnât succeed, here is the
error
message:
Â
mor - moran.test(x, res2)
Â
Avis dans
moran.test(x, res2) :
 Out-of-range p-value:
Hi Alex,
cov() uses the formula for a sample covariance. The denominator is N
- 1 instead of N. However, the formula you used from Wikipedia is for
a population covariance. To move between the two in this case:
cov(sr, sr) * (4/5) # should equal E[XX] - E[x] * E[x]
Cheers,
Josh
On Wed, Jan
Dear Josh,
I would like to thank you for your reply.
I think that it is clear that I miss a lot of theory. I have tried in google to
read and study more about covariance but I have this feeling that the term
'covariance' is used for many different things in different topics.
You mentioned
You could include a call to the source function that will read and run an R
script (I have not tried this, but don't see any reason that it would not work).
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
greg.s...@imail.org
801.408.8111
-Original
Hi Sandy,
It's difficult to know what's going wrong without a small reproducible
example (https://github.com/hadley/devtools/wiki/Reproducibility) -
could you please provide one? You might also have better luck with an
email directly to the ggplot2 mailing list.
Hadley
On Wed, Jan 19, 2011 at
That's a Bloomberg issue, not RBloomberg:
use bdp() (point) instead of bdh() (history) or bds() (set)
bdp(conn,ticker,c(PX_LAST,CPN))
PX_LAST CPN
912828JV Govt NA 0.875
912828JY Govt 100.0234 0.875
HTH,
David L. Reiner
XR Trading LLC
-Original Message-
From:
Hi,
On Wed, Jan 19, 2011 at 9:59 AM, Alaios ala...@yahoo.com wrote:
Dear Josh,
I would like to thank you for your reply.
I think that it is clear that I miss a lot of theory. I have tried in google
to read and study more about covariance but I have this feeling that the term
'covariance'
I just installed R 2.12.1, and when I went to run a few old programs with it,
nothing worked.
I got a ton of error messages saying such and such package was built before R
2.10.0 and needed to be reinstalled.
These were not just warning messages, but error messages that prevent the
programs
On Jan 19, 2011, at 9:49 AM, Joseph Boyer wrote:
I just installed R 2.12.1, and when I went to run a few old programs with it,
nothing worked.
I got a ton of error messages saying such and such package was built before R
2.10.0 and needed to be reinstalled.
These were not just warning
Hi,
I have missing values in my time series. na.action = na.pass works
for acf and pacf. Why do I get the following error for the ccf?
ts(matrix(c(dev$u[1:10],dev$q[1:10]),ncol=2),start=1,freq=1)
Time Series:
Start = 1
End = 10
Frequency = 1
Series 1 Series 2
1 68.0
You can use the updateusr function in the TeachingDemos package to get them to
match.
Sent from my iPod
On Jan 18, 2011, at 10:05 AM, Andy Aldersley ajalders...@hotmail.co.uk
wrote:
Hello all,
Hoping that there is a fairly simple solution to my query...
I'm trying to overlay a
Dear Joe,
On Wed, Jan 19, 2011 at 7:49 AM, Joseph Boyer joseph.g.bo...@gsk.com wrote:
I just installed R 2.12.1, and when I went to run a few old programs with it,
nothing worked.
I got a ton of error messages saying such and such package was built before R
2.10.0 and needed to be
population is all individual.
sample is a part of population.
If the sample size is enough large, mean(sr*sr)-mean(sr)*mean(sr) will
almost equal to cov(sr, sr).
On Wed, Jan 19, 2011 at 3:59 PM, Alaios ala...@yahoo.com wrote:
Dear Josh,
I would like to thank you for your reply.
I think that
Dear R users
Is there a way to obtain the residuals from a model fitted by quantile
regression? Thank you.
Thanaset
--
View this message in context:
http://r.789695.n4.nabble.com/Quantile-Regression-Extracting-Residuals-tp3225423p3225423.html
Sent from the R help mailing list archive at
I forgot to mention the plot looks like this :-(
http://img402.imageshack.us/i/volcanoi.jpg/
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting
I don't seem to have that method:
zts = as.ts.zoo(z)
Error: could not find function as.ts.zoo
I'm finding forum posts from people trying to use stl directly with zoo[reg]
functions (despite the documentation specifying ts) but discovering that stl
gives errors with NA entries. I encountered
Hello All,
I need to create a three-way interaction plot and I was wondering if anyone
could recommend a package I could use to design these plots.
Thank you!
PHRI DISCLAIMER
This information is directed in confidence solely to the person named above and
may
Hello all,
I was trying to make a volcano plot with some real data, using log2(ratio)
vs. Z-value significance.
However the scatter of the points is too less contrary to 'normal' volcano
plots and I'm getting a sharp 'V' shaped plot.
*Am I doing something wrong here?*
The data(ratio) is
Hi all,
I was trying to install a package called ei from this
http://gking.harvard.edu/eiR link . With GCC, r-base-dev and some other R
package dependencies installed it works perfect in my Ubuntu Linux system.
Now I am trying to install the same in Windows 7 environment. So far my
setup is
Hi all,
I understand that rmvdc generates random number from mvdc object. But the
mvdc object can only be used if we define the marginals! So my question is
suppose we don't find any distribution which fit marginals so we use the
Canonical Maximum Likelihood method (This approach uses the
Try looking here:
?rq.object
?residuals
R has excellent documentation that can answer many such questions with
less than a minute of reading.
--
Jonathan P. Daily
Technician - USGS Leetown Science Center
11649 Leetown Road
Kearneysville WV, 25430
(304)
Ah! I was always trying collapse with sep and other options. Not by itself.
Perfect!
And yes, that was my bad example.
Robin Jeffries
MS, DrPH Candidate
Department of Biostatistics
UCLA
530-624-0428
On Tue, Jan 18, 2011 at 10:27 PM, Joshua Wiley jwiley.ps...@gmail.comwrote:
Hi Robin,
I have a loop that I would like to use to extract the stddev for
each itteration so I can average the stddev for all the runs. It
would be helpful to know how to extract the stddev for each run from
the VarCorr. Thanks
MCruns-1000
sighatlvec-rep(NA,MCruns)
sighatbvec-rep(NA,MCruns)
Hello everyone,
I'm quite new to R and am trying to run a logistic model to look at how various
measures of boldness in individual animals influences probability of capture,
however I also want to include random effects and I'm not sure how to construct
a model that incorporates both of these
Did you do:
update.packages(checkBuilt=TRUE)
Rob
--
From: Joshua Wiley jwiley.ps...@gmail.com
Sent: Wednesday, January 19, 2011 11:43 AM
To: Joseph Boyer joseph.g.bo...@gsk.com
Cc: r-help@r-project.org
Subject: Re: [R] how to get old packages to
On Wed, Jan 19, 2011 at 1:10 PM, Adam Oliner oli...@gmail.com wrote:
I don't seem to have that method:
zts = as.ts.zoo(z)
Error: could not find function as.ts.zoo
I'm finding forum posts from people trying to use stl directly with zoo[reg]
functions (despite the documentation specifying ts)
Dear Stephanie,
I'm not sure what kind of plot you have in mind, but you might take a look
at the effects package, which can plot interactions for a variety of models.
I hope this helps,
John
John Fox
Senator William McMaster
Professor of Social Statistics
On Wed, Jan 19, 2011 at 3:07 PM, Adam Oliner oli...@gmail.com wrote:
This still does not give me the desired behavior:
z = zoo(x, order.by=t, frequency=24)
t = as.ts(z)
... inserts 23 NA values between every actual value. This is not correct;
the original data has a frequency of 24 and
Dear R People:
I'm building R-2.12.1 from scratch on a Windows XP machine.
The make all recommended worked fine, but when I did the make in the
bitmap directory, I got the following error:
C:\R\R-2.12.1\src\gnuwin32\bitmapmake
make
make CC='gcc -std=gnu99' AR='ar' \
CFLAGS=-O3
Cant you just save that column of values to an object
stddev-VC[,2]#if it's the second column that has the values
also I wouldn't save your model name in the same name as the command, when
you have lmer-lmer(y~1+(1|sampfac)+(1|labfac)), may be best to have
Hello again.
Libpng version 1.5.0 does not have a scripts/makefile.mingw.
However, Libpng version 1.2.44 does. That seems to be the snag.
Presumably I should stick with the older version.
Erin
--
Erin Hodgess
Associate Professor
Department of Computer and Mathematical Sciences
University
Is this what you want?
attr(VarCorr(merobject)[[1]], stddev)
Alternately, look at:
getMethod(VarCorr, signature = mer)
For what is happening behind the scenes and why I chose what I did to extract
the standard deviations
Cheers,
Josh
On Jan 19, 2011, at 11:32, Stephanie Cooke
This still does not give me the desired behavior:
z = zoo(x, order.by=t, frequency=24)
t = as.ts(z)
... inserts 23 NA values between every actual value. This is not correct;
the original data has a frequency of 24 and doesn't need one forced upon it
during coercion.
z = zoo(x, order.by=t)
t =
Thanks Greg,
work i'm lost a { after function ().
Alessandro
2011/1/19 Greg Snow greg.s...@imail.org:
You could include a call to the source function that will read and run an R
script (I have not tried this, but don't see any reason that it would not
work).
--
Gregory (Greg) L. Snow
Gabor,
If you follow those steps on a series with missing values, the resulting tt
will contain tons of new NAs. (See my original email and the call to zrt =
as.ts(zr).)
At any rate, I finally managed to finagle a ts with the correct frequency by
coercing to ts _twice_ and specifying the
On Wed, Jan 19, 2011 at 4:17 PM, Adam Oliner oli...@gmail.com wrote:
Gabor,
If you follow those steps on a series with missing values, the resulting tt
will contain tons of new NAs. (See my original email and the call to zrt =
zoo and zooreg classes handle series which may not be regularly
Hi
Still having problems in that when I use geom_hline and facet_grid together I
get two extra empty panels
A reproducible example can be found at:
https://gist.github.com/786894
Sandy Small
From: h.wick...@gmail.com [h.wick...@gmail.com] On Behalf Of
I suspect you have libpng 1.5.0 (assuming you installed libpng at
all). So
1) That version has a CERT advisory, and should not be used. 1.5.1 is
due tomorrow.
2) R 2.12.1 cannot work with future versions of external software.
libpng 1.4.5 does work, but not 1.5.x.
(And BTW, questions
Hi all,
I had a weird problem with xyplot and I am wondering if anyone can help me
figure out what's wrong here. Basically if I call xyplot() inside a function
and I call quartz() after, then the previous xyplot() does not do anything. A
short demo of the issue is like this:
library(lattice)
On Wed, Jan 19, 2011 at 4:27 PM, Adam Oliner oli...@gmail.com wrote:
My series was regularly spaced, it simply contained missing values. I don't
see why the new values were added.
Besides, my previous email showed exactly how to convert a zoo object with
Please provide a minimal
LAURA WEIR 0611792W at student.gla.ac.uk writes:
Hello everyone,
I'm quite new to R and am trying to run a logistic model
to look at how various measures of boldness in
individual animals influences probability of capture,
however I also want to include random effects and
I'm not
On 11-01-19 4:54 PM, Saiwing Yeung wrote:
Hi all,
I had a weird problem with xyplot and I am wondering if anyone can help me
figure out what's wrong here. Basically if I call xyplot() inside a function
and I call quartz() after, then the previous xyplot() does not do anything. A
short demo
Hi all
I'm plotting colored contour maps using filled.contour. My levels
are very unevenly spaced, with, say, high resolution in the small
numbers but ranges that can be an order of magnitude or two larger in
absolute value compared to where the action takes place. Aside from
transforming the
Hi Sandy:
I can reproduce your problem given the data provided. When I change
ecd_rhythm from character to factor, it works as you intended.
str(lvefeg)
List of 4 ### Interesting...
$ cvd_basestudy: chr [1:10] CBP05J02 CBP05J02 CBP05J02 CBP05J02 ...
$ ecd_rhythm : chr [1:10] AF AF AF AF
I have an mcmc object and I'm trying to plot the quantiles of the variables -
and not as a function of the iterations as in cumuplot.
I cannot seem to find the right combination of indexing to access the
variables; after which I'm sure I can plot all the statistics I could hope for.
Any hints
Hello,
I
am using the Kinhom function to plot envelopes from my fitted
inhomogeneous models in spatstat. But the results I obtain do not have a
theoretical value (theo). They only give me the hi, lo, r and
mmean values but no theo. Please can anyone help?
Mike.
[[alternative
My series was regularly spaced, it simply contained missing values. I don't
see why the new values were added.
Besides, my previous email showed exactly how to convert a zoo object with
missing values into a ts object with missing values and a frequency other
than 1, so I'm not sure what you
The objective is to get the stl calls to work without making up values for
the NAs:
library(zoo)
data = c(50, 32, 22, 13, 14, 6, 12, 3, 6, 21, NA, NA, 15, 13, 17, 9, 11, 34,
63)
dataz = zoo(data, 1:length(data))
t = as.ts(dataz)
datat = as.ts(dataz)
datatt = ts(t, frequency=2)
frequency(datatt)
This is a version of FAQ 7.22: you need to tell R to print lattice graphics.
Without the next statement, the result would be auto-printed.
Duncan Murdoch
That fixed it. Thank you again. Sorry for not reading the FAQ.
Saiwing
__
On Wed, Jan 19, 2011 at 5:21 PM, Adam Oliner oli...@gmail.com wrote:
The objective is to get the stl calls to work without making up values for
the NAs:
library(zoo)
data = c(50, 32, 22, 13, 14, 6, 12, 3, 6, 21, NA, NA, 15, 13, 17, 9, 11, 34,
63)
dataz = zoo(data, 1:length(data))
t =
Dear Greg Snow,
I'm a biologist trying to write a mathematical formula for a doubly
truncated normal distribution to be used in the language R. I realise this
is simple stuff for a mathematician but I'm stumped.
Wikipedia gives what seems a fairly simple formula - with function = maths
with
Thank you Santosh.
I am so new to R that I don't even have droplevels:
1 tcc2 = droplevels(tcc)
Error: could not find function droplevels
I looked up ?droplevels:
1 ?droplevels
No documentation for 'droplevels' in specified packages and libraries:
you could try '??droplevels'
1 ??droplevels
Help
Gabor,
On Wed, Jan 19, 2011 at 3:36 PM, Gabor Grothendieck ggrothendi...@gmail.com
wrote:
stl does not support series with NAs
Thank you. This is what I needed to know.
Cheers,
--
- Adam J. Oliner
[[alternative HTML version deleted]]
Whoops, fixed a typo:
library(zoo)
data = c(50, 32, 22, 13, 14, 6, 12, 3, 6, 21, NA, NA, 15, 13, 17, 9, 11, 34,
63)
dataz = zoo(data, 1:length(data))
datat = as.ts(dataz)
datatt = ts(datat, frequency=2)
frequency(datatt)
length(datatt)
# these all give errors
stl(datatt, s.window=per)
# this next
On Wed, 19 Jan 2011, JClark wrote:
Dear Greg Snow,
I'm a biologist trying to write a mathematical formula for a doubly
truncated normal distribution to be used in the language R. I realise this
is simple stuff for a mathematician but I'm stumped.
Wikipedia gives what seems a fairly simple
On Jan 19, 2011, at 6:03 PM, JClark wrote:
Dear Greg Snow,
I'm a biologist trying to write a mathematical formula for a doubly
truncated normal distribution to be used in the language R. I
realise this
is simple stuff for a mathematician but I'm stumped.
Wikipedia gives what seems a
Hello everyone,
I'm stumped. I'd like to create a scatterplot matrix with circular
reference lines. Here is an example in 2d:
library(ellipse)
set.seed(1)
dat - matrix(rnorm(300), ncol = 3)
colnames(dat) - c(X1, X2, X3)
dat - as.data.frame(dat)
grps - factor(rep(letters[1:4], 25))
panel.circ -
Hi,
It seems simple but I could not find the solution on the R site.
I can't save properly .RData and .Rhistory in Windows 7 folders. I get
files with no extension.
Thanks in advance for your help.
JP
Jean-Philippe Puyravaud
9A Frederic Osanam Street
Colas Nagar
Thank you Santosh.
I am so new to R that I don't even have droplevels:
1 tcc2 = droplevels(tcc)
Error: could not find function droplevels
I looked up ?droplevels:
1 ?droplevels
No documentation for 'droplevels' in specified packages and libraries:
you could try '??droplevels'
1 ??droplevels
Help
Ben -
droplevels was introduced in R version 2.12.0.
Since you are running version 2.10.1, it won't be
available. You might want to upgrade to a more
recent version.
- Phil Spector
Statistical Computing
I think you get files with extensions, but you have not told Windows to show
them to you. Adjust your file display settings in a file manager window.
Jean-Philippe Puyravaud jp.puyrav...@gmail.com wrote:
Hi,
It seems simple but I could not find the solution on the R site.
I can't save
100 matches
Mail list logo
