similar to: R equivalent to `estimate' in SAS proc mixed

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,]