[R] Why does my R compiler repeat my program whenever it is compiled?

2012-04-30 Thread jpm miao
Hi, I am using RStudio as my R editor. After someday I accidentally hit something, the whole program is repeated in the Console whenever I compile it. How can I fix it so that the whole program won't be repeated in the future? Thanks, miao [[alternative HTML version deleted]]

[R] Using GEE with sample weights

2012-04-30 Thread Stefán Hrafn Jónsson
Dear R community I am using the gee package to run logistics regression on paired cases from a panel sample. We are getting request from a reviewer to use sample weights to account for different level of attrition. After searching the documentation I am unable to find a way to incorporate

Re: [R] Why does my R compiler repeat my program whenever it is compiled?

2012-04-30 Thread Steve Lianoglou
Hi, This looks like an rstudio question that's not specific to R itself There are support forums at rstudio.org which you should ask your question on, ani think you'll get better help there On Monday, April 30, 2012, jpm miao wrote: Hi, I am using RStudio as my R editor. After

Re: [R] parallel package

2012-04-30 Thread Uwe Ligges
On 29.04.2012 12:01, ya wrote: Hi Uwe and Indrajit, Thank you very much for the response. Since we are talking about this package, can I ask more about how to use it to deal with statistical issues please? I have this data with categorical missing values in it. Now I am trying to impute

[R] parallel package

2012-04-30 Thread ya
Hi Uwe, Thank you very much for the suggestions. ya ya From: Uwe Ligges Date: 2012-04-30 15:18 To: ya CC: Indrajit Sengupta; r-help@r-project.org Subject: Re: [R] parallel package On 29.04.2012 12:01, ya wrote: Hi Uwe and Indrajit, Thank you very much for the response. Since we are

Re: [R] Error in if (nuhat 2) stop(The degrees of freedom must be greater than or equal to 2) : missing value where TRUE/FALSE needed

2012-04-30 Thread Uwe Ligges
On 29.04.2012 16:28, Kyriakos Antoniou wrote: Hi, i am trying to run an ANCOVA and a bootstrapped ANCOVA analysis on a specific data set. I am using the ancova and ancboot functions as in the following code: setwd(C:/Users/User/Desktop/Rdatabilingualstudy2012)

Re: [R] Problem extracting enough coefs from gam (mgcv package)

2012-04-30 Thread Martijn Wieling
Dear useRs, Just to let you know, the results using mgcv 1.7-13 are OK. @Simon: could you perhaps still let me know if the procedure I apply w.r.t. assigning the labels to random slopes and intercepts is correct? With kind regards, Martijn -- *** Martijn

Re: [R] need help with avg.surv (Direct Adjusted Survival Curve), Message-ID:

2012-04-30 Thread Terry Therneau
Well, I would suggest using the code already in place in the survival package. Here is my code for your problem. I'm using a copy of the larynx data as found from the web resources for the Klein and Moeschberger book. larynx - read.table(larynx.dat, skip=12,

[R] Generate Dendrogram

2012-04-30 Thread farscape2012
Hi I have a distance matrix which is computed by user defined method. I would like to plot the dendrogram. I would like to use different color and want the leaves laying down bottom. The script like this. I am not familiar with R. I followed the example shown in

[R] HOw compare 2 models in logistic regression

2012-04-30 Thread santoshdvn
Hi, I have created 2 models in logistic regression and got the predictor values are significance on output. Here are 2 summaries of 2 model. HOw can we compare 2 models by what factor or coefficient and say which model is best Please help

[R] Area between 2 curves

2012-04-30 Thread sappy
Hello, i have a question calculating the shaped area between the two curves (see image). http://r.789695.n4.nabble.com/file/n4597813/test.png I try to use a Simpson-Integral but it doesn't work. R doesn't know the command! It is possible, that i need a library? Curve 1: y Curve 2: z

[R] Different varable lengths

2012-04-30 Thread Saint
Hi! I'm trying to do a lm() test on three objects. My problem is that R protests and says that the variable lengths differ for one of the objects (Sweden.GDP.gap). But I have double checked that the number of observations are the same. All three objects should contain 9 observations but R only

[R] R2 in multilevel modelling

2012-04-30 Thread klai...@libero.it
Goodmorning everybody, i'm an italian statistician and i'm using R for research. Could someone tell me some indices to see the goodness of fit in multilevel modelling? I'm using the lmer function, and I want to know if my model fit well my data. I actually want to justify the use of

Re: [R] Using GEE with sample weights

2012-04-30 Thread Haynes, Maurice (NIH/NICHD) [E]
See the function geeglm() in package:geepack Maurice O. Maurice Haynes, Ph.D. Statistician, Child and Family Research Section Laboratory for Comparative Ethology Eunice Kennedy Shriver National Institute of Child Health and Human Development NIH, DHHS 6705 Rockledge Drive, Suite 8030 Bethesda,

Re: [R] Gls function in rms package

2012-04-30 Thread Frank Harrell
Hi Mark, I'm sorry I forgot about this. I found that Gls was out of date with respect to gls. It's now fixed and you can source('http://biostat.mc.vanderbilt.edu/tmp/Gls.s') to override the old version with the fixed version until the next release of rms. You won't need the workaround below.

Re: [R] Error in if (nuhat 2) stop(The degrees of freedom must be greater than or equal to 2) : missing value where TRUE/FALSE needed

2012-04-30 Thread Kyriakos Antoniou
Sorry!! Data now attached. Note that this is a .dat file while in the R code i import data directly from an SPSS file (.sav). Best, Kyriakos Date: Mon, 30 Apr 2012 15:20:04 +0200 From: lig...@statistik.tu-dortmund.de To: antonio...@hotmail.com CC: r-help@r-project.org Subject: Re: [R]

Re: [R] Error in if (nuhat 2) stop(The degrees of freedom must be greater than or equal to 2) : missing value where TRUE/FALSE needed

2012-04-30 Thread Kyriakos Antoniou
Dear Uwe. Many thanks for your reply. I attach the data set that i use for the analyses. Hopefully this will help you identify the problem. I'm really new in using r so please bear with me. All the best, Kyriakos Date: Mon, 30 Apr 2012 15:20:04 +0200 From:

Re: [R] Different varable lengths

2012-04-30 Thread Petr PIKAL
Hi! I'm trying to do a lm() test on three objects. My problem is that R protests and says that the variable lengths differ for one of the objects (Sweden.GDP.gap). But I have double checked that the number of observations You shall check it again. BTW How did you checked it? are the

Re: [R] Error in if (nuhat 2) stop(The degrees of freedom must be greater than or equal to 2) : missing value where TRUE/FALSE needed

2012-04-30 Thread Uwe Ligges
On 30.04.2012 15:52, Kyriakos Antoniou wrote: Sorry!! Data now attached. Note that this is a .dat file while in the R code i import data directly from an SPSS file (.sav). Oh, come on, we do not have endless time to import your data, so inlcuding code please. Moreover, you have not

Re: [R] Min , Max

2012-04-30 Thread Jonsson
Hi again, Sorry for connecting you again(I am still getting this error:subscripts out of bounds). the problem is that I am trying to apply R script for one file to many files. It succeeded to do that for the first file but not for the others.I think the problem in my code is that ( d ) (in the

Re: [R] HOw compare 2 models in logistic regression

2012-04-30 Thread ONKELINX, Thierry
Dear Santos, Since your models are nested you can apply a likelihood ratio test. M0 - glm(formula = Spend_bucket ~ Freq + Address_is_res + last_update_days_ago, data = qdataset, family = binomial) M1 - glm(formula = Spend_bucket ~ Freq + Address_is_res, data = qdataset, family = binomial)

Re: [R] CForest Error Logical Subscript Too Long

2012-04-30 Thread Achim Zeileis
On Sun, 29 Apr 2012, Abigail Clifton wrote: Hi, This is my code (my data is attached): library(languageR) library(rms) library(party) OLDDATA - read.csv(/Users/Abigail/Documents/OldData250412.csv) OLDDATA$YD - factor(OLDDATA$YD, label=c(Yes, No))  OLDDATA$ND - factor(OLDDATA$ND, label=c(Yes,

Re: [R] Error in if (nuhat 2) stop(The degrees of freedom must be greater than or equal to 2) : missing value where TRUE/FALSE needed

2012-04-30 Thread Kyriakos Antoniou
Dear Uwe, really sorry about this. The code i use is the following (data attached again): setwd(C:/Users/User/Desktop/Rdatabilingualstudy2012) bilingualdataforconferences2012-read.delim(bilingualdataforconferences2012.dat, header=TRUE)

Re: [R] Error in if (nuhat 2) stop(The degrees of freedom must be greater than or equal to 2) : missing value where TRUE/FALSE needed

2012-04-30 Thread Uwe Ligges
As far as I understand you want: lmo - lm(Simon ~ TypeBiling + FAS, data = bilingualdataforconferences2012) anova(lmo) summary(lmo) for the ordinary ancova. Not sure if there is a proper implementation for that kind of bootstrapped robust ancova you are using. You may have to fix the code

Re: [R] Different varable lengths

2012-04-30 Thread Keith Jewell
Missing comma between -1.180620213 and -0.525964648 Saint bahman.za...@comhem.se wrote in message news:1335779012548-4597768.p...@n4.nabble.com... Hi! I'm trying to do a lm() test on three objects. My problem is that R protests and says that the variable lengths differ for one of the

Re: [R] Min , Max

2012-04-30 Thread Jonsson
Yes, But what is your opinion about d in sprintf(C:\\Users\\aalyaari\\Desktop\\New folder (6)\\Yar_%00d.bin, d), because will go from one to the number of rows(some files have 248 and others have 224) and my results will change accordingly. for example i will get from the first file 31 files

Re: [R] R2 in multilevel modelling

2012-04-30 Thread Bert Gunter
As you have not provided a clue of what your models are, one can only guess. But if you mean using lots of fixed effects vs a random effect, the answer is that there is no such animal. They are two different non-nested models, and should be chosen based on subject matter considerations.

Re: [R] R2 in multilevel modelling

2012-04-30 Thread Joshua Wiley
Hi Chiara, If you just want to compare model fit, you could use a LRT between models where you do and do not estimate the variance/covariance matrix of random effects. R^2 in mixed models do not have the same nice properties they do in fixed effects models. Cheers, Josh On Mon, Apr 30, 2012

[R] log= argument to pairs()?

2012-04-30 Thread Szumiloski, John
Dear useRs, I would like to be able to smoothly use the log base graphics argument as an argument to pairs(). Currently in my hands, pairs() accepts the argument, makes transformed plots accordingly, but then stumbles over some annotation commands. Example: begin session

[R] question on jitter in plot.Predict in rms

2012-04-30 Thread Mike Babyak
Dear colleagues, I have a question regarding controlling the jitter when plotting predictions in the rms package. Below I've simulated some data that reflect what I'm working with. The model predicts a continuous variable with an ordinal score, a two-level group, and a continuous covariate. Of

Re: [R] Area between 2 curves

2012-04-30 Thread R. Michael Weylandt
You can't just decide a command exists and try to use it by your arbitrarily chosen name did you look at the function integrate() which does exist? It uses a smarter algorithm than Simpson's rule. Alternatively, Simpson's rule is very easy to code with vectorization -- you should be able to

Re: [R] Subtract days to dates in POSIXct format

2012-04-30 Thread Jeff Newmiller
On Mon, 30 Apr 2012, Santiago Guallar wrote: Hello, I'm having problems working with date values in POSIXct format. Indeed you are. Here is what I got (eg.lig attached): x - read.table(eg.txt, sep = ',', col.names=c(ok,time,secs,lig)) # it gives time as factor z -

Re: [R] Area between 2 curves

2012-04-30 Thread Keith Jewell
A quick Google suggests that sintegral is in the Bolstad2 package: http://cran.r-project.org/web/packages/Bolstad2/ R. Michael Weylandt michael.weyla...@gmail.com wrote in message news:CAAmySGOCtNxWNHJDka=68criphxnppwap6fm7ufkhtx7xsz...@mail.gmail.com... You can't just decide a command exists

[R] [R-pkgs] RExcelInstaller update

2012-04-30 Thread Neuwirth Erich
New version RExcel 3.2.7-1 Download wrapper for new version of RExcel. ___ R-packages mailing list r-packa...@r-project.org https://stat.ethz.ch/mailman/listinfo/r-packages __ R-help@r-project.org mailing

Re: [R] R2 in multilevel modelling

2012-04-30 Thread Kevin Wright
If you want to know if your model fit will your data, then looking at residual-type plots is useful, as can be plotting the model-predictions and observed data together. You might also find interesting R2 statistics for mixed models by Matthew Kramer. Beware--as others have indicated, there is

Re: [R] Generate Dendrogram

2012-04-30 Thread Mercier Eloi
hc.obj is not a dendrogram. plot(dend1, xlab='',ylab='',sub='', nodePar=list(pch = c(1,NA), cex=0.8, lab.cex = 0.8)) See ?plot.hclust and ?plot.dendrogram. They are not using the same arguments. Regards, Eloi On 12-04-30 06:36 AM, farscape2012 wrote: Hi I have a distance matrix which

[R] Hmisc::pstamp, mfcol, and spacing

2012-04-30 Thread David Reiner
I am trying to make a nice 2x1 plot and add a timestamp with comment. The pstamp function from Hmisc works nicely when mfcol=c(1,1), but when mfcol=c(2,1), the stamp winds up in the wrong place: require('Hmisc') opar - par(mfcol=c(2,1)) plot(1:10) title(main=MAIN Title) plot(1:20)

Re: [R] Generate Dendrogram

2012-04-30 Thread farscape2...@live.com
Thank. I made a stupid mistake. Now it is working. On 4/30/2012 8:30 PM, Mercier Eloi wrote: hc.obj is not a dendrogram. plot(dend1, xlab='',ylab='',sub='', nodePar=list(pch = c(1,NA), cex=0.8, lab.cex = 0.8)) See ?plot.hclust and ?plot.dendrogram. They are not using the same arguments.

Re: [R] RODBC Error Code 202 on Mac OS X 10.6

2012-04-30 Thread MacQueen, Don
This question might better be asked on r-sig-db or r-sig-mac (I am not sure which one would better, but please not both). -Don -- Don MacQueen Lawrence Livermore National Laboratory 7000 East Ave., L-627 Livermore, CA 94550 925-423-1062 On 4/29/12 7:08 AM, julia.jacob...@arcor.de

Re: [R] Exporting plots generated by a 'for' loop

2012-04-30 Thread Greg Snow
A couple of things that I did not see mentioned by the others: Generally statistics plots work better in .png files than in .jpg files due to the type of compression each uses (detailed image plots and some surface plots may be the exception), though that may be out of your control (some journals

Re: [R] Different varable lengths

2012-04-30 Thread John Kane
length(Sweden.GDP.gap.adjust) [1] 8 Sweden.GDP.gap.adjust [1] 0.673792123 1.196706756 1.196131539 0.646944002 -0.312886525 [6] -1.706584861 -0.369401194 -0.003280389 R is reading Sweden.GDP.gapp.adjust as having 8 elements. Why I don't know. John Kane Kingston ON Canada -Original

Re: [R] Using GEE with sample weights

2012-04-30 Thread Thomas Lumley
On Mon, Apr 30, 2012 at 10:05 PM, Stefán Hrafn Jónsson shj...@gmail.com wrote: Dear R community I am using the gee package to run logistics regression on paired cases from a panel sample. We are getting request from a reviewer to use sample weights to account for different level of

Re: [R] Different varable lengths

2012-04-30 Thread Sarah Goslee
Because of the missing comma. On Mon, Apr 30, 2012 at 4:56 PM, John Kane jrkrid...@inbox.com wrote: length(Sweden.GDP.gap.adjust) [1] 8 Sweden.GDP.gap.adjust [1]  0.673792123  1.196706756  1.196131539  0.646944002 -0.312886525 [6] -1.706584861 -0.369401194 -0.003280389 R is reading

Re: [R] Subtract days to dates in POSIXct format

2012-04-30 Thread Jeff Newmiller
The sunrise and sunset are calculated for each time value on the input vector, but all times are treated in the TZ timezone. I can see that the Canary Islands are in the GMT timezone, so the example TZ is right. I can see that you entered the latitude and longitude backward (longitude is x),

[R] for loop problem

2012-04-30 Thread Schreiber, Stefan
Hi all, I was wondering if you can help me with the following situation: I have a data frame that includes weather station data for 30 years in the form: YEAR, MONTH, DAY, TEMP 1970, 01, 01, -15 ... 1999, 12, 31, -21 I would like to add another variable JULIAN that assigns the integers 1 to

Re: [R] for loop problem

2012-04-30 Thread R. Michael Weylandt
Perhaps this is easier: convert your data to Date objects and use the yday (year day) function from the lubridate package: i.e., something like this: jd - apply(DAT, 1, function(x) as.Date(paste(x, collapse = ), format = %Y %m %d) # Just a guess you might have to debug it a little yday(jd) It

Re: [R] for loop problem

2012-04-30 Thread Tyler Rinker
I don't really work with dates but thought I'd pass a solution on.  I think that there some great packages for handling dates though (lubridate) and you may want to convert your data to a true date instead of separate columns. # FUNCTION TO INDEX DATES date.int - function(month, year, day){

Re: [R] for loop problem

2012-04-30 Thread Schreiber, Stefan
Hey Michael and Tyler! Thanks for the input and the effort! The lubridate package is really awesome! What I did now is just this here: dat$DATE-as.Date(paste(dat$YEAR,dat$MONTH,dat$DAY,sep='-'),format='%Y-%m-%d') require(lubridate) dat$JULIAN-yday(dat$DATE) That's it. Nice and easy :)

Re: [R] Min , Max

2012-04-30 Thread Rui Barradas
Hello, again. Try to use the values in 'listfile' directly: # Wrong!!! for(h in length(listfile)) # Just last value ?!? {b=file.info(listfile[h])$size/67420/4 wind - file(listfile[h]), rb) # Parenthesis do not match. # Right for(h in listfile) {b=file.info(h)$size/67420/4 wind - file(h, rb)

[R] Help with loess Standard Error of the Residuals

2012-04-30 Thread John Haslett
Dear All I'm having trouble working out what exactly loess means by its Standard Error of the Residuals denoted s and in particular when the weights argument is invoked. For example, if the weights are weights are all =1, then s^2 is nearly sum sq res/ (n -1 - 'equiv num paras') If the

Re: [R] Generate Dendrogram

2012-04-30 Thread cberry
farscape2012 farscape2...@live.com writes: Hi I have a distance matrix which is computed by user defined method. I would like to plot the dendrogram. I would like to use different color and want the leaves laying down bottom. The script like this. I am not familiar with R. I followed the

[R] Calling Rscript from Matlab on OS X

2012-04-30 Thread Kent Conover
This posting is for the help-file record in that the problem has already been resolved via a misplaced discussion in the R-bug forum (Bug 14894). I place it here so make it easier for others to find. The problem involved an error in calling Rscript from Matlab on Mac OS X 10.6 Snow Leopard. The

[R] Optim (fct): Parameters=LowerBounds!?

2012-04-30 Thread barb
Hey, i am trying to do the MLE for Garch and have a problem with the optim function. Initally i tried optim with Method=BFGS. Reading trhough the forum i found out i would neet bounds. So i went on with Method=L-BFGS-B. But now my parameters equal the lower bounds. out - optim(par=initial,

[R] Robust regression with clusters

2012-04-30 Thread Aldo Córdova Palomera
Hi, I need to perform a robust regression with several dummy independent variables (i.e., binary), which are paired (clustered). The outcome (dependent variable) is continuous. I have tried using the robcov function with the next syntaxis: out1 - ols(y ~ x1 + x2 + x3 , x=T) robustclust -

[R] Extracting week number starting from a specific date

2012-04-30 Thread irene
Hello, I have a list of dates, such as dates- as.Date(c(1996-03-29,1996-05-30,1996-09-28,1996-05-09)) from which I would like to extract the week number for each date, with week n°1 being the week going from Dec 30th 1995 to Jan 6th 1996 (1995-30-12 to 1996-06-01). Any suggestion for a simple

Re: [R] kiteChart to show real values with scalebar

2012-04-30 Thread bodiless
On 04/28/2012 02:32 PM, Jim Lemon [via R] wrote: On 04/28/2012 04:12 AM, bodiless wrote: Dear R-users I hope someone could help me on this problem. I want to create a multiple kiteChart showing the real values with a scalebar on each indicating the scale . Here are some

[R] Clustering analysis with ordination plots

2012-04-30 Thread borinot
Hello to all, I'm new to R so I have a lot of problems with it, but I'll only ask the main one. I have clustered an environmental matrix with 2 different methods, and I'd like to plot them in a PCA and a db-RDA. I mean, I want see these clusters in the plots like points of differents colours,

[R] [R-pkgs] New package for morphology and smoothing in any number of dimensions: mmand

2012-04-30 Thread Jon Clayden
Dear all, I'm pleased to announce the release on CRAN of the mmand package (for Mathematical Morphology in Any Number of Dimensions). It provides functions for performing mathematical morphology (erode, dilate, etc.), smoothing, and other kernel-based operations on array-like objects of any

Re: [R] Different varable lengths

2012-04-30 Thread Saint
Thanks a lot! One gets blind sometimes. -- View this message in context: http://r.789695.n4.nabble.com/Different-varable-lengths-tp4597768p4598762.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

[R] Regression

2012-04-30 Thread Saint
Trying to do a regression for four variables. I get information only for two. The rest is marked NA. Why? And what does NA mean? Thanks! Sweden.infl.dev - c(0.2, 0.6, -1, -0.2, -1, -1.5, -0.7, 0.7, -0.1, 0.3) Sweden.GDP.gap - c(-0.024662769, -0.01519859, -0.027251109, -0.027129556,

[R] Matrix transposition

2012-04-30 Thread Filoche
Hi everyone. I want to transpose a data frame. Lets say the following DF: df = data.frame(matrix(ncol=4, nrow = 10)) df[,1] = c(rep(1,5),rep(2,5)) df[,2] = c(rep('a',4),rep('b',3),rep('c',3)) df[,3] = c(letters[5:14]) df[,4] = runif(10) I would like to form a data frame with each line

Re: [R] Gls function in rms package

2012-04-30 Thread alice.jones
Hello, I was wondering whether there had been any progress on this bug yet? I have checked for updates to the package and it doesn't seem that there have been any, but I might be mistaken! Similarly to Mark, I am getting the following error message when running a Gls with two (factor) covariates

Re: [R] Extracting week number starting from a specific date

2012-04-30 Thread Rui Barradas
Hello, Try # # Adapted from several posts to R-Help, # mostly by David Winsemius, also by Gabor Grothendiek. # nweek - function(x, format=%Y-%m-%d, origin){ if(missing(origin)){ as.integer(format(strptime(x, format=format), %W)) }else{ x -

Re: [R] Matrix transposition

2012-04-30 Thread Rui Barradas
Hello, Filoche wrote Hi everyone. I want to transpose a data frame. Lets say the following DF: df = data.frame(matrix(ncol=4, nrow = 10)) df[,1] = c(rep(1,5),rep(2,4), 3) df[,2] = c(rep('a',4),rep('b',3),rep('c',3)) df[,3] = c(letters[c(5:13,13)]) df[,4] = runif(10) I would

Re: [R] compiling an R script

2012-04-30 Thread strycker
I'm interested in this question, too, not so I can hide code, but so I can run R code faster. It's been my experience that compiled code always runs faster than interpreted code. Can you explain further how to compile a front-end? -- View this message in context:

Re: [R] Matrix transposition

2012-04-30 Thread Filoche
Thank you sire, this is exactly what I was looking for. Many thanks, Phil -- View this message in context: http://r.789695.n4.nabble.com/Matrix-transposition-tp4599219p4599310.html Sent from the R help mailing list archive at Nabble.com. __

[R] 95% confidence interval of the coefficients from a bootstrap analysis

2012-04-30 Thread baconbeach
Hello, I am doing a simple linear regression analysis that includes few variables. I am using a bootstrap analysis to obtain the variation of my variables to replacement. I am trying to obtain the coefficients 95% confidence interval from the bootstrap procedure. Here is my script for the

[R] Extracting coefficients values with bootstrap

2012-04-30 Thread baconbeach
Hello fellow R users, I am trying to extract the coefficient values during a bootstrap operation. Here is the list of my variables that I would like to extract the coefficient values from: (Intercept) LogRds_25k GeoRockbimodal

Re: [R] 95% confidence interval of the coefficients from a bootstrap analysis

2012-04-30 Thread baconbeach
Hi Rui, Thanks for your help!! It works perfectly, but when I call stor.confint, I obtain the list of confidence interval (10,000 times). Is there an easy way to summarize the results and getting only one 2.5% and 97.5% values for each variable? Thanks again for your help Steeve -- View this

[R] Questions

2012-04-30 Thread haugb002
Hello, I am having difficulty putting a frequency table on my clipboard in R studio. I do not know what the error message is. I was thinking it might have to do with changing a setting on my clipboard but I do not know how to do this. I will copy and paste the exact command and error message:

Re: [R] Subtract days to dates in POSIXct format

2012-04-30 Thread Santiago Guallar
Thanks a lot, Jeff Yes, I applied colsplit in package reshape. Sorry, I forgot to state it in the piece of code I attached. Actually, these are archival data from a geolocator attached to a seabird. The origin lies on Gran Canaria and data were gathered all along its migratory journey for one

[R] read multiple dataset at one time

2012-04-30 Thread ya
Hi, Is it possible to read, say 100 dat files in one time? I want to combine these 100 data set into one single dat file. Any suggestions? Thank you very much. ya [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] Hmisc::pstamp, mfcol, and spacing

2012-04-30 Thread Taimur Sajid
You could try legend(bottomright, legend = Sys.time()) to get the timestamp in the right place (or alternatively, using text()). The documentation for pstamp mentions that it has a problem under Windows when mfrow is used. The mar argument allows you to adjust the plot margins. See ?par.

[R] Usage line in .pdf version of manual is truncated

2012-04-30 Thread Zeeberg, Barry (NIH/NCI) [E]
I am creating a package, and a rather long usage line for one function becomes truncated in the generated .pdf manual. Is there a way that I can prevent that truncation? __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

[R] IV estimation

2012-04-30 Thread kebrab67
Hello, I have a set of 100 variables with 1560 observations. I did an O.L.S regression of three of these variables on a fourth. But there are problems of endogeneity... So I look in my dataset for instruments to do an IV. I can't find a good instrument because their correlation with my

Re: [R] 95% confidence interval of the coefficients from a bootstrap analysis

2012-04-30 Thread Rui Barradas
Hello, baconbeach wrote Hello, I am doing a simple linear regression analysis that includes few variables. I am using a bootstrap analysis to obtain the variation of my variables to replacement. I am trying to obtain the coefficients 95% confidence interval from the bootstrap

Re: [R] 95% confidence interval of the coefficients from a bootstrap analysis

2012-04-30 Thread Rui Barradas
Hello, baconbeach wrote Hi Rui, Thanks for your help!! It works perfectly, but when I call stor.confint, I obtain the list of confidence interval (10,000 times). Is there an easy way to summarize the results and getting only one 2.5% and 97.5% values for each variable? Thanks

Re: [R] Scatter plot / LOESS, or LOWESS for more than one parameter

2012-04-30 Thread David Doyle
First let me say thank you for all the help everyone gave me. Below is the code I came up with to look at Rainfall pH reading for Kentucky and Tennessee. Also there is a more *generic* ver of the code on my web page if anybody ever needs it in the future.

Re: [R] 95% confidence interval of the coefficients from a bootstrap analysis

2012-04-30 Thread baconbeach
Thanks again for your swift response!! With your last line, I get rowMeans(sapply(stor.confint, colMeans)) 2.5 %97.5 % 0.3256882 0.4604677 I need the values (2.5% and 97.5%) for each variable of my model. I don't think this what I am getting. This is what my script looks like now,

Re: [R] Usage line in .pdf version of manual is truncated

2012-04-30 Thread Peter Langfelder
On Mon, Apr 30, 2012 at 12:34 PM, Zeeberg, Barry (NIH/NCI) [E] zeebe...@mail.nih.gov wrote: I am creating a package, and a rather long usage line for one function becomes truncated in the generated .pdf manual. Is there a way that I can prevent that truncation? How long is rather long? In

Re: [R] IV estimation

2012-04-30 Thread Bert Gunter
Whatever your post means, it's about statistics, not R. Post to a statistics list like stats.stackexchange.com instead. -- Bert On Mon, Apr 30, 2012 at 3:24 PM, kebrab67 selamgetac...@gmail.com wrote: Hello, I have a set of 100 variables with 1560 observations. I did an O.L.S regression of