Displaying 20 results from an estimated 5000 matches similar to: "R equivalent to `estimate' in SAS proc mixed"
2006 May 19
1
make error for R-2.3.0
Hello,
I'm trying to install R on a linux machine running Red Hat 8. I ran
./configure
make
and get the following error. I've installed several versions of R (2.2.1
most recently) on this machine and haven't had any problems until now. I
wondered if the outdated compiler (gcc version 3.2) was the problem and
attempted to install my own, more recent version. I tried gcc versions 4.1.0
2002 Aug 29
3
tru64 patch: openssh-SNAP-20020826.tar.gz does not contain 'configure', so how to build?
Hi-
Since the tru64 patch was designed for -current, I thought I would try
to build it with a recent snapshot before backporting to 3.4p1.
So I downloaded openssh-SNAP-20020826.tar.gz frpm the portable snapshots, but it
does not contain the 'configure' script. I tried copying the 'configure' from
3.4p1, but that does not create a Makefile from the Makefile.in.
Where are the
2002 Sep 11
1
tru64 sia: move call of session_setup_sia() to do_setusercontext(), letting grantpty() and friends handle pty perms
Hi-
Under privsep, I experimented with moving the session_setup_sia() out of
do_child() and into do_setusercontext(), which is where the uids/gids are set
to the final execution user. The call is made with a NULL tty, and this
is functional provided that any later pty allocation uses grantpty() to
set the device permissions. Logging in with this method shows that a utmp
entry does get made for
2010 Jul 29
1
R Equivalent of SAS Datastep Line-Hold (@@) Specifier?
Hello Everyone,
Below is some SAS code that uses a "line hold specifier" to read multiple observations from each of several input lines of data. There are 3 patients per line in the in-stream data.
Is there a way in R to read this kind of data? I've looked in my books and online but haven't found anything
Thanks,
Paul
DATA EXAMP.TRIAL;
INPUT TRT $ CENTER PAT SEX $ AGE
2010 Jul 25
2
R equivalent of SAS proc freq
Dear R-users,
I am looking for a R function that would be the equivalent of the SAS proc
freq (
http://support.sas.com/documentation/cdl/en/procstat/63104/HTML/default/viewer.htm#/documentation/cdl/en/procstat/63104/HTML/default/procstat_freq_sect006.htm).
The table, ftable, xtabs functions are close but do not quite offer the same
capabilities (e.g. they just return counts and no %ages as far as
2011 Aug 08
1
mixed model fitting between R and SAS
Hi al,
I have a dataset (see attached), which basically involves 4 treatments for a chemotherapy drug. Samples were taken from 2 biopsy locations, and biopsy were taken at 2 time points. So each subject has 4 data points (from 2 biopsy locations and 2 time points). The objective is to study treatment difference.?
I used lme to fit a mixed model that uses "biopsy.site nested within pid"
2017 Oct 11
0
Converting SAS Code
I have no problem setting up my mixed model, or performing anova or lsmeans
on my model?s outputs. However, performing lsd mean separation is giving me
fits.
So I do not have a problem when using two-way anova model. When using the
code:
fit.yield.add <- lm(data = ryzup, Yield ~ Rep + Nitrogen + Treatment)
LSD.test(fit.yield.add, trt = "Nitrogen", alpha = 0.1, console = TRUE)
2002 Sep 04
2
uid transition and post-auth privsep (WAS Re: possible fundamental problem with tru64 patch) (fwd)
What do we loose by not having post-auth privsep?
What code is executed between authorization and actual setting of the
effective uid?
On Tue, 3 Sep 2002, Chris Adams wrote:
> Once upon a time, Toni L. Harbaugh-Blackford <harbaugh at nciaxp.ncifcrf.gov> said:
> > It appears that the integration of the sia session setup will either
> > have to be rethought or abandoned
2007 Jul 14
1
return() in nested functions
Dear WizaRds,
After consulting different sources I am still unable to understand the
correct use of return() in nested functions. To illustrate the problem:
f <- function(x,y,type){
est1<-function(x,y){
z=x+y
out(x,y,z)}
est2<-function(x,y){
z=x*y
out(x,y,z)}
out<-function(x,y,z)
return(x,y,z)
if (type=="est1") est1(x,y)
if (type=="est2") est2(x,y)
}
2004 Feb 08
1
a newbie question.
Hi,
I'm a newbie.Sorry for my question.
I have RH9 & downloaded rpm of dovecot for RH9.
I installed it:
rpm -ivh nameOfthatrpm
& when i ran
#service dovecot start ,I recieved following lines:
Starting IMAP daemon (dovecot): Fatal: Can't use SSL
certificate /usr/share/ssl/certs/dovecot.pem: No such
file or directory
2006 Apr 26
1
help using tapply
Dear R-mates,
# Here's what I am trying to do. I have a dataset like this:
id = c(rep(1,8), rep(2,8))
dur1 <- c( 17,18,19,18,24,19,24,24 )
est1 <- c( rep(1,5), rep(2,3) )
dur2 <- c(1,1,3,4,8,12,13,14)
est2 <- rep(1,8)
mydata = data.frame(id,
estat=c(est1, est2),
durat=c(dur1, dur2))
# I want to one have this:
id = c(rep(1,8), rep(2,8))
2004 Apr 19
2
Can't resolve address 127.0.0.1
Hi ,
I recently installed Fedora Core 1.
Installed postfix on it and that is running fine.
I am trying to get dovecot to run and I have changed only the following line
from the default dovecot.conf :
imap_listen = 127.0.0.1 since I'd like to install squirrelmail later on.
Now when I try to start dovecot I get "Starting Dovecot Imap: Fatal: Can't
resolve address
2007 Feb 27
0
Optimizing the loop for large data
Rusers:
I am trying to apply a quadratic discriminant function to find the best
classification outcomes.
1 is assigned to the values greater than a threshold value; and 0 otherwise.
I would like to see how the apparent error rates and the optimal error
rate change with increasing threshold values.
I have a 1000*10 data matrix: n=1000 and p=10.
Here is what I wrote so far, but seems to be
2004 Feb 10
4
UW-IMAP -> Dovecot : patch to remove nasty "internal messages" from UW-IMAP ?
Hello there,
We want to move from UW-IMAP to Dovecot (mbox format), works like a
charm, but UW-IMAP has that bloody habbit to make a stupid message in
all his mailbox files :
From MAILER-DAEMON Thu Apr 19 10:30:01 2001
Date: 19 Apr 2001 10:30:01 +0200
From: Mail System Internal Data <MAILER-DAEMON at mailhub2-vil.isdnet.net>
Subject: DON'T DELETE THIS MESSAGE -- FOLDER INTERNAL DATA
2002 Aug 28
0
patch almost works on 5.1A openssh 3.4p1 - get in, but get kicked out (fwd)
Hi-
I applied the privsep patch to Tru64 5.1A openssh 3.4p1 and it
*almost* works.
I get in from the client side and xauth is run, but in the meantime
the server side disconnects. Running sshd in debug mode level 3 gives
the following output:
.
.
.
debug1: session_input_channel_req: session 0 req shell
debug1: fd 5 setting TCP_NODELAY
debug1: channel 0: rfd 13
2007 Feb 12
1
Width of a plotting point (in inches) in grid package
Hello,
I'm trying to determine the width of a plotting point (in inches) in the
grid package. I naively thought I could create a pointsGrob with only one
point and get the width (as tried below), but this results in an object with
a size of 0inches (changing cex has no effect). Does anyone have a better
approach? Of course, it would be dependent upon the graphics parameters and
viewport...
2008 Oct 15
2
Network meta-analysis, varConstPower in nlme
Dear Thomas Lumley, and R-help list members,
I have read your article "Network meta-analysis for indirect treatment
comparisons" (Statist Med, 2002) with great interest. I found it very
helpful that you included the R code to replicate your analysis;
however, I have had a problem replicating your example and wondered if
you are able to give me a hint. When I use the code from the
2006 Aug 24
1
how to constrast with factorial experiment
Hello, R users,
I have two factors (treat, section) anova design experiment where
there are 3 replicates. The objective of the experiment is to test if
there is significant difference of yield between top (section 9 to 11)
and bottom (section 9 to 11) of the fruit tree under treatment. I
found that there are interaction between two factors. I wonder if I
can contrast means from levels of
2006 Aug 29
0
how to contrast with factorial experiment
Hello, R experts,
If I understand Ted's anwser correctly, then I can not contrast the
mean yields between sections 1-8 and 9-11 under "Trt" but I can
contrast mean yields for sections 1-3 and 6-11 because there exists
significant interaction between two factors (Trt:section4,
Trt:section5). Could I use the commands below to test
the difference between sections 1-3 and 6-11 ?
2007 Mar 14
1
How to transform matrices to ANOVA input datasets?
Hello, R experts,
I have a list called dataHP which has 30 elements (m1, m2, ..., m30).
Each element is a 7x6 matrix holding yield data from two factors
experimental design, with treatment in column, position in row. For
instance, the element 20 is:
dataHP[[20]]
col1 col2 col3 trt1 trt2 trt3
[1,] 22.0 20.3 29.7 63.3 78.5 76.4
[2,]