[R] Using Sampling Weights in R

2007-04-11 Thread Thomas W. Volscho
Dear List, I have a dataset that provides sampling weights (National Survey of Family Growth 2002). I want to produce a cross-tabulation and use the provided sampling weights to obtain representative population estimates. (I believe they are simply frequency weights but codebook is

Re: [R] Reasons to Use R

2007-04-11 Thread Wensui Liu
Greg, As far as I understand, SAS is more efficient handling large data probably than S+/R. Do you have any idea why? On 4/10/07, Greg Snow [EMAIL PROTECTED] wrote: -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Bi-Info

Re: [R] Reasons to Use R

2007-04-11 Thread Greg Snow
I think SAS has the database part built into it. I have heard 2nd hand of new statisticians going to work for a company and asking if they have SAS, the reply is Yes we use SAS for our database, does it do statistics also? Also I heard something about SAS is no longer considered an acronym, they

Re: [R] read.spss (package foreign) and SPSS 15.0 files

2007-04-11 Thread John Kane
--- Thomas Lumley [EMAIL PROTECTED] wrote: On Thu, 5 Apr 2007, John Kane wrote: Heck. I'd be happy to get an answer to what is happening here: mac - spss.get(H:/ONTH/Raw.data/Follow.sav) Warning message: H:/ONTH/Raw.data/Follow.sav: Unrecognized record type 7, subtype 16

Re: [R] read.spss (package foreign) and SPSS 15.0 files

2007-04-11 Thread John Kane
--- Prof Brian Ripley [EMAIL PROTECTED] wrote: On Thu, 5 Apr 2007, Michael Conklin wrote: Not being the developer I cannot answer definitively but, as a frequent user of SPSS files I can give you my experience. 1) The unrecognized coding is perhaps due to the locale of the SPSS

Re: [R] read.spss (package foreign) and SPSS 15.0 files

2007-04-11 Thread Thomas Lumley
On Thu, 5 Apr 2007, John Kane wrote: Heck. I'd be happy to get an answer to what is happening here: mac - spss.get(H:/ONTH/Raw.data/Follow.sav) Warning message: H:/ONTH/Raw.data/Follow.sav: Unrecognized record type 7, subtype 16 encountered in system file It means that your file had a

Re: [R] plot log scale, axis original scale

2007-04-11 Thread Petr PIKAL
Hi Dean Sonneborn [EMAIL PROTECTED] napsal dne 10.04.2007 18:28:43: Petr, This is great! Thank you so much for responding. Could I get one more point clarified. My A values range from 1 to 35. I would really like to use something like AT=1 to 35 by 5 instead of AT=log(a). at=log(a)

Re: [R] read.spss (package foreign) and SPSS 15.0 files

2007-04-11 Thread Prof Brian Ripley
On Thu, 5 Apr 2007, Michael Conklin wrote: Not being the developer I cannot answer definitively but, as a frequent user of SPSS files I can give you my experience. 1) The unrecognized coding is perhaps due to the locale of the SPSS installation. I have had success reading in files from

Re: [R] read.spss (package foreign) and SPSS 15.0 files

2007-04-11 Thread Prof Brian Ripley
On Thu, 5 Apr 2007, Michael Conklin wrote: Not being the developer I cannot answer definitively but, as a frequent user of SPSS files I can give you my experience. 1) The unrecognized coding is perhaps due to the locale of the SPSS installation. I have had success reading in files from

Re: [R] read.spss (package foreign) and SPSS 15.0 files

2007-04-11 Thread John Kane
--- Prof Brian Ripley [EMAIL PROTECTED] wrote: On Thu, 5 Apr 2007, Michael Conklin wrote: Not being the developer I cannot answer definitively but, as a frequent user of SPSS files I can give you my experience. 1) The unrecognized coding is perhaps due to the locale of the SPSS

Re: [R] read.spss (package foreign) and SPSS 15.0 files

2007-04-11 Thread John Kane
--- Thomas Lumley [EMAIL PROTECTED] wrote: On Thu, 5 Apr 2007, John Kane wrote: Heck. I'd be happy to get an answer to what is happening here: mac - spss.get(H:/ONTH/Raw.data/Follow.sav) Warning message: H:/ONTH/Raw.data/Follow.sav: Unrecognized record type 7, subtype 16

Re: [R] read.spss (package foreign) and SPSS 15.0 files

2007-04-11 Thread Thomas Lumley
On Thu, 5 Apr 2007, John Kane wrote: Heck. I'd be happy to get an answer to what is happening here: mac - spss.get(H:/ONTH/Raw.data/Follow.sav) Warning message: H:/ONTH/Raw.data/Follow.sav: Unrecognized record type 7, subtype 16 encountered in system file It means that your file had a

[R] hypothesis testing

2007-04-11 Thread genomenet
Hi there, Can anyone tell me how to test a null hypothesis such as: h0:u1=u2=u3=u4 AND u5u6 means 'not equal' Thank you very much. Fan __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] Rserve and R to R communication

2007-04-11 Thread Ramon Diaz-Uriarte
On 4/11/07, AJ Rossini [EMAIL PROTECTED] wrote: On Tuesday 10 April 2007 23:17, Ramon Diaz-Uriarte wrote: Of course, you are right there. I think that might still be the case. At the time we made our decision, and decided to go for MPI, MPI 2 was already out, and MPI seemed more like the

[R] Boxplot names format

2007-04-11 Thread Jose Sierra
Hello. I need help. I create a boxplot but the names are too longs and i can´t see them complete. How can i format its size and font? How can i resize the boxplot? Thanks and regards. Ramon Diaz-Uriarte escribió: On 4/11/07, AJ Rossini [EMAIL PROTECTED] wrote: On Tuesday 10 April 2007

Re: [R] sas.get problem

2007-04-11 Thread Charles Dupont
John Kane wrote: I have 3 SAS files all in the directory F:/sas, two data files and a format file : form.ea1.sas7bdat form.ea2.sas7bdat sas.fmts.sas7bdat F is a USB. I am trying import them to R using sas.get. I have not used SAS since I was downloading data from mainframe and

Re: [R] Boxplot names format

2007-04-11 Thread Peter Danenberg
I create a boxplot but the names are too longs and i cant see them complete. If you're referring to labels on the x-axis, Jose, I'll sometimes rotate them and increase the bottom margin: # Increase bottom margin by 0.1 par(omd=c(0,1,0.1,1)) boxplot(...) # Add

Re: [R] Boxplot names format

2007-04-11 Thread Jose Sierra
Thank you very much Peters. it runs Peter Danenberg escribió: I create a boxplot but the names are too longs and i cant see them complete. If you're referring to labels on the x-axis, Jose, I'll sometimes rotate them and increase the bottom margin: # Increase

Re: [R] Computing fundamental harmonics from a periodogram

2007-04-11 Thread Uwe Ligges
hadley wickham wrote: Dear all, I'm trying to finding the fundamental harmonics (ie. peaks in a periodogram) from a time series (extracted from an mp3). For example, if I look at spectrum(fdeaths, spans = c(3,3)) A heuristical procedure for finding fundamentals of a monophonic sound

[R] enable-R-shlib ?

2007-04-11 Thread Bernardo Rangel Tura
Hi R Masters! Recently I migrated to Ubuntu Linux [I a former windows user]. Now I think compile a 64-bit R version for my computer [Turion AMD], but I not sure if using the configure option --enable-R-shlib. I note this option is usefull for some GUI like gnomeGui and JGR, but this will go

[R] Random Forest Imputations

2007-04-11 Thread Ricky Jacob
Dear All, I am not able to run the random forest with my dataset.. X- 280 records with satellite data(28 columns) - B1min, b1max, b1std etc.. y- 280 records with 3 columns - TotBasal Area, Stem density and Volume yref - y[1:230,] #Keeping 1st 230 records as reference records want to set 0

Re: [R] enable-R-shlib ?

2007-04-11 Thread Prof Brian Ripley
You'll find the definitive discussion in the R-admin manual. On Wed, 11 Apr 2007, Bernardo Rangel Tura wrote: Hi R Masters! Recently I migrated to Ubuntu Linux [I a former windows user]. Now I think compile a 64-bit R version for my computer [Turion AMD], but I not sure if using the

Re: [R] p-values and significance

2007-04-11 Thread Bernd Dittmann
Hi Paul, here's a lm model to illustrate this: summary(lm(y~x.1+x.2)) Call: lm(formula = y ~ x.1 + x.2) Residuals: Min 1Q Median 3QMax -0.0561359 -0.0054020 0.0004553 0.0056516 0.0515817 Coefficients: Estimate Std. Error t value Pr(|t|)

Re: [R] Random Forest Imputations [Broadcast]

2007-04-11 Thread Liaw, Andy
Please provide the information the posting guide asks (version of R, packages used, version of package used, etc). There are no yaImpute() or yai() functions in the randomForest package. Andy From: [EMAIL PROTECTED] on behalf of Ricky Jacob Sent: Wed

Re: [R] Random Forest Imputations [Broadcast]

2007-04-11 Thread Ricky Jacob
I am currently using R 2.4.1 version. Am using the yaImpute package for k-NN imputation.. http://forest.moscowfsl.wsu.edu/gems/yaImpute.pdf In yaImpute, i am using the yai function which uses randomForest as a method for finding out the k-Nearest Neighbours..

[R] ASYMMETRIC DCC MGARCH

2007-04-11 Thread James Raoul
Hello, Does any one know any implementation of the Asymmetric (in volatility) DCC MGARCH. I used the toolbox of sheppard but I have some problems in generalizing the GARCH-DCC-MGARCH to the Asymmetric one (EGARCH -DCC-MGARCH, GJR- DCC-MGARCH) Would you please correct to me the code that I

[R] Boxplot with quartiles generated from different algorithms

2007-04-11 Thread Pietrzykowski, Matthew \(GE, Research\)
R users: I am trying to replicate the boxplot output I achieve with Minitab in R. I realize that R gives the user many more options on the algorithm used to calculate the IQR than Minitab, so I concentrated on type=6 when using the quantile() function in R. The problem I am having is setting the

Re: [R] Memory management

2007-04-11 Thread yoooooo
I guess I have more reading to do Are there any website that I can read up on memory management, or specifically what happen when we 'pass in' variables, which strategy is better at which situation? Thanks~ - y Prof Brian Ripley wrote: On Tue, 10 Apr 2007, yoo wrote: Hi

[R] negative variances

2007-04-11 Thread Tu Yu-Kang
Dear R experts, I had a question which may not be directly relevant to R but I will be grateful if you can give me some advices. I ran a two-level multilevel model for data with repeated measurements over time, i.e. level-1 the repeated measures and level-2 subjects. I could not get

Re: [R] Boxplot names format

2007-04-11 Thread Charilaos Skiadas
On Apr 11, 2007, at 5:31 AM, Jose Sierra wrote: Thank you very much Peters. it runs Peter Danenberg escribió: I create a boxplot but the names are too longs and i cant see them complete. If you're referring to labels on the x-axis, Jose, I'll sometimes rotate them and

Re: [R] Reasons to Use R

2007-04-11 Thread Charilaos Skiadas
A new fortune candidate perhaps? On Apr 10, 2007, at 6:27 PM, Greg Snow wrote: Remember, everything is better than everything else given the right comparison. -- Gregory (Greg) L. Snow Ph.D. Haris Skiadas Department of Mathematics and Computer Science Hanover College

Re: [R] negative variances

2007-04-11 Thread Dimitris Rizopoulos
In fact this implies that random-effects might not be the way to go for your data. When you're using random-effects the marginal covariance matrix is of the form: V = Z D Z^t + Sigma, where Z is the design matrix for the random-effects, D their covariance matrix and Sigma is the covariance

Re: [R] Multivariate GARCH model in R

2007-04-11 Thread chiya sharma
Heard that I can't use multivariate GARCH model in R because R has only univariate GARCH models So, how can I run a multivariate GARCH model in R? There is a package in R for bivariate GARCH You can try search on MGARCHBEKK. Hope it helps, Gaurav [[alternative HTML version

[R] Programming Problem (for loop, random # control, 3 dimentional graph)

2007-04-11 Thread Mohammad Ehsanul Karim
Dear List, This is just a programming problem which i cannot seem to figure out. I am trying to get a set of power from a test (say, kolmogorov smirnov) out of a distribution (say, G-K distribution) as follows. I am trying to reduce to pain of writing the whole set of data points (p# below) using

Re: [R] Memory management

2007-04-11 Thread Prof Brian Ripley
Start with the 'R Internals' manual. R has 'call by value' semantics, but lazy copying (the idea is to make a copy only when an object is changed and there are still references to the original version, but that idea is partially implemented). 'which strategy is better at which situation' is

Re: [R] Memory management

2007-04-11 Thread Charilaos Skiadas
Before you go down that road, I would recommend first seeing if it is really a problem. Premature code optimization is in my opinion never a good idea. Also, reading the Details on ?attach you will find this: The database is not actually attached. Rather, a new environment is

Re: [R] Random Forest Imputations [Broadcast] [Broadcast]

2007-04-11 Thread Liaw, Andy
The package has a doc/ subdirectory (in the pre-compiled package, or inst/doc in the source package), which contains yaImputePaper.pdf. Page 9 of that document may be of some help to you. This is the first time I've seen this package, so can't help you much there. It looks like the package

[R] Kmeans cluster analysis

2007-04-11 Thread Monica Pisica
Hi Nataniel, As far as i know there is a package called clustTool which has a very nice interface with the capability to do different cluster analyses. It also prodused a plot of each cluster and the mean for each cluster of each variable - and i guess this is what you are after! But depending

[R] graphs superimposed on pictures?

2007-04-11 Thread Robert Biddle
Hi: I am doing some work that involves plotting points of interest superimposed on photographs and maps. I can produce the plots fine in R, but so far I have had to do the superimposition externally, which makes it tedious to do exploratory work. I have looked to see if there is some capability

[R] MAFA and DFA for short time series

2007-04-11 Thread Chabot Denis
Hi, recently I heard a talk where MAFA (min/max autocorrelation factor analysis) and DFA (dynamic factor analyses) were used to analyse short time series typical of fisheries data. I searched on internet to learn more about these techniques and (hopefully) to learn how to perform them in

Re: [R] Computing fundamental harmonics from a periodogram

2007-04-11 Thread hadley wickham
On 4/11/07, Uwe Ligges [EMAIL PROTECTED] wrote: hadley wickham wrote: Dear all, I'm trying to finding the fundamental harmonics (ie. peaks in a periodogram) from a time series (extracted from an mp3). For example, if I look at spectrum(fdeaths, spans = c(3,3)) A heuristical

Re: [R] Programming Problem (for loop, random # control, 3 dimentional graph)

2007-04-11 Thread Kuhn, Max
Wow. Some suggestions are below. 1. You should vectorize the for loop in the function g. I would do something like this: z1 - rnorm(20 * 100, 0, 1) Q1=0+1*z1*(1+.83*(1-exp(-G*z1))/(1+exp(-G*z1)))*(1+z1**2)**K Q1Matrix - matrix(Q1, ncol = 100) ksStat - function(x) ks.test(x, pnorm, mean =

[R] Why warnings using lmer-model with family=binomial

2007-04-11 Thread Åsa Granberg
Hi all! My question is why, and what I can do about that I sometimes, but not always, get warning-messages like nlminb returned message singular convergence (7) in: LMEopt(x = mer, value = cv) or IRLS iterations for PQL did not converge when trying to fit a model looking like this:

Re: [R] Kmeans cluster analysis

2007-04-11 Thread paulandpen
Hi Nataniel, A quick and easy way to do this is to run a decision tree analysis (like chaid) or discriminant analysis using the cluster variables you used as the predictors and the cluster membership as the outcome variable. The decision tree will highlight in a monotonic sense which

Re: [R] Programming Problem (for loop, random # control, 3 dimentional graph)

2007-04-11 Thread Kuhn, Max
Sorry, there was a typo in this line: function(x) ks.test(x, pnorm, mean = mean(x), sd = sd(x))$p.value Max -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Kuhn, Max Sent: Wednesday, April 11, 2007 9:55 AM To: [EMAIL PROTECTED];

Re: [R] sas.get problem

2007-04-11 Thread John Kane
--- [EMAIL PROTECTED] wrote: John, According to the sas.get documentation, the format library option must specify the directory containing the file formats.sct (I don't think it can be a sas dataset). Ah yes. I think I had it in one attempt and then decided I had misread the

Re: [R] Reasons to Use R

2007-04-11 Thread Douglas Bates
On 4/10/07, Wensui Liu [EMAIL PROTECTED] wrote: Greg, As far as I understand, SAS is more efficient handling large data probably than S+/R. Do you have any idea why? SAS originated at a time when large data sets were stored on magnetic tape and the only reasonable way to process them was

Re: [R] crosstabs in R?

2007-04-11 Thread Sarah Goslee
I don't know what the S-Plus crosstabs() does, but there's a simple crosstab() function in the ecodist package. Sarah On 4/10/07, A.R. Criswell [EMAIL PROTECTED] wrote: Hello All, I looked through the archives and I've seen this question asked before, but with no answer: Is there an R

Re: [R] graphs superimposed on pictures?

2007-04-11 Thread Gabor Grothendieck
Check out these posts: http://finzi.psych.upenn.edu/R/Rhelp02a/archive/94605.html http://finzi.psych.upenn.edu/R/Rhelp02a/archive/66704.html On 4/11/07, Robert Biddle [EMAIL PROTECTED] wrote: Hi: I am doing some work that involves plotting points of interest superimposed on photographs

Re: [R] MAFA and DFA for short time series

2007-04-11 Thread Stéphane Dray
Not sure about mafa and brodgar but ade4 has multispati. This performs a multivariate analysis which maximizes the product of autocorrelation by variance. I have a manuscript submitted to JVS (in revision) where the method is presented that I could send; Chabot Denis wrote: Hi, recently I

Re: [R] negative variances

2007-04-11 Thread Douglas Bates
On 4/11/07, Tu Yu-Kang [EMAIL PROTECTED] wrote: Dear R experts, I had a question which may not be directly relevant to R but I will be grateful if you can give me some advices. I ran a two-level multilevel model for data with repeated measurements over time, i.e. level-1 the repeated

[R] Reshape and column order

2007-04-11 Thread Ian Kennedy
If I try to reshape a data frame into a long format with more than one column it seems to mix up the column headings (or the columns, depending on how you look at it). For example: d - data.frame(row=1:2,b.1=c(1b1,2b1),b.2=c(1b2,2b2),a.1=c(1a1,2a1),a.2=c(1a2,2a2)) d row b.1 b.2 a.1 a.2 1 1

Re: [R] sas.get problem

2007-04-11 Thread John Kane
--- Tim Churches [EMAIL PROTECTED] wrote: John Kane wrote: I have 3 SAS files all in the directory F:/sas, two data files and a format file : form.ea1.sas7bdat form.ea2.sas7bdat sas.fmts.sas7bdat F is a USB. I am using R2.4.1 under Windows XP SAS files were created with

Re: [R] Computing fundamental harmonics from a periodogram

2007-04-11 Thread Uwe Ligges
hadley wickham wrote: On 4/11/07, Uwe Ligges [EMAIL PROTECTED] wrote: hadley wickham wrote: Dear all, I'm trying to finding the fundamental harmonics (ie. peaks in a periodogram) from a time series (extracted from an mp3). For example, if I look at spectrum(fdeaths, spans =

Re: [R] graphs superimposed on pictures?

2007-04-11 Thread Earl F. Glynn
This example isn't quite what you're asking for, but perhaps it's a starting point: http://addictedtor.free.fr/misc/gabor/ Earl F. Glynn Stowers Institute Robert Biddle [EMAIL PROTECTED] wrote in message news:[EMAIL PROTECTED] I am doing some work that involves plotting points of interest

Re: [R] Reasons to Use R

2007-04-11 Thread Mike Prager
Certainly true. In particular, SAS was designed from to store data items on disk, and to read into core memory the minimum needed for a particular calculation. The kind of data SAS handles is (for the most part) limited to rectangular arrays, similar to R data frames. In many procedures they can

Re: [R] graphs superimposed on pictures?

2007-04-11 Thread David Forrest
On Wed, 11 Apr 2007, Robert Biddle wrote: Hi: I am doing some work that involves plotting points of interest superimposed on photographs and maps. I can produce the plots fine in R, but so far I have had to do the superimposition externally, which makes it tedious to do exploratory work.

[R] persp Error. x,y,z with the same length

2007-04-11 Thread Felix Wave
Hello, I have a problem in drawing a 3d graphic. I simplified the program to show you the problem. I have a text file with hundreds of entries in 3 columns. I scaned these columns with matrix(scan). Then, I transformed some datas. Now I have the datas in h,x,y. I created a new matrix to use

[R] Reasons to Use R

2007-04-11 Thread Alan Zaslavsky
Right: SAS objects (at least in the base and statistics components of the system -- there are dozens of add-ons for particular markets) are simple databases. the predominant model for data manipulation and statistical calculation is a row by row operation that creates modified rows and/or

Re: [R] graphs superimposed on pictures?

2007-04-11 Thread Robert Biddle
Thanks; I don't have the actual book, but had looked at the online examples and not seen anything. Now I go back and see he does have one example in chapter 3 (using a photo of the moon) and I just hadn't noticed it. I think I got frustrated because I was typically searching for image, and R

Re: [R] graphs superimposed on pictures?

2007-04-11 Thread Christos Hatzis
Take a look at Figure 3.26 from Paul Murrell's book and associated code posted here: http://www.stat.auckland.ac.nz/~paul/RGraphics/chapter3.html -Christos -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Robert Biddle Sent: Wednesday, April 11,

Re: [R] Reasons to Use R

2007-04-11 Thread Bi-Info (http://members.home.nl/bi-info)
I certainly have that idea too. SPSS functions in a way the same, although it specialises in PC applications. Memory addition to a PC is not a very expensive thing these days. On my first AT some extra memory cost 300 dollars or more. These days you get extra memory with a package of

Re: [R] Using Sampling Weights in R

2007-04-11 Thread Thomas Lumley
On Tue, 10 Apr 2007, Thomas W. Volscho wrote: Dear List, I have a dataset that provides sampling weights (National Survey of Family Growth 2002). I want to produce a cross-tabulation and use the provided sampling weights to obtain representative population estimates. (I believe they are

Re: [R] graphs superimposed on pictures?

2007-04-11 Thread Greg Snow
Here are some examples that may be what you are looking for: library(rimage) x - read.jpeg(system.file(data, cat.jpg, package=rimage)) plot(x) par(new=TRUE) plot( runif(100), rnorm(100), col='green' ) x2 - x/4 + 3/4 plot(x2) locator() # click on corners to get extent of picture

Re: [R] graphs superimposed on pictures?

2007-04-11 Thread David McArthur
Robert Biddle robert_biddle at carleton.ca writes: Hi: I am doing some work that involves plotting points of interest superimposed on photographs and maps. I can produce the plots fine in R, but so far I have had to do the superimposition externally, which makes it tedious to do

Re: [R] persp Error. x,y,z with the same length

2007-04-11 Thread Duncan Murdoch
Felix Wave wrote: Hello, I have a problem in drawing a 3d graphic. I simplified the program to show you the problem. I have a text file with hundreds of entries in 3 columns. I scaned these columns with matrix(scan). Then, I transformed some datas. Now I have the datas in h,x,y. I

Re: [R] Reasons to Use R

2007-04-11 Thread Rajarshi Guha
On Wed, 2007-04-11 at 11:06 -0400, Alan Zaslavsky wrote: I have thought for a long time that a facility for efficient rowwise calculations might be a valuable enhancement to S/R. The storage of the object would be handled by a database and there would have to be an efficient interface for

Re: [R] Reasons to Use R

2007-04-11 Thread Duncan Temple Lang
Rajarshi Guha wrote: On Wed, 2007-04-11 at 11:06 -0400, Alan Zaslavsky wrote: I have thought for a long time that a facility for efficient rowwise calculations might be a valuable enhancement to S/R. The storage of the object would be handled by a database and there would have to be an

Re: [R] persp Error. x,y,z with the same length

2007-04-11 Thread Greg Snow
x and y should be vectors with the coordinates, z should be a matrix with the heights at the combination of those coordinates. If x and y are properly gridded in your dataset, then try something like: h - Measure3[,4] x - sort(unique(Measure3[,3])) y - sort(unique(Measure3[,2])) z -

[R] creating a path diagram in sem

2007-04-11 Thread John Smith
Hello, I finally run my measurement model in sem - successfully. Now, I am trying to print out the path diagram that is based on the results - but for some reason it's not working. Below is my script - but the problem is probably in my very last line: # ANALYSIS OF ANXIETY, DEPRESSION, AND

Re: [R] Reasons to Use R

2007-04-11 Thread Greg Snow
-Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Alan Zaslavsky Sent: Wednesday, April 11, 2007 9:07 AM To: R-help@stat.math.ethz.ch Subject: [R] Reasons to Use R [snip] I have thought for a long time that a facility for efficient rowwise

[R] Significance p.values for partial.cor

2007-04-11 Thread e . rapsomaniki
Dear R users Is there an easy way to compute the significance values for the correlation matrix obtained from partial cor(in Rcmdr) using method=kendall? Any help or suggestions would be much appreciated! Many thanks Eleni Rapsomaniki __

Re: [R] Reasons to Use R

2007-04-11 Thread Marc Schwartz
On Wed, 2007-04-11 at 17:56 +0200, Bi-Info (http://members.home.nl/bi-info) wrote: I certainly have that idea too. SPSS functions in a way the same, although it specialises in PC applications. Memory addition to a PC is not a very expensive thing these days. On my first AT some extra memory

Re: [R] sas.get problem

2007-04-11 Thread Cody_Hamilton
John, I believe the format file must have a .sc or .sas7bcat file extension (not .sas7bdat which is a sas dataset extension). I think that's why you're getting the error F:/sas/formats.sc? or formats.sas7bcat not found. Also, is the carriage return after 'C:/Program in the code below caused

Re: [R] Reasons to Use R

2007-04-11 Thread Alan Zaslavsky
thanks, I will take a look. __ R-help@stat.math.ethz.ch 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

[R] Error with corCompSymm and lme fit for repeated measures

2007-04-11 Thread Francois Teste
Dear R Friends, I need help with an error associated with corCompSymm in an lme fit. I am using a mixed effects model to analyze a split-plot with repeated measures and would like to fit with the compound symmetry correlation structure. This problem doesn't occur when using corAR1 or any of the

[R] Questions about statistics: ARIMA modeling seasonal MA term close to -1

2007-04-11 Thread sj
Hello, I am trying to fit a seasonal ARIMA model. I am using the ARIMA package in R. when I look at the model coefficients I find that that the first order seasonal moving average term is very close to -1(this is after first order seasonal differencing), I have tried to look at the literature on

Re: [R] creating a path diagram in sem

2007-04-11 Thread John Fox
Dear Dmitri, path.diagram() takes a fitted sem object as its first argument, not a model-specification object: path.diagram(sem.anxiety1,minrank='a1,a2,a3,d1,d2,d3,f1,f2,f3', maxrank='ANXIETY,DEPRESS,FEAR') digraph sem.anxiety1 { rankdir=LR; size=8,8; node [fontname=Helvetica fontsize=14

Re: [R] Reasons to Use R [Broadcast]

2007-04-11 Thread Liaw, Andy
From: Douglas Bates On 4/10/07, Wensui Liu [EMAIL PROTECTED] wrote: Greg, As far as I understand, SAS is more efficient handling large data probably than S+/R. Do you have any idea why? SAS originated at a time when large data sets were stored on magnetic tape and the only

Re: [R] sas.get problem

2007-04-11 Thread John Kane
--- [EMAIL PROTECTED] wrote: John, I believe the format file must have a .sc or .sas7bcat file extension (not .sas7bdat which is a sas dataset extension). I think that's why you're getting the error F:/sas/formats.sc? or formats.sas7bcat not found. I looked at that about 4 times

Re: [R] Reasons to Use R

2007-04-11 Thread Marc Schwartz
On Wed, 2007-04-11 at 11:26 -0500, Marc Schwartz wrote: On Wed, 2007-04-11 at 17:56 +0200, Bi-Info (http://members.home.nl/bi-info) wrote: I certainly have that idea too. SPSS functions in a way the same, although it specialises in PC applications. Memory addition to a PC is not a very

Re: [R] graphs superimposed on pictures?

2007-04-11 Thread Michael Lawrence
I think the previous suggestions have solved your problem, but just in case you need even more flexibility, you could try using the cairo graphics device (cairoDevice) in conjunction with the GdkPixbuf bindings provided by RGtk2. This allows fairly advanced image manipulation operations (scaling,

Re: [R] sas.get problem

2007-04-11 Thread Cody_Hamilton
John, Do you not have sas on your machine? (That would definitely keep sas.get from executing.) Regards, -Cody John Kane [EMAIL

[R] Help with pexp( )

2007-04-11 Thread Jeann S
Dear all, Sorry for bringing up an old issue: pexp(50, 0.5) [1] 1 In some cases, pexp() gives CDF=1. I read some discussion in 2002 saying it has been patched. However it's not working in R2.4.1Patched. Could anyone help me out? Thanks a lot, Jeann

Re: [R] sas.get problem

2007-04-11 Thread John Kane
--- [EMAIL PROTECTED] wrote: John, I believe the format file must have a .sc or .sas7bcat file extension (not .sas7bdat which is a sas dataset extension). I think that's why you're getting the error F:/sas/formats.sc? or formats.sas7bcat not found. Current state of play: Files : SAS

Re: [R] sas.get problem

2007-04-11 Thread John Kane
--- [EMAIL PROTECTED] wrote: John, Do you not have sas on your machine? (That would definitely keep sas.get from executing.) Regards, -Cody No I don't have SAS on my machine. Nobody in my group uses it and we seldom use SPSS. I think one person has STATA. Most of the work is more

[R] sem is not taking the model

2007-04-11 Thread John Smith
A strange problem with sem: I downloaded the sem library and then, I specified my simple measurement model (below). I highlighted it and ran it. It ran, but it did NOT tell me 22 lines read. And nothing works after that - it looks like it runs, but it does not produce anything... Did I make a

[R] User defined grid on map

2007-04-11 Thread Adriana Bejarano
This may be a trivial question for some of you: Is there a way to add defined grid sizes on specific areas of a map?. I want to add 10kmx10km grids to all California coastal counties: california,alameda, california,contra costa, california,del norte, california,humboldt, california,los

Re: [R] Help with pexp( )

2007-04-11 Thread Douglas Bates
On 4/11/07, Jeann S [EMAIL PROTECTED] wrote: Dear all, Sorry for bringing up an old issue: pexp(50, 0.5) [1] 1 In some cases, pexp() gives CDF=1. I read some discussion in 2002 saying it has been patched. However it's not working in R2.4.1Patched. Could anyone help me out? And

Re: [R] Positioning in xyplot

2007-04-11 Thread Sundar Dorai-Raj
Michael Kubovy said the following on 4/10/2007 5:54 PM: On Apr 10, 2007, at 7:07 PM, Deepayan Sarkar wrote: On 4/10/07, Michael Kubovy [EMAIL PROTECTED] wrote: Dear R-helpers, I have an xyplot of the following kind: xYplot(y ~ x | p, groups = factor(gg, levels = c('1', '5', '2', '4',

Re: [R] Help with pexp( )

2007-04-11 Thread Prof Brian Ripley
On Wed, 11 Apr 2007, Douglas Bates wrote: On 4/11/07, Jeann S [EMAIL PROTECTED] wrote: Dear all, Sorry for bringing up an old issue: pexp(50, 0.5) [1] 1 In some cases, pexp() gives CDF=1. I read some discussion in 2002 saying it has been patched. However it's not working in

Re: [R] Help with pexp( )

2007-04-11 Thread Jeann S
Hi Douglas, Thank you very much. Why I need the exact value is that I need apply qnorm() after. If it pexp() gives 1, then qnorm produces Inf. Your suggestion is very helpful. In my case, I'll jsut simply use -qnorm(pexp(..., lower.tail=FALSE)). All the best, Jeann From: Douglas Bates

Re: [R] Help with pexp( )

2007-04-11 Thread Ted Harding
On 11-Apr-07 15:59:09, Jeann S wrote: Dear all, Sorry for bringing up an old issue: pexp(50, 0.5) [1] 1 In some cases, pexp() gives CDF=1. I read some discussion in 2002 saying it has been patched. However it's not working in R2.4.1Patched. Could anyone help me out? Thanks

Re: [R] Why warnings using lmer-model with family=binomial

2007-04-11 Thread Douglas Bates
On 4/11/07, Åsa Granberg [EMAIL PROTECTED] wrote: Hi all! My question is why, and what I can do about that I sometimes, but not always, get warning-messages like nlminb returned message singular convergence (7) in: LMEopt(x = mer, value = cv) or IRLS iterations for PQL did not

Re: [R] sas.get problem

2007-04-11 Thread Marc Schwartz
On Wed, 2007-04-11 at 14:38 -0400, John Kane wrote: --- [EMAIL PROTECTED] wrote: John, Do you not have sas on your machine? (That would definitely keep sas.get from executing.) Regards, -Cody No I don't have SAS on my machine. Nobody in my group uses it and we seldom

Re: [R] sas.get problem

2007-04-11 Thread John Kane
--- Marc Schwartz [EMAIL PROTECTED] wrote: On Wed, 2007-04-11 at 14:38 -0400, John Kane wrote: --- [EMAIL PROTECTED] wrote: John, Do you not have sas on your machine? (That would definitely keep sas.get from executing.) Regards, -Cody No I don't have

Re: [R] Positioning in xyplot

2007-04-11 Thread Deepayan Sarkar
On 4/11/07, Sundar Dorai-Raj [EMAIL PROTECTED] wrote: Seems like you may get a workaround (albeit kludgey) by using ?print.trellis. Here's another example: library(lattice) z - expand.grid(x = 1:10, p = 1:5, r = 1:10) z$y - rnorm(nrow(z)) z$p - factor(z$p, levels = c(1, 5, 2, 4, 3)) bot -

Re: [R] Reasons to Use R

2007-04-11 Thread Robert Duval
So I guess my question is... Is there any hope of R being modified on its core in order to handle more graciously large datasets? (You've mentioned SAS and SPSS, I'd add Stata to the list). Or should we (the users of large datasets) expect to keep on working with the present tools for the time

Re: [R] sas.get problem

2007-04-11 Thread Tim Churches
John Kane wrote: How do I make this change? I naively have tried by a) list sas.get and copy to editor b) reload R without loading Hmisc c) made recommended changes to sas.get d) stuck a sas.get - in front of the function and ran it. Here is what I do, until Frank fixes the problem in the

Re: [R] Reasons to Use R

2007-04-11 Thread Wensui Liu
I think the reason that stata is fast is because it only keeps 1 work table in ram. if you just keep 1 data frame in R, it will run fast too. But ... On 4/11/07, Robert Duval [EMAIL PROTECTED] wrote: So I guess my question is... Is there any hope of R being modified on its core in order to

  1   2   >