On Thu, 8 Jan 2009, Neil Beddoe wrote:
You can also use subscripts to get at things with a bit of playing around.
You can, but it is easier not to fight R to do so. Much more transsparent
is:
cf - coef(summary(lm fit)
cf[2,2] # index as a matrix.
You are
a) indexing a matrix as a vector
Hello,
I am unable to update packages. I am running R version 2.8.1 (2008-12-22)
under VISTA. When I start R I get the following error:
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
Error in
Dear R Help,
I wonder the way to show the source code of [vcov] command. Usually, it
can show the source code after input the command and enter. But for
[vcov], it shows
function (object, ...)
UseMethod(vcov)
environment: namespace:stats
I appreciate for your help. Best wishes.
On Thu, 2009-01-08 at 10:39 +0530, Kishore wrote:
Hi,
I have good understanding on Econometrics and statistical techniques.
However, I am new to R. What would be the best way to learn R as I would
be one of the few in my team started exploring R in your team. I have got a
few downloads on
Hello,
You can do
stats:::vcov.lm
to see the source code for that particular method. In order to see which
are the methods supported by vcov, write
methods(vcov)
Best regards,
Carlos J. Gil Bellosta
http://www.datanalytics.com
On Wed, 2009-01-07 at 21:37 -0600, Yang Wan wrote:
Dear R
Hello,
You are not very precise there. Do you mean that the rows in your text
file do not all have the same number of separators (commas, in your
case)?
Best regards,
Carlos J. Gil Bellosta
http://www.datanalytics.com
On Thu, 2009-01-08 at 04:38 -0500, rahul-a.agar...@ubs.com wrote:
I have a
Simon Blomberg wrote:
Some people familiar with R describe it as a supercharged version of
Microsoft's Excel spreadsheet software...
Maybe it's my dry Australian humour, but I think this should go into the
fortunes package.
It is somewhat humourous even in damper London.
But I think
Hi Rahul,
try to see
?read.table
you may substitute commas with a NA
and then sum the NA.
I hope that it can help.
Cheers Anna
Anna Freni Sterrantino
Ph.D Student
Department of Statistics
University of Bologna, Italy
via Belle Arti 41, 40124 BO.
Please do read the rw-FAQ (as requested in the posting guide): it will
enlighten you.
On Wed, 7 Jan 2009, Greg Blevins wrote:
Hello,
I am unable to update packages. I am running R version 2.8.1 (2008-12-22)
under VISTA. When I start R I get the following error:
I presume that is
Dear R helpers
Suppose I have a portfolio of securities with exposure to Equity, Bonds and
Forex (say $ 100 each).
Is there any fucntion in R that will help me calculate Value at Risk (VaR)
using Monte carlo Simulation , Historic simulation and Variance - Covariance
Simulation.
With
Yes, there are: replicate and quantile are your friends.
You will find better support in the R-Finance list, though.
Best regards,
Carlos J. Gil Bellosta
http://www.datanalytics.com
On Thu, 2009-01-08 at 01:36 -0800, Maithili Shiva wrote:
Dear R helpers
Suppose I have a portfolio of
On Tue, Apr 29, 2008 at 7:56 AM, Ted Harding
ted.hard...@manchester.ac.uk wrote:
Certainly for Windows users, the suggested change could be
wise, since Windows (for reasons best known to someone who
should have known better) conceals the extension from the
I think the reason for hiding
Dear all,
I found that printing with 'cat' is very slow.
For example in my machine this snippet
__BEGIN__
# I need to resolve to use this type of loop.
# because using write(), I need to create a matrix which
# consumes so much memory. Note that foo, bar, qux object
# is already very large
On Thu, Jan 8, 2009 at 3:12 AM, Gundala Viswanath gunda...@gmail.com wrote:
Dear all,
I found that printing with 'cat' is very slow.
For example in my machine this snippet
__BEGIN__
# I need to resolve to use this type of loop.
# because using write(), I need to create a matrix which
#
Barbara.Spillmann at agrar.uni-giessen.de writes:
dear R users,
I have troubles performing Moran.I test as suggested on
http://www.ats.ucla.edu/stat/r/faq/morans_i.htm
Using the ape package which is written for a different application area is
not necessarily a good idea - the web page
What exactly is the problem you are trying to solve. What is going to
be done with the data? Is it going to be read by some other program?
How much physical memory do you have on your machine? Is there paging
occuring due to the size of the objects? Have you consider creating a
structure with
Dear Jim and Henrik,
What exactly is the problem you are trying to solve.
Is it going to be read by some other program?
I simply want to print the data out. Surely, this data
will be manipulated (with Excel or other
programming languages) by other people suit to their purpose.
Typically the
Hi
r-help-boun...@r-project.org napsal dne 08.01.2009 14:26:32:
Dear Jim and Henrik,
What exactly is the problem you are trying to solve.
Is it going to be read by some other program?
I simply want to print the data out. Surely, this data
will be manipulated (with Excel or other
on 01/07/2009 09:47 PM Gabor Grothendieck wrote:
On Wed, Jan 7, 2009 at 10:26 PM, Dirk Eddelbuettel e...@debian.org wrote:
On 7 January 2009 at 18:24, Gabor Grothendieck wrote:
| By running the code below we see that the:
| - sum of the three seems to be rising at a constant rate
| - S is
I doubt that you have a dataframe with those features, since R would
not allow such an event; more likely you have data in a file. If your
goal is to determine the number of items, then you should definitely
look at:
?count.fields
count.fields(filename, sep=,) # would give you 1 + the
Here is one way of doing it. To write out 1 million rows on my system
took 21 seconds.
# create some data
dataSize - 1e6
foo - runif(dataSize)
bar - runif(dataSize)
n - 1000 # number of items to write out each time
output - file('/output.txt', 'w')
# now split the indices into groups of
Dear R-help,
I am reading strings (varchar) from MySQL database
and instead of some characters get garbled chars.
Strings are encoded in UTF-8. I have tried to change
encoding in MySQL to latin7, but it does not help.
I was not able to find solution, how to get strings
in proper encoding/form.
Browse[1] j - c(1,2,NA)
Browse[1] j[j==1][-1]
[1] NA
Browse[1] j[j==1][-2]
[1] 1
Browse[1] j[j==1][-2] - 2
Error during wrapup: NAs are not allowed in subscripted assignments
As far as I can see, I have no NA in the lhs (not after the second
subscript anyway).
Besides, I have a single value on
More commentary on Slashdot:
http://developers.slashdot.org/article.pl?sid=09/01/07/2316227
__
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
Hi all,
I would like to upload some plots I create wth R via ftp or something
similar but I don t really understand which command/syntax I should
use:
should I go for make.socket + write.socket or try and create
environment variables like frp_proxy_user and then ftp my files? how?
many thanks
Thanks Duncan, I'll try that with my data.
--- On Wed, 1/7/09, Duncan Mackay mac...@northnet.com.au wrote:
From: Duncan Mackay mac...@northnet.com.au
Subject: Re: [R] xtable-longtable question
To: mazatlanmex...@yahoo.com
Cc: r-help@r-project.org
Date: Wednesday, January 7, 2009, 1:15 PM
I
The open-source mentality is invaluable, as most on this list know. That
is what keeps the R evolution progressing at a pace that SAS cannot keep
up with.
On a side note (a very side note), I am a zealot for an exercise program
called Crossfit. Crossfit has adopted the same open-source mentality
Try system() with curl or a decent ftp client (I don't see that package
RCurl covers this, but it might despite its description only mentioning
HTTP). From 'man curl'
curl offers a busload of useful tricks like proxy support, user authen-
tication, ftp upload, HTTP post, SSL
Readers,
I have looked at various documents hosted on the web site; I couldn't
find anything on interpolation. So I started r and accessed the help
(help.start()). (by the way is it possible to configure r to open help
in opera instead of firefox?) Initially I read the help for the akima
package
Daren Tan dare...@hotmail.com wrote:
I have developed some scripts that basically ask for input
tab-limited format files, do some processing, and output several
pictures or csv. Now I need to have some gui to wrap on top of
the scripts, [...] Please advice me if there is any tools or
project
On 1/7/09, Gabor Grothendieck ggrothendi...@gmail.com wrote:
Here is the same number of messages/posts data
for each of S, SAS, R:
- reworked into a 3 column ts class time series
- with Jan 2009 removed since its not complete
- leading and trailing NA rows removed
My software of choice
Many thanks for that.
I am running R under windows and I am not sure what you mean by this curl?
Is'nt there a way to ftp via an intrinsic R function?
Thanks again
Thomas
2009/1/8 Prof Brian Ripley rip...@stats.ox.ac.uk:
Try system() with curl or a decent ftp client (I don't see that package
Dear:
I am trying to plot the histogram graph for my observed data. Then plot
fitted distribution on the same graph of histogram plot in R.
1. histogram plot y.
2. based on 1, plotting y1 v. x;
3. based on 1, plotting y2 v. x;
4. based on 1, plotting y3 v. x;
All of
Hi!
I have been working several years with R but it's my first public question.
I hope I'll be clear :) .
This question is related to obtaining letter-based representation of
non-parametric pairwise comparisons.
I have a dataframe with this structure (but with quite more rows and cols):
A B C
Thank you. It worked!
Quoting Henrique Dallazuanna www...@gmail.com:
Try this:
reshape(d, timevar = x, idvar = y, direction = wide)
On Tue, Jan 6, 2009 at 6:29 AM, poas...@umich.edu wrote:
Dear Sir or Madam,
I have the following data frame (which is just a toy example of my larger
Hallo,
I didn´t found any facilities for Halbergs cosinor analysis in R. This analysis
is well known in the Chronobiology as the least square approximation of time
series using cosine function of known period (in my case of 24hours-period). I
tried to write a script but crashed...
Can you give
Hi,
You may want to try the double exponential transformation on the numerator and
the denominator on this one.
The method is described in detail here:
http://projecteuclid.org/DPubS?service=UIversion=1.0verb=Displayhandle=euclid.prims/1145474600
If you want to give it a shot outside R
Hi all,
I would like to upload some plots I create wth R via ftp or something
similar but I don t really understand which command/syntax I should
use:
should I go for make.socket + write.socket or try and create
environment variables like frp_proxy_user and then ftp my files? how?
many thanks
If you are on Windows (you didn't say!) see also ?winMenuAdd, ?winDialog
?choose.files etc.
-- Bert Gunter
Genentech
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Mike Prager
Sent: Thursday, January 08, 2009 8:29 AM
To:
Some people familiar with R describe it as a supercharged version of
Microsoft's Excel spreadsheet software...
It is easy to ridicule this line from the NYT article. But this is not only
a very sensible comment by a smart reporter, but also one that is good for
R:
It is good for R because
On Thu, Jan 8, 2009 at 5:42 PM, Thomas Loridan thomas.lori...@gmail.com wrote:
Many thanks for that.
I am running R under windows and I am not sure what you mean by this curl?
http://curl.haxx.se/
Is'nt there a way to ftp via an intrinsic R function?
That I don't know, perhaps some else
Using indexing and putting a minus sign in front of a vector of column
names that you want to exclude would be a typical approach:
df - data.frame(a=LETTERS[1:4], b= rnorm(4), c=rnorm(4), d=
letters[5:9])
apply(df[ , -c(a,d)], 2, sum)
(Pretty sure this will run properly but don't have R
On Thu, 8 Jan 2009, Gábor Csárdi wrote:
On Thu, Jan 8, 2009 at 5:42 PM, Thomas Loridan thomas.lori...@gmail.com wrote:
Many thanks for that.
I am running R under windows and I am not sure what you mean by this curl?
http://curl.haxx.se/
And it comes as standard on Mac OS X and many Linux
On Thu, Jan 8, 2009 at 5:40 AM, retama retama1...@gmail.com wrote:
Hi!
I have been working several years with R but it's my first public question.
I hope I'll be clear :) .
This question is related to obtaining letter-based representation of
non-parametric pairwise comparisons.
I have a
Hi:
Could someone try to run this example? I can't make it work.
## Not run:
\begin{small}
\setlongtables
\begin{longtable}{
results=tex,fig=FALSE=
cat(paste(c('c', rep('cc', 34/2-1), 'c'), collapse='@{\hspace{2pt}}'))
@
}
\hline
\endhead
\hline
\endfoot
results=tex,fig=FALSE=
On Thu, 8 Jan 2009, Rune Schjellerup Philosof wrote:
Browse[1] j - c(1,2,NA)
Browse[1] j[j==1][-1]
[1] NA
Browse[1] j[j==1][-2]
[1] 1
Browse[1] j[j==1][-2] - 2
Error during wrapup: NAs are not allowed in subscripted assignments
As far as I can see, I have no NA in the lhs (not after the second
You are plotting the histogram in the frequency scale. A quick look
at the doc page for hist() would reveal the freq option:
hist(x,freq=FALSE)
then you can add the densities with lines()
albyn
Quoting Xin Shi jasonshi...@hotmail.com:
Dear:
I am trying to plot the histogram
In addition to the other excellent examples already given, I maintain
a list of resources for newcomers to R here:
http://www.revolution-computing.com/community/resources.php#TOC-Resources-for-Newcomers-to-R
I can also suggest this blog category (which will grow further with time):
Hello - and happy newyear to all of you!
I've got some data that I'm plotting with bwplot, a 3x2x3 design where
the observable decreases with the principle independent factor, but at
different rates.
I'd like to get lattice to impose not a single set of axes ranges
identical for all panels, but
Even using the VBA back of Excel to create interfaces with R would
make a lot of sense. Suppose I could have access to VBA macros that
import and export data into R , it would be great.
The R GUI series like Rattle come even closer to Excel...so a VBA
_R_ExCel package might be useful to ordinary
On 1/8/09, Maithili Shiva maithili_sh...@yahoo.com wrote:
Is there any fucntion in R that will help me calculate Value at Risk (VaR)
using Monte carlo Simulation , Historic simulation and Variance - Covariance
Simulation.
There are some un-published Crystal Ball functions for R [1], which
Prof Brian Ripley wrote:
Try system() with curl or a decent ftp client (I don't see that package
RCurl covers this, but it might despite its description only mentioning
HTTP).
It does support FTP, and all of the protocols that are supported in
the installed libcurl, so it depends the
There is RExcel (available by downloading the CRAN package
RExcelInstaller. It allows to transfer data between R and Excel,
and run R code from within Excel. So you can start with your data in
Excel, let R do an analysis, and transfer the results back to Excel.
You can write VBA macros which do
If you want to just linearly interpolate, then use the functions approx or
approxfun from the stats package (one of those that is loaded by default). See
the function TkApprox in the TeachingDemos package for an interactive way to
plot the approximations with the interpolations plotted.
# Unfortunately due to a small bug in the current version of ggplot,
you need to use
# jpeg with an e to get jpegs
ggsave(file = ggsave_qplot_output.jpeg)
# But I'd recommend using png instead - much better quality
ggsave(file = ggsave_qplot_output.png)
Here's a cute illustration of why I
On Thu, Jan 8, 2009 at 5:52 AM, Christian Kamenik
christian.kame...@giub.unibe.ch wrote:
'Apply' is a great thing for running functions on rows or columns of a
matrix:
X - rnorm(20, mean = 0, sd = 1)
dim(X) - c(5,4)
apply(X,2,sum)
Is there a way to use apply for excluding rows or
Hi,
I have an vector object that looks like
DA - c(1991q1, 1993q2)
(first quarter of 1991 etc) and I want to convert it into a date object
using as.Date(). I did this for montly data but am stumped when it comes to
dealing with quarterly data and as.Date.
Would anyone be able to help?
Hi Erich,
I saw that it uses a remote server ( which can be the same machine ) to
compute.
Here is the question-
What is the remote server is Amazon EC2 which has upscalaing and downscaling
facillity for RAM and CPU...
Will it work ?
is there a SaaS version of this?
Regards,
Ajay
On Thu,
Hi Erich,
I would like to share and embed the RExcel Training video (just like youtube
allows me to) . How can I do that ?
Regards,
Ajay
www.decisionstats.com
On Fri, Jan 9, 2009 at 12:20 AM, Ajay ohri ohri2...@gmail.com wrote:
Hi Erich,
I saw that it uses a remote server ( which can be
Can't you just set a link to our server?
I do not know how to embed it.
I used Camtasia to produce this video, and it
does all the embedding.
Ajay ohri wrote:
Hi Erich,
I would like to share and embed the RExcel Training video (just like
youtube allows me to) . How can I do that ?
Here is a direct quote from the help page for hist:
if 'FALSE', probability densities, component 'density', are plotted (so that
the histogram has a total area of one).
That looks like it may apply to your problem, just figure out what the above
quote is referring to. The last example on that
On Thu, 8 Jan 2009, Anne Berger wrote:
Hallo,
I didn?t found any facilities for Halbergs cosinor analysis in
R. This analysis is well known in the Chronobiology as the least
square approximation of time series using cosine function of known
period (in my case of 24hours-period). I
Well, testing would have been wise.. The last variable in the sample
dataframe had 5 elements and the use negation of a character vector
is not proper.
df[ , -c(a,d)]
Error in -c(a, d) : invalid argument to unary operator
So if you limit yourself to negative indexing of numeric references
hi,
i was wondering if i have been giving ANOVA table in R:
Response:MPG
DF Sum SqMean Sq F vaule
Pr(F)
Model 1 216.750216.7506.1272
0.04811*
Model.Mixture4
As the product manager for S+, I'd like to comment as well. I think the
burgeoning interest in R demonstrates that there's demand for analytics
to solve real, business-critical problems in a broad spectrum of
companies and roles, and that some of the incumbent analytics offerings,
in particular
On Thu, 2009-01-08 at 10:42 -0600, Stas Kolenikov wrote:
A really good measure for R will be the total # of the downloads of
r-base for all platforms from all CRAN mirrors (and I would expect
that # can be found from the servers' logs).
Hello,
You obviate here that many of us are downloading
Yes I think R as a package can really learn from SAS and SPSS in
making GUI more user friendly , even at the risk of dumbing down some
complexity..
also as a consultant I know that selling software requires a lot of
marketing follow ups..which is why R has lagged behind in actual
implementation
I believe R as a package has everything people with little knowledge of
programming can handle quite easily. Moreover even if someone has no
programming knowledge can learn R without much effort.
I also believe if people in corporate world start using R instead of
other complex software which are
You can use zoo to do this:
library(zoo)
DA - c(1991q1, 1993q2)
as.Date(as.yearqtr(DA))
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
On Behalf Of Shruthi Jayaram
Sent: Thursday, January 08, 2009 1:46 PM
To: r-help@r-project.org
Subject: [R]
on 01/08/2009 01:12 PM Andrew Choens wrote:
On Thu, 2009-01-08 at 10:42 -0600, Stas Kolenikov wrote:
A really good measure for R will be the total # of the downloads of
r-base for all platforms from all CRAN mirrors (and I would expect
that # can be found from the servers' logs). Given that it
On Thu, 2009-01-08 at 13:52 -0600, Marc Schwartz wrote:
Reading the posts on SAS-L since yesterday via Google RSS, where the
NYT
article was also posted, some have noted that SAS itself offers online
support forums (http://support.sas.com/forums/index.jspa). From a
quick
review, it looks
Thank you all! In future I will include examples of my code to make things
simpler for you. This is what I settled on:
Sp.presence-Data[,14:31] ##The subset of my data set I'm interested in
(the presence/absence data)
Sp.presence$Species-apply(Sp.presence,1,function(x)
Hi:
I am going through some of the xtable examples and I can't make the one below
work. I need to create a longtable on the fly keeping the column headers for
all the pages and I thought this example could give some ideas on how to do it.
I am using Sweave and xtable to create my tables and
Thank you again for your response.
This worked great.
Quick question about the legend for qplot. Instead of being outside the plot,
is it possible to move the location of the legend to the upper left or right
corner of the plot? Could you possibly provide an example.
Thank you again
Hello,
I would like to add the 'equivalent' sign (i.e., three horizontal lines,
almost like an equal sign) to a plot. Does R have anything like the LaTeX
command '\equiv' that I could put into a plot title?
Thanks
--
View this message in context:
Try
plot(1:10, main=expression(x%==%y))
For more details see
?plotmath
or
demo(plotmath)
HTH ...
Peter Alspach
-Original Message-
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On Behalf Of jjh21
Sent: Friday, 9 January 2009 10:07 a.m.
To:
Dont be so American
Cut the paper work
lets get going
[[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
I'm having trouble with alignment of a trend line overlayed onto an image
plot. The two should be plotted on the same x-axis (time-series). However,
the trend line begins about an inch into the image plot x-axis and ends
about an inch off of end of the image plot. Once I have the alignment
The updateusr function in the TeachingDemos package can be used to change the
user coordinates of the first plot to match what you think they should before
adding to the plot.
Hope this helps,
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
greg.s...@imail.org
Cheers Marc, it works perfectly now. Thanks for your help!
Marc Schwartz wrote:
You can generalize the approach by using something like:
...
ylim = c(0, max(DF$TACC, DF$Catch) * 1.1)
...
That would allow you to use the max value of the two columns, multiplied
by a fudge
I have a function that is called from a sourced script, whenever
certain conditions are met.
I have this inside it: if(P 0.9)save( file = today() , myLMSmodel );
This gives me files that look like this: Thu Jan 8 14:04:43 2009
That works fine on my Mac, but it would be much more convenient
Unfortunately a little long
Try this (I have used Miktex 2.7 on win XP and R 2.8 for this)
I have not gone into the detail of the first chunk I just created the
header strings in R and copied to make the Sweave file lines
With a little thought it can be done I suppose
## Sweave script
Gregg
You could either manipulate the character returned by today(), or try
using Sys.time() instead. For example, on Windows,
paste(strsplit(date(), ' ')[[1]][c(2,3,5,4)], collapse='_')
gives
Jan_09_2009_13:40:57
which approaches what you want. Alternatively
format(Sys.time(),
Hi, I created the script below to plot multiple graphs. It seems to work
pretty well (see BYX7) except when the max value for catch or TACC is very
small. In this case the x-axis labels are removed from the x-axis (see
BYX10). I can't figure out why given that y co-ordinate for the text
Your easiest starting point is probably
file = format(Sys.time())
and then look at the online help for the strptime function to learn
about formatting options.
For example
format( Sys.time() , '%m%d%y' )
Use paste() and some formatting of Sys.time() to construct the form
you are
Is anyone in the leadership of the R-project going to contact the New York
Times and clarify that the article gave remarkably short shrift to the
people who designed the user interface for R, to a large extent ATT
researchers from an earlier generation? It would be the appropriate thing to
do.
Is anyone in the leadership of the R-project going to contact the New
York Times and clarify that the article gave remarkably short shrift
to the people who designed the user interface for R, to a large extent
ATT researchers from an earlier generation? It would be the
appropriate thing to do.
It depends on how the data are arranged
##
x-matrix(c(1,2,3,2,8,2,4,5,6),nrow=3)
y-matrix(c(10,2,13,0,8,4,4.2,5.2,6.2),nrow=3)
q-mapply(t.test,as.data.frame(x),as.data.frame(y))
q
## The ith column of q contain the results of applying t.test
Hi, if i understood correctly, you probably want the something like:
## R
format(Sys.time(), LMSMOD_%b_%d_%y_%H%M)
(works on windows, don't know about Mac)
Hope that helps a little,
Tony Breyal
On 8 Jan, 23:47, greggal...@gmail.com wrote:
I have a function that is called from a sourced
what does R have to compare with , say , proc surveymeans, estimate survey
means/proportions with standard errors, using Taylor methods?
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
A procedure that , after adjusting for sampling weights, also
explicitly does an age adjustment to conform with an age distribution
of an older census?
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
Dear all,
The basename() function returns the extension also:
myfile - path1/path2/myoutput.txt
basename(myfile)
[1] myoutput.txt
Is there any other function where it just returns
plain base:
myoutput
i.e. without 'txt'
- Gundala Viswanath
Jakarta - Indonesia
It has worked wonders for me over the last years. It clunkily and
reliably delivers messages to my inbox that are either someone elses
question about R or an answer to one of my questions.
clunkily yours
Stephen Sefick
On Thu, Jan 8, 2009 at 8:58 PM, Robert Wilkins irishhac...@gmail.com wrote:
You can use 'sub' to get rid of the extensions:
sub(^([^.]*).*, \\1, 'filename.extension')
[1] filename
sub(^([^.]*).*, \\1, 'filename.extension.and.more')
[1] filename
sub(^([^.]*).*, \\1, 'filename without extension')
[1] filename without extension
On Thu, Jan 8, 2009 at 9:10 PM, Gundala
I have been corresponding with Harald Baayen about this problem, and it
appears that something is broken with the mcmcsamp function that pvals.fnc
calls. Running pvals.fnc on an lmer object actually changes the contents of
the object (get a summary before and after running pvals.fnc and you'll
On 1/8/09, Robert Wilkins iwriteco...@gmail.com wrote:
what does R have to compare with , say , proc surveymeans, estimate survey
means/proportions with standard errors, using Taylor methods?
survey package by Thomas Lumley has pretty much everything you would
need. Way more than SAS, at any
I am new to R but have spent a little time messing around with it. I
am running version 2.8.1 under Leopard OS 10.5 on a power PC and also
on my laptop, which is Intel. Unless, I am wrong, this version is just
a bit buggy.
1. In the R Commander, there is no direct entry of data. The Edit
Mostly I am interested in using R for statistics.
I am also interested in being able to look at source code.
I hope to be able to write extensions.
I tried the suggestion below but was unable to access vcov.lm.
methods(vcov)
[1] vcov.Arima* vcov.glm* vcov.lm*vcov.mlm* vcov.nls*
stats::vcov.lm
Error: 'vcov.lm' is not an exported object from 'namespace:stats'
stats:::vcov.lm
function (object, ...)
{
so - summary.lm(object, corr = FALSE)
so$sigma^2 * so$cov.unscaled
}
environment: namespace:stats
/Henrik
On Thu, Jan 8, 2009 at 1:44 AM, Carlos J. Gil Bellosta
library(StreamMetabolism)
day - (structure(c(10.08, 10.08, 10.05, 10.03, 10, 9.98, 9.96, 9.95,
9.95, 9.96, 9.96, 9.98, 10.01, 10.05, 10.06, 10.09, 10.11, 10.11,
10.13, 10.13, 10.15, 10.15, 10.13, 10.14, 10.11, 10.13, 10.14,
10.13, 10.12, 10.13, 10.14, 10.16, 10.18, 10.19, 10.23, 10.27,
10.33,
1 - 100 of 106 matches
Mail list logo