Re: [R] Help with loess "Standard Error of the Residuals"

On May 1, 2012, at 07:27 , Bert Gunter wrote: > John: > > 1. This is not really an R question. You should post on statistics > help list , e.f. stats.stackexchange.com > If it is about the concrete implementation, then it is actually an R question. It probably won't harm to take a look at the

Re: [R] Help with loess "Standard Error of the Residuals"

John: 1. This is not really an R question. You should post on statistics help list , e.f. stats.stackexchange.com Also, inline below. -- Bert On Mon, Apr 30, 2012 at 8:44 AM, John Haslett wrote: > Dear All > > I'm having trouble working out what exactly loess means by its "Standard > Error of

Re: [R] Questions

This is not the RStudio support forum. I use it, I like it, but you need to learn the difference between R and RStudio. In this case, the fact that you are using RStudio is probably irrelevant. What is more relevant is that you are not on windows, though I can't be sure what OS you actually are

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

Break the line :) On Mon, Apr 30, 2012 at 12:34 PM, Zeeberg, Barry (NIH/NCI) [E] 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? > > ___

Re: [R] read multiple dataset at one time

Without more information on the format of the files, or any other thing, it is hard to give a full answer. But shortly: yes, it is possible. With a bit more details: You can combine: readLines (or read.table, or what ever will work for your file), with list.files And run through them with a for lo

Re: [R] IV estimation

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 wrote: > Hello, > I have a set of 100 variables with 1560 observations. I did an O.L.S > regression of three of these variable

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

On Mon, Apr 30, 2012 at 12:34 PM, Zeeberg, Barry (NIH/NCI) [E] 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 our WGCNA package w

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

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] Scatter plot / LOESS, or LOWESS for more than one parameter

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. http://sites.google.com/site/davidsstatis

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

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? > > Th

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

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 > bootstra

[R] IV estimation

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 endogeneous

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

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

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

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. T

[R] read multiple dataset at one time

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 ht

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

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 f

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

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

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:

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

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 boots

[R] Extracting coefficients values with bootstrap

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 volc

Re: [R] Matrix transposition

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-help@r-p

Re: [R] compiling an R script

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: http://r.789695.n4.nabble.com/c

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

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 <- as

Re: [R] Matrix transposition

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) > >

Re: [R] Gls function in rms package

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 (

[R] Matrix transposition

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 correspo

[R] Regression

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

Re: [R] Different varable lengths

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 https

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

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 dime

[R] Clustering analysis with ordination plots

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

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

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 . > > Her

[R] Extracting week number starting from a specific date

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 fo

[R] Robust regression with clusters

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 <- ro

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

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] Calling Rscript from Matlab on OS X

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". Th

Re: [R] Generate Dendrogram

farscape2012 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 example > shown in

[R] Help with loess "Standard Error of the Residuals"

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 wei

Re: [R] Min , Max

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,

Re: [R] for loop problem

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 :) Tha

Re: [R] for loop problem

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

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)

[R] for loop problem

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 3

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

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"),

Re: [R] Different varable lengths

Because of the missing comma. On Mon, Apr 30, 2012 at 4:56 PM, John Kane 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 Sweden.GDP.gapp.ad

Re: [R] Using GEE with sample weights

On Mon, Apr 30, 2012 at 10:05 PM, Stefán Hrafn Jónsson 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 attrition. > > Aft

Re: [R] Different varable lengths

> 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] Exporting plots generated by a 'for' loop

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] RODBC Error Code 202 on Mac OS X 10.6

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" wrote:

Re: [R] Generate Dendrogram

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.

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

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) > t

Re: [R] Generate Dendrogram

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 i

Re: [R] R2 in multilevel modelling

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

[R] [R-pkgs] RExcelInstaller update

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 list

Re: [R] Area between 2 curves

A quick Google suggests that sintegral is in the Bolstad2 package: http://cran.r-project.org/web/packages/Bolstad2/ "R. Michael Weylandt" wrote in message news:CAAmySGOCtNxWNHJDka=68criphxnppwap6fm7ufkhtx7xsz...@mail.gmail.com... > You can't just decide a command exists and try to use it by you

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

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 <- cbind(x,cols

Re: [R] Area between 2 curves

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 i

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

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

[R] log= argument to pairs()?

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 > s

Re: [R] R2 in multilevel modelling

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 at

Re: [R] R2 in multilevel modelling

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. Standar

Re: [R] Min , Max

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 w

Re: [R] Different varable lengths

Missing comma between -1.180620213 and -0.525964648 "Saint" 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 objects > (Sweden.GD

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

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] Error in if (nuhat < 2) stop("The degrees of freedom must be greater than or equal to 2") : missing value where TRUE/FALSE needed

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) bilingualdataforconferences2012$TypeBiling<-facto

Re: [R] CForest Error Logical Subscript Too Long

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

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

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) ano

Re: [R] Min , Max

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 tw

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

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

Re: [R] Different varable lengths

> 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 th

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

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: lig...@statistik.tu-dortmund.

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

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:

Re: [R] Gls function in rms package

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. Fr

Re: [R] Using GEE with sample weights

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

[R] R2 in multilevel modelling

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 multilevel

[R] Different varable lengths

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 acc

[R] Area between 2 curves

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 q<-abs(

[R] HOw compare 2 models in logistic regression

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] Generate Dendrogram

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 http://stat.ethz.ch/R-manual/R

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

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, col.names=

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

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 W

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

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") bilingualismdata<-read.sps

[R] "parallel" package

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

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 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 them

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

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 s

[R] Using GEE with sample weights

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 sample

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

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