Hello,
I was under the impression that in a call to gls, if i specify
gls(.., cor=AR1(form= ~ 1 | subject))
and suppose there are N subjects, then i would have obtained estimates
for N Phi (s) (the AR1 coefficient) - but I only get 1 estimate.
I then checked "Linear Mixed-Effects Models Using R
telemStats() should return to the prompt and the
print in telemStats should not appear
On Mon, Mar 16, 2015 at 4:02 PM, David Winsemius wrote:
>
> On Mar 16, 2015, at 3:08 PM, Saptarshi Guha wrote:
>
>> Hello,
>>
>> I would like a function X to return to the place that c
Hello,
I would like a function X to return to the place that called the
function XParent that called this function X.
Y calls XParent
Y = function(){
XParent()
print("hello")
}
XParent calls X
XParent = function(){
X()
print("H")
}
X returns to the point just after the call to XParen
Hello,
Couldn't think of a better subject line. Rather than a matrix like
x,y
..,..
.,..
I have a matrix like
x,y,n,
..,..,..,
..,..,..
and so on. Also, sum(n) is roughly few hundred million. The number of rows
is <1MM
Are they routines to fit a 2d kde estimate to data provided in this form?
I
hello,
This
rbindlist(list(list(a=NA,b=NA),list(a=20,b=FALSE)))
returns
a b
1: NANA
2: TRUE FALSE
as per the documentation ?rbindlist
is there a way to specify the column class of 'a' to be numeric?
In actual usage, i wont be able to re-order the 2nd list entry to be the
firs
't have to write two statements (a sequential
> and a parallel statement) -- if a parallel backend is running, it will use
> that, otherwise it'll execute in sequential mode.
>
> --j
>
> From: r-help-boun...@r-project.org [r-help-boun
Hello,
I find Map to be nice interface to mapply. However
Map calls mapply which in turn calls mapply via .Internal.
Is there a parallel version of mapply (like mcapply) or do I need to write
this myself?
Regards
Saptarshi
[[alternative HTML version deleted]]
___
Hello,
Consider the following transcript
> x=NULL
> x$date=10
> x
$date
[1] 10
or
> x=NULL
> x[10]=10
> x
[1] NA NA NA NA NA NA NA NA NA 10
Wouldn't one expect that x remain NULL despite the further additions (i.e
the x$date, x[10]) etc? Is coercing appropriate here?
Cheers
Saptarshi
Hello,
I was looking at the R (installed on RHEL6) shell script and saw
R_HOME_DIR=/usr/lib64/R. Nowhere (and I could have got it wrong) does
it read in the environment value R_HOME_DIR. I have the need to rsync
the entire folder below /usr/lib64/R to another computer into another
directory locati
Hello,
If i were to embed R (in C), i could over ride
ptr_R_ShowMessage = Re_ShowMessage;
ptr_R_WriteConsoleEx =Re_WriteConsoleEx;
where Re_ShowMessage,Re_WriteConsoleEx are my own functions.
How would i do this in a package?
I would like the flexibility of having my function be called /every-t
quot; "base"
> boot "boot" "recommended"
> ...
> mgcv "mgcv" "recommended"
> ...
> tools "tools" "base"
> utils "utils" "base"
>
> Bi
Hello,
Is there an R function that tells me the packages that come with R e.g.
c("base","boot","methods","mgcv",...)
Thank you
Saptarshi
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/lis
thanks much.
Saptarshi
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Hello,
How would i use bquote to create names for a list.
e.g. to create list(a=10)
bquote(list(.(X)=10), list(X="a"))
does not work.
The best i could come up with is
bquote({ a=list(10);names(a)=.(X); a}, list(X="a"))
which is quite ugly.
Is there an elegant way to solve this?
Regards
Sa
Hello,
[ Sorry, I sent the last email as HTML, this time it's in plain text ]
Suppose my variable,S, (time for something to start) is a skew-normal
distribution [1]. Can glm and mgcv handle this type of distribution for the
dependent variable?
Regards
Saptarshi
[1] http://azzalini.stat.unipd.it/S
I'm trying to install a package on OS X (Snow Leopard) using the following
command, but instead of respecting the arch directive it install i386. How
can this be resolved?
R --arch=x86_64 CMD INSTALL
* installing to library
â/Library/Frameworks/R.framework/Versions/2.15/Resources/library
Hello,
(Sorry for the repost, i am resending in plain text)
I have a function 'ewrap' (see below for definition).
It takes 3 expressions and returns another expression e.g.
map <- ewrap({
len <- length(r$addon)
rhcollect(len,1)
})
becomes:
expression({
NULL
result <- mapply(fu
Hello,
I have a function 'ewrap' (see below for definition).
It takes 3 expressions and returns another expression e.g.
map <- ewrap({
len <- length(r$addon)
rhcollect(len,1)
})
becomes:
expression({
NULL
result <- mapply(function(.index, k, r) {
{
len <- l
Hello,
Is there a way to add ylab on the right hand side also (in lattice)?
Different from the left hand side?
Cheers
Saptarshi
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-h
Hello,
I am lattice scatterplot that has 2 panels (could be a few more). Both
panels have a y-axis label that is different on the
left and right sides. However the right hand side axis labelling depends on
which panel i'm in .
(I am plotting two curves in one panel, the left y-axis has the scale
Hello,
I have the following code
asTheEconomist(bwplot(ver~quans,data=z,cex=0.5,subset=abs(quans)<1.1,scales=list(x=list(tick.num=20,cex=0.55,alternating=3))
,xscale.components = xscale.components.CF
,main=list('FOOO"',cex=0.8)))
The xscale.components displays the exponent of the
Hello,
I would like to write a function where substitute operates on the
parameter, but ...
> Expression = function(o,l) substitute(o, l)
> Expression({x=.(FOO)}, list(FOO=2))
o
How do i get substitute to work on the contents of o.
Regards
Saptarshi
_
Hello,
I'm trying to use rjson to parse a JSON object. The object can be found here
http://pastebin.com/np0s5hgM
(you'll probably need to add quotes around the content)
fromJSON returns the error
glibc detected home/sguha/lib64/R/bin/exec/R: realloc(): invalid next
size: 0x0209e210 **
Hello,
I notice that e.g /home/sguha/lib64 is hard coded into the /bin/R file .
I nstalled R as ./configure --prefix=$HOME ...
What i need to do is ship the entire R distribution to remote nodes,
and run R. These are shipped to ephemeral directories
so I dont know the path ahead of time.
R_HOME
Hello,
Given the following display
library(maps)
library(mapproj)
m <- map('world',plot=FALSE)
map('world',proj='mollweide',bg=bgcolor[2],col='white')
map.grid(col=2,lim=c(-175,-175,-180,180),label=FALSE,lty=2)
Is there a way to color the region outside the dotted red border?
Thank you
Saptars
Hello,
I should see how to mirror a google group (if indeed possible).
There are few things stopping from having RHIPE work on the remote
while the client sends queries from a Windows client.
But as of now, it does not work.
Regards
Saptarshi
On Fri, Mar 18, 2011 at 7:50 AM, XiaoboGu wrote:
(\\d+)\\s\"(.*?)\"\\s\"(.*?)\"\\s\"(.*?)\"$"
> > regexpr(pat, test, perl=TRUE)
> [1] 1
> attr(,"match.length")
> [1] 436
>
> 3. Consider a different approach, e.g. scan(textConnection(test),
> what=character(0))
>
> Hope this
Hello R users,
I have this regex see [1] for apache log lines. I tried using R to parse
some data (only because I wanted to stay in R).
A sample line is [2]
(a) I saved the line in [1] into "~/tmp/a.txt" and [2] into "/tmp/a.txt"
pat <- readLines("~/tmp/a.txt")
test <- readLines("/tmp/a.txt")
te
Hello,
Consider this little function (not concerned about efficiency here)
m=function(r) { a<-serialize(r,NULL); zz<-rawConnection(raw(0),"r+");
writeBin(length(a),zz); c(rawConnectionValue(zz),a)}
a=m(10)
p=rawConnection(a,"r+"); u=readBin(p,"int")
unserialize(p)
Error in unserialize(p) : unk
Dirk, thanks.
On Thu, Oct 14, 2010 at 9:42 AM, Dirk Eddelbuettel wrote:
>
> On 14 October 2010 at 11:21, Dirk Eddelbuettel wrote:
> | On 14 October 2010 at 08:49, Saptarshi Guha wrote:
> | | Hello,
> | | I have this script which will be invoked as
> | |
> | | Rscript a.r
Yes, thanks much.
On Thu, Oct 14, 2010 at 9:36 AM, jim holtman wrote:
> ?chartr
>
> Is this what you are looking for?
>
> On Thu, Oct 14, 2010 at 12:10 PM, Saptarshi Guha
> wrote:
>> Hello,
>>
>> Would one use gsub to perform pythons translate?
>>
&
Hello,
Would one use gsub to perform pythons translate?
in python
table = string.maketrans("abc","xyz")
"abc".translate(table)
"xyz"
how does one do this in R?
Thank you
Saptarshi
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/lis
Hello,
I have this script which will be invoked as
Rscript a.r < a.r
a.r follows
#!/usr/bin/Rscript --vanilla
f=file("stdin")
while(TRUE){
y=readLines(f,n=1,warn=TRUE)
if(length(y)==0) break else print(y)
}
But it only reads one line from a.r
How can I read line by line from standard inp
Hello,
I take this back, wrong implementation of the experiment.
Sorry for the noise.
Regards
Saptarshi
On Sun, Aug 15, 2010 at 12:22 PM, Saptarshi Guha
wrote:
> Hello,
> I have question about the overhead in lapply.
> x is a list of 3000 lists. Each of the i (1<=i<=3000) l
Hello,
I have question about the overhead in lapply.
x is a list of 3000 lists. Each of the i (1<=i<=3000) list elements is
pair of two elements: a string vector and a data frame
x is roughly 235MB.
> gc()
##
> z <- system.time(y <- lapply(x,function(r){
system.time(serialize(r,NULL))['elapsed
Hello,
I have the following data frame, using lattice and qqmath, how do I display the
quantile plot (e.g distribution = qunif), but highlight the f==TRUE and f==FALSE
differently?
library(lattice)
pp <- function(n) ((1:n)-0.5)/n
x=data.frame(x=runif(20),f=sapply(runif(20),function(r) r>0.5))
I
Hello,
Do you need /all/ the data in memory at one time? Is your goal to
divide the data (e.g according to some factor /or/ some function of
the columns of data set ) and then analyze the divisions? And then,
possibly, combine the results ?
If so, you might consider using Rhipe. We have analyzed (e
Dear R Users,
Using codetools I obtained the text representation of the parse tree for this
snippet
z=quote({x[1]<-2})
showTree(z)
> ("{" (<- ([ x 1) 2)) (A)
If I understand correctly, x[1]<-2 ought to be "[<-"(x,1,2), so shouldn't i see
("{" ( [<- x 1 2 ) )
If indeed the parse tree in
Hello,
I would like to compose a series of functions each taking one
argument. Let us assume the values returned by the functions are
contained within the domains.
So
b<-function(r) r[-1]
compose("b","b")
b(b(x))
This is my solution, is there something succinct ?
Regards
Saptarshi
compose <- fu
ust want : exportPattern("^rh")
>
> What did you expect the backslash to do ?
>
> On 12/01/2009 05:18 PM, Saptarshi Guha wrote:
>>
>> Hello,
>> I wrote a package, which in the NAMESPACE file exports functions like
>> this:
>>
>> exportPattern(
As a followup, i can do things Package:::rh*, there just no exported.
On Tue, Dec 1, 2009 at 11:18 AM, Saptarshi Guha
wrote:
> Hello,
> I wrote a package, which in the NAMESPACE file exports functions like this:
>
> exportPattern("^\\rh")
>
> On R-2.8 (Linux, 64), upo
Hello,
I wrote a package, which in the NAMESPACE file exports functions like this:
exportPattern("^\\rh")
On R-2.8 (Linux, 64), upon loading the package I have the rh functions present.
On R-2.10, Mac OS X, (32 bit), it builds, loads, but the functions are not
loaded, i.e the only function is rh
> The substitute() function can do it, but the bquote() function is the most
> convenient:
>
>> foo <- function(s) {
> + bquote(expression({
> + v <- .(s)
> + print(v)
> + }))
> + }
>> foo(10)
> expression({
> v <- 10
> print(v)
> })
>
> Duncan Murdoch
>
Thanks, after sending this em
when eval'd e.g
x <- foo(10)
eval(x)
Error in eval(expr, envir, enclos) : object 's' not found
Q: How do I force a substitution so that the returned expression has the value
of 's' 'hardcoded' in it e.g
> foo(10)
returns
expression({
v <- 10
I have some beginner's questions regarding local, in the docs, it says that
"local evaluates an expression in a local environment".
Q1: why is B different from A? In B, is a<-a+1 getting evaluated
before eval proceeds?
#A
a=0
eval(quote(a<-a+1),new.env())
a # 0
#B
a=0
eval(a<-a+1,new.env())
a #
Hello,
(Previous message was scrubbed)
I apologize if this is not the list for package announcements. Having
said that I've released a package called rprotos which implements a
simple
serialization using Googles protocol buffers[1].
The package also includes some miscellaneous functions for
writin
it is possible to
extend this to serialize functions, expressions, environments and
several other objects.
However that is some time in the future.
Regards
Saptarshi
[1] http://code.google.com/apis/protocolbuffers/docs/overview.html
Download:
http://ml.stat.purdue.edu/rproto_1.0.tar.gz
Saptarshi
Hello,
I have used socketConnection to connect to a TCP server. I havent
figured out a way to do the same with a UDP server.
i.e I have a server listening on 9000, communicating via UDP. I would
like to , from R, send packets to this server,
This does not work
u <- socketConnection('localhost',9
assignments made by this operator are located.
>
>
> Bill Venables
> http://www.cmis.csiro.au/bill.venables/
>
>
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
> Behalf Of Saptarshi Guha
> Sent: Friday, 17 July 2009
Hello,
I tried this pseudo-generator style experiment (see below). The "<<-"
operator assigns to in the calling environment which would be the
environment of "getN".
Yet when the function incr is returned, isn't this environment lost?
Also the print displays GlobalEnv, yet the globalenv does not h
Hello,
On Fedora 8, yum install R , installs R-2.6.
I dont have much experience with yum, how do change the repos(?) to
install R-2.9?
Thank you
Saptarshi
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read t
Hello,
I am using EC2 to launch several instances. On each instance, I
perform the following commands
yum -y install R (installs R-2.6)
R CMD javareconf
sleep 10
wget rpackage.tgz
R CMD INSTALL rpackage.tgz
Now, rpackage.tgz needs to be built with jni libraries. When the
instance is fully up, i l
> You've hit on the reason -- representation of NA
>
> src/main/arithmetic.c:157: R_NaInt = INT_MIN;
Aha! Yes, I missed that.
>
> (R-devel seems more appropriate for this question).
Yes, i was not sure.
Thank you
Saptarshi
__
R-help@r-project.org
Hello,
I'm using R-2.9 and have a question about integer ranges. On C, (OS X
platform),
LONG_MIN and LONG_MAX(from limits.h) correspond to -2147483648(-2^31)
and 2147483647 (2^31-1)
In Java(1.4), Integer.MIN_VALUE and Integer.MAX_VALUE have the same values.
On R, as.integer(2^31-1) is the maximum
educe=r,input.folder="X",output.folder="Y")
URL: http://ml.stat.purdue.edu/rhipe
There are some downsides to RHIPE which are described at
http://ml.stat.purdue.edu/rhipe/install.html#sec-5
Regards
Saptarshi Guha
__
R-help@r-project.
fortunately, I
haven't
had the chance to compare with rmpi and snow.
TODO(will appear soon ): another function: rhmr - mapreduce using R.
Regards
Saptarshi Guha
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
http
Hello R users
I have question about the time involved in list assignment.
Consider the following code snippet(see below). The first line creates
a reader object,
which is the interface to 1MM key-value pairs (serialized R objects) spanning 50
files (a total of 50MB). rhsqstart initiates the reading
.jclassPath()
does contain an entry to this class.
Moreover, when run at the command line using the *same* classpath, I
don't get any such error messages.
Is there a possible reason why?
Regards
Saptarshi
Saptarshi Guha
__
R-help@r-project.org
quite understand it. What does the
option do and why did I initially run into an error ?
>From man R
--max-nsize
Set max number of cons cells to N
Why should there be a maximum?
Thanks for your time
Saptarshi Guha
__
R-help@r-project.org mail
)
The two figures are not base aligned. I would like them share the same
the baseline and same height, if necessary the paper width and height
can be adjusted
( i tried setting the paper width and height to no avail).
Is the way to base align the two figur
more information and examples please visit this page: http://www.stat.purdue.edu/~sguha/rhipe
.
If anyone would like to contribute to this project, please email me
directly - any help is welcome.
Regards
Saptarshi Guha
__
R-help@r-project.org ma
to this instead of iterating over all objects E
and assigning into .GlobalEnv?
Thank you
Saptarshi
--
Saptarshi Guha - saptarshi.g...@gmail.com
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
Bin(length(x),con=conn,endian="big") #conn opened in binary
But the response on the server side is wrong, eg. 1 is received as
4294967297
How do I then write long across the network from R?
Thank you
Saptarshi
Saptarshi Guha | saptarshi.g
und this?
Regards
Saptarshi
R.version = 2.8
>
> On Nov 30, 2008, at 5:53 PM, Saptarshi Guha wrote:
>
>> Hello,
>> I have some code which generates lattice objects. The function
>> recieves serialized forms of the lattice objects which it
o interpret
[EMAIL PROTECTED]
*** REXPFactory unable to interpret
[EMAIL PROTECTED] *
** REXPFactory unable to interpret
[EMAIL PROTECTED]
How can I in general,force the unserialize to return a list?
Saptarshi Guha | [EMAIL PROTECTED] | http://www.stat.purdue.edu/~sguha
Marijuana
dataset, the y-axis
labels overlap so I would like to split it across several panels(1
panel per page) with about 25-50 per panel.
Is there a way to do this?
Regards
Saptarshi
Saptarshi Guha | [EMAIL PROTECTED] | http://www.stat.purdue.edu/~sguha
[[alternative HTML
e of the dummy variable pulling in data from a postgres database
setting my limits for each panel based on the pulled in data and
plotting that.
Saptarshi Guha | [EMAIL PROTECTED] | http://www.stat.purdue.edu/~sguha
[[alternative HTML version deleted]]
__
tles.)
Is there a way to solve this? Also, each level of C appears in exactly
one
combination of A and B (the combinations that appear in the data set,
not the ones produced by lattice)
Kind regards Saptarshi
Saptarshi Guha | [EMAIL PROTECTED] | http://www.stat.purdue.edu/~sguha
[[a
Forgot to send to the list.
Begin forwarded message:
> From: Saptarshi Guha <[EMAIL PROTECTED]>
> Date: September 19, 2008 1:43:50 PM GMT-04:00
> To: Gabor Grothendieck <[EMAIL PROTECTED]>
> Subject: Re: [R] Lines between panels in lattice
>
> Certainly.
>
&g
Got it. It should be plot1.panel.1.2.off.vp and likewise.
Thanks
Saptarshi
On Sep 19, 2008, at 12:08 PM, Saptarshi Guha wrote:
> Hello,
> I have a multi-page display each consisting of two-panels above each
> other.
> I need to draw a line from the top panel to bottom
lot1.panel.1.1.vp")
grid.line.to(unit(1.2,"native"),unit(3.5,"native"))
However, I can see the line in the bottom panel, but nothing in the
top. Is there some
clipping going on?
Thanks
Saptarshi
Saptarshi Guha | [EMAIL PROTEC
ght happen (as pointed out in the help page)
but how can I get around it? Should I specify a width to the device?
Thank you
Saptarshi
Saptarshi Guha | [EMAIL PROTECTED] | http://www.stat.purdue.edu/~sguha
[[alternative HTML version deleted]]
__
(Mac) device, there is very little padding
between the bottom of the curve and the lower horizontal axis and
between the maximum values of the curve and the upper horizontal axis.
How can i increase the padding - top,bottom,left and right?
Thanks in advance for your time.
Hello,
Is there a way to add minor tick marks to the Y-axis of a lattice
plots?
Thank you
Saptarshi
Saptarshi Guha | [EMAIL PROTECTED] | http://www.stat.purdue.edu/~sguha
[[alternative HTML version deleted]]
__
R
dding between the endpoints and
>> the
>> vertical axes.
>>How can I remove this padding?
>
> A reproducible example would have helped, but perhaps you want
>
> scales = list(x = list(relation = "free", axs = "i"))
>
>
Hello,
I need to plot a line graph using lattice so that the endpoints are
flush with the left and right vertical axes.
So the scale end points and horizontal axes length should coincide.
This works
u=newdata[1:100,] #u[,1] is from 1700-1795
xyplot
Thank you.
On Feb 17, 2008, at 1:18 PM, Deepayan Sarkar wrote:
> See ?print.trellis.
>
> p <- xyplot(1 ~ 1, aspect = 0.5)
> p
> plot(p, panel.width = list(2, "inches"), panel.height = list(1,
> "inches"))
>
> This overrides 'aspect=' though, that is, you will need to specify
> both width and heig
same scales.
I suppose i need to do something in prepanel but what?
Is there any function i can call to set the actual length of the axes?
Thank you
Saptarsho
Saptarshi Guha | [EMAIL PROTECTED] | http://www.stat.purdue.edu/~sguha
[[alternative HTML ve
thing in prepanel but what?
Is there any function i can call to set the actual length of the axes?
Saptarshi Guha | [EMAIL PROTECTED] | http://www.stat.purdue.edu/~sguha
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/lis
Hello,
I have an optimization problem which I would like to solve in R but
am not sure how to approach this.
I have a cost function C which is quadratic in rx and ry (if it
helps, the cost function is separable in rx and ry) and is always
positive for all rx,ry.
I wish
Hello,
I downloaded R-2.5.1 and configured it on a AMD x86-74 running
Redhat. During compilation i get this error
configure --prefix=$HOME/mine
make
mkdir -p -- ../../../../library/methods/libs
make[5]: Leaving directory `/home/sguha/tmp/R-2.5.1/src/library/
methods/src'
make[4]: Leaving direc
81 matches
Mail list logo