Re: [R] ANOVA 1 too few degrees of freedom

On May 5, 2011, at 23:30 , Rovinpiper wrote: > Thanks slre, > > I seem to be making some progress now. > > Using a colon instead of an asterisk in the code really changes things. I > had been getting residual SS and MS of zero. Which is ridiculous. Now I get > much more plausible values. > >

[R] Bug in lattice that shipped with R 2.13.0

Hi all, I had meant to make this announcement earlier but had forgotten (a recent bug report reminded me). The version of lattice that ships with R 2.13.0 has a fairly serious bug in panel.abline (which would neglect to draw many negative slope lines). If you use lattice with R 2.13.0, you should

Re: [R] memory and bootstrapping

(Regarding bootstrapping logistic regression.) If the number of rows with Y=1 is small, it doesn't matter that n is huge. If both number of successes and failures is huge, then ad Ripley notes you can use asymptotic CIs. The mean difference in predicted probabilities is a nonlinear function of t

Re: [R] Vermunt's LEM in R

I guess LEM is a software for latent class analysis. If so, you may want to have a look at poLCA package. Regards Ronggui On 5 May 2011 23:34, David Joubert wrote: > > Hello- > > Does anyone know of packages that could emulate what J. Vermunt's LEM does ? > What is the closest relative in R ? >

Re: [R] MacKinnon critical value

On May 5, 2011, at 9:10 PM, Lee Schulz wrote: Hello, I am doing an Engle Granger test on the residuals of two I(1) processes. I would like to get the MacKinnon (1996) critical value, say at 10%. I have 273 observations with 5 integrated explanatory variables , so that k=4. Could som

Re: [R] Installing rgdal in R: correct -configure flags for GDAL install on Linux Redhat

Hi Katrina, the error message below is actually pretty explicit. Have you installed the PROJ.4 library? If not, then you need to install it. When I had to do this I think that I used macports: sudo port install proj Then you need to tell configure where to find it, using the protocol suggested

Re: [R] help with a survplot

Hi Marco, You're welcome. The number at risk at given time points is a fairly standard thing to add to survival plots. I don't think many people will confuse it with numbers of events. Personally I find shaded confidence bands a bit more helpful but both are useful. Frank Marco Barbàra-2 wro

[R] MacKinnon critical value

Hello, I am doing an Engle Granger test on the residuals of two I(1) processes. I would like to get the MacKinnon (1996) critical value, say at 10%. I have 273 observations with 5 integrated explanatory variables , so that k=4. Could someone help me with the procedure in R? Thank you

[R] Installing rgdal in R: correct -configure flags for GDAL install on Linux Redhat

Hi, I'm installing rgdal but I keep having failures because I have not been able to find a good source of information for the correct configuration settings when installing GDAL. My error from the R install.packages("rgdal") is below. Can someone point me to a good source to tell me how to set af

Re: [R] two-way group mean prediction in survreg with three factors

Even then, I think that there's a problem. If C is in the model, then the response varies by C. The simplest way is to pick a value for C, and then evaluate the group mean estimates of A and B (and C). Something in my brain keeps asking whether another way to marginalize C for the purposes of pr

Re: [R] Vermunt's LEM in R

Hi David, you might have more luck with your request if you tell us what Vermunt's LEM *does*, and provided some links to introductory reading material ... Cheers Andrew On Thu, May 05, 2011 at 03:34:02PM +, David Joubert wrote: > > Hello- > > Does anyone know of packages that could emula

Re: [R] functions pandit and treebase in the package apTreeshape

Hi Arnau, please send the output of sessionInfo() and the exact commands and response that you used to install and load apTreeshape. Cheers Andrew On Thu, May 05, 2011 at 04:42:58PM +0200, Arnau Mir wrote: > Hello. > > I'm trying to use the functions pandit and treebase. They are in the packag

Re: [R] Using GSUB to obtain a printing "\"

gsub does this because the string "dsff\nfsd" does not contain a backslash - the 5th character is a newline. The deparsed representation (used for printing the string) of a newline is "\n" but the string itself has not backslash. You can feed the output of deparse into message (or cat) so they sh

Re: [R] Draw a nomogram after glm

Thanks Frank I will try rms package and give after the result. Komine -- View this message in context: http://r.789695.n4.nabble.com/Draw-a-nomogram-after-glm-tp3498144p3499771.html Sent from the R help mailing list archive at Nabble.com. __ R-he

[R] Using GSUB to obtain a printing "\"

Hello. I'd like to be able to print variable strings which contain "\" as-is, without interpreting (for example) "abcde\nuvxyz" as having an embedded newline (or whatever other escaped instruction). To do this, I've tried gsub, and here's some of my output (I've tried all kinds of variations t

Re: [R] factors

as.character() will enable me to do this in a loop (so thanks!). It looks like you guys are suggesting I accomplish this in some kind of vectorized fashion (David: yes, the number of rows are not equal - and I need to compare ALL pairs, not just adjacent pairs, so to vectorize, I suppose I'd have

Re: [R] factors

On Thu, May 5, 2011 at 6:28 PM, J wrote: > Hi, I'm requesting you don't berate me for asking this question: > > I clearly don't have the gist of factors. > > I have two dataframes, A and B. > > Each of them has a column containing strings (they're labels). > > I want to, one-by-one in a loop, comp

[R] [caret package] [trainControl] supplying predefined partitions to train with cross validation

Hi all, I run R 2.11.1 under ubuntu 10.10 and caret version 2.88. I use the caret package to compare different models on a dataset. In order to compare their different performances I would like to use the same data partitions for every models. I understand that using a LGOCV or a boot type re-sam

Re: [R] hurdle, simulated power

Dave: We are planning an intervention study for adolescent alcohol use, and I am planning to use simulations based on a hurdle model (using the hurdle() function in package pscl) for sample size estimation. The simulation code and power code are below -- note that at the moment the "power" c

Re: [R] factors

Hi Jonathan, On Thu, May 5, 2011 at 6:28 PM, J wrote: > Hi, I'm requesting you don't berate me for asking this question: > > I clearly don't have the gist of factors. > > I have two dataframes, A and B. > > Each of them has a column containing strings (they're labels). > > I want to, one-by-one i

Re: [R] factors

On May 5, 2011, at 6:28 PM, J wrote: Hi, I'm requesting you don't berate me for asking this question: I clearly don't have the gist of factors. I have two dataframes, A and B. Each of them has a column containing strings (they're labels). I want to, one-by-one in a loop, compare the particu

[R] factors

Hi, I'm requesting you don't berate me for asking this question: I clearly don't have the gist of factors. I have two dataframes, A and B. Each of them has a column containing strings (they're labels). I want to, one-by-one in a loop, compare the particular string in an entry from dataframe A

Re: [R] Using $ accessor in GAM formula

On 06/05/11 07:05, P Ehlers wrote: Gene, David has given you the preferred code. I just want to point out that the $-accessor is often not the best thing to use. Both dat[["y"]] and dat[, "y"] will work just fine. Admittedly one should use the preferred code, i.e. gam(y ~ s(x),data=dat) and av

Re: [R] Looking for equivalent for "strstr"

On Thu, May 5, 2011 at 4:55 PM, Worik R wrote: > Friends > > This is an elementary question.  Is there is a built in R function for > finding a sub-string in another string?  Like strstr in C. > > I can easily roll my own, but if there is a built in that is one less thing > I can do wrong! > Try

Re: [R] quick question : interpolating file name in pipe command

On May 5, 2011, at 5:15 PM, Abhishek Pratap wrote: Hi Guys I am trying to read a bunch of files in the loop but pipe function which I use to cut few columns is somehow unable to interpolate the file variable. eg: file="check.txt" data <- read.table(pipe("cut -f 2,3 file"), sep="\t", col

Re: [R] quick question : interpolating file name in pipe command

You can ignore my question I was able to figure out the way. I guess when I touch R after couple of weeks I am rusty. -Abhi On Thu, May 5, 2011 at 2:15 PM, Abhishek Pratap wrote: > Hi Guys > > I am trying to read a bunch of files in the loop but pipe function > which I use to cut few columns i

[R] reading a column as a character vector

Hi! I have 2 columns (even though the data looks like there's more columns than just two) of data in the following format: 0,"58905313R0EOL" "229742002R0EOL" "58905312R0EOL" 1,"58905317R0DBL" "58905303R0DBL" "58905313R0IL" "58905313R0VH" "58905313R0EOL" "223354003R0IL" "223354003R0VH" "58905308R0

Re: [R] ANOVA 1 too few degrees of freedom

Thanks slre, I seem to be making some progress now. Using a colon instead of an asterisk in the code really changes things. I had been getting residual SS and MS of zero. Which is ridiculous. Now I get much more plausible values. Also, When I used an asterisk instead of a colon It wouldn't give

Re: [R] Insert values to histogram

Histograms plot data in bins - you don't get the exact value, because each bin contains a range of values. Do you want to plot the range of values the bin contains? Also, check ?hist to see how to set the values of the breaks between the bins. __ R-help@

[R] quick question : interpolating file name in pipe command

Hi Guys I am trying to read a bunch of files in the loop but pipe function which I use to cut few columns is somehow unable to interpolate the file variable. eg: > file="check.txt" > data <- read.table(pipe("cut -f 2,3 file"), sep="\t", > col.names=c('pos','cov') ) cut: file: No such file or

[R] Looping over graphs in igraph

Hi, I'm trying to do some basic social network analysis with igraph in R, but I'm new to R and haven't been able to find documentation on a couple basic things: I want to run igraph's community detection algorithms on a couple thousand small graphs but don't know how to automate igraph looking a

Re: [R] Looking for equivalent for "strstr"

On May 5, 2011, at 4:55 PM, Worik R wrote: Friends This is an elementary question. Is there is a built in R function for finding a sub-string in another string? Like strstr in C. I can easily roll my own, but if there is a built in that is one less thing I can do wrong! I have no acqua

Re: [R] Averaging uneven measurements by time with uneven numbers of measurements

On 2011-05-05 14:20, Schatzi wrote: I do not want smoothing as the data should have jumps (it is weight left in feeding bunker). I was thinking of maybe using a histogram-like function and then averaging that. Not sure if this is possible. (It would be useful to include your original request -

Re: [R] Averaging uneven measurements by time with uneven numbers of measurements

In your first request for help you said, "How do I average the measurements by every 30 minutes?" With 25000 readings over three days, it looks as if you are getting readings just about every second. Okay, why don't you use the first reading as your initial weight, w0. Then subtract each su

[R] Looking for equivalent for "strstr"

Friends This is an elementary question. Is there is a built in R function for finding a sub-string in another string? Like strstr in C. I can easily roll my own, but if there is a built in that is one less thing I can do wrong! cheers Worik [[alternative HTML version deleted]] __

[R] problem with cor() using bigmemory

Hello, I have a rather large set of data I need to analyze, currently I need to work with a 20 by 20 matrix, I'm using the package bigmemory but so far I can only allocate a 66000 by 66000 matrix, when I increase those values I get the following error:   > AdjMat <- big.matrix(nrow=68

Re: [R] Averaging uneven measurements by time with uneven numbers of measurements

I do not want smoothing as the data should have jumps (it is weight left in feeding bunker). I was thinking of maybe using a histogram-like function and then averaging that. Not sure if this is possible. - In theory, practice and theory are the same. In practice, they are not - Albert Einstei

Re: [R] Averaging uneven measurements by time with uneven numbers of measurements

I'd be tempted to do a robust fit (loess?) to the data with a relatively small span (I'm assuming that there are errors in the measurements and some degree of smoothing is acceptable) then predict the fit at a regular interval (e.g., every 30 minutes). -- Clint BowmanINTERN

[R] Averaging uneven measurements by time with uneven numbers of measurements

I have a new device that takes measurements anywhere from every second, to every 15 minutes (depending on changes). The matrix has a date, time and Y column (Y is the measurement). For three days it is 25,000 rows. How do I average the measurements by every 30 minutes so my matrix is 48 rows per da

Re: [R] R CMD check warning

On 05/05/2011 11:06 AM, swaraj basu wrote: Dear All, I am trying to build a package for a set of functions. I am able to build the package and its working fine. When I check it with R CMD check I get a following warning : no visible global function def

Re: [R] R CMD check warning

On 6/05/2011 6:06 a.m., swaraj basu wrote: Dear All, I am trying to build a package for a set of functions. I am able to build the package and its working fine. When I check it with R CMD check I get a following warning : no visible global function def

Re: [R] Instrumental variable quantile estimation of spatial autoregressive models

This is a very detailed question, and possibly ought to have been posted on R-sig-geo. Please re-post there, and see if anyone steps forward. As far as I am aware, no implementation exists. I suggest you also refer to Kostov (2009) in Spatial Economic Analysis volume 4. Roger Marie-Line Glaesen

Re: [R] help with a survplot

Il giorno Mon, 2 May 2011 07:30:10 -0700 (PDT) Frank Harrell ha scritto: > Please elaborate. It is simply that, generally speaking, i don't like adding numbers to a plot. I eventually realized that riskset cardinality may be a useful indication. However, do not discriminate between events and ce

Re: [R] Using $ accessor in GAM formula

Gene, David has given you the preferred code. I just want to point out that the $-accessor is often not the best thing to use. Both dat[["y"]] and dat[, "y"] will work just fine. Peter Ehlers On 2011-05-05 12:06, David Winsemius wrote: On May 5, 2011, at 1:08 PM, Gene Leynes wrote: This is

Re: [R] ANOVA

On May 5, 2011, at 2:31 PM, Asan Ramzan wrote: Hello R-Help How can i exctact and store the "within group mean squared difference" from an anova summary table into a varible. In the absence of an example and code it's speculation, but something along the lines of: anova(fit)$`Mean Sq`

Re: [R] ANOVA 1 too few degrees of freedom

>>> Rovinpiper 04/05/2011 22:43 >>> >So this seems to indicate that I have what I want. I have two >respiration data points at each plot on each day. Yes; if you had only Plot+Day you'd have a completely balanced full factorial ... for Plot and Day. But I think I now see an answer to your puzz

[R] ANOVA

Hello R-Help How can i exctact and store the "within group mean squared difference" from an anova summary table into a varible. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-he

[R] R CMD check warning

Dear All, I am trying to build a package for a set of functions. I am able to build the package and its working fine. When I check it with R CMD check I get a following warning : no visible global function definition for ‘biocLite’ I have use

Re: [R] Using $ accessor in GAM formula

On May 5, 2011, at 1:08 PM, Gene Leynes wrote: This is not mission critical, but it's bothering me. I'm getting inconsistent results when I use the $ accessor in the gam formula *In window #1:* library(mgcv) dat=data.frame(x=1:100,y=sin(1:100/50)+rnorm(100,0,.05)) str(dat) gam(dat$y~s(dat$x)

Re: [R] Using functions/loops for repetitive commands

On May 5, 2011, at 1:45 PM, dereksloan wrote: Thanks a lot, I understand what you say but I'm having problems - maybe with the syntax or the specific command. You are right - I have a dataframe to store the data and want to automate the analysis. i.e. I want do a chisq.test with to kno

Re: [R] cross-correlation table with subscript or superscript to indicate significant differences

Here is an example for my earlier question. Say you have a 3x3 correlation matrix:corrs <- matrix(c(0.25,0.32,0.66,0.14,0.24,0.34,0.44,0.34,0.11), nrow=3, ncol=3, dimnames = list(c('varA','varB', 'varC'), c('varA','varB', 'varC')))And another matrix for the sample size of each correlation:s

[R] Using $ accessor in GAM formula

This is not mission critical, but it's bothering me. I'm getting inconsistent results when I use the $ accessor in the gam formula *In window #1:* > library(mgcv) > dat=data.frame(x=1:100,y=sin(1:100/50)+rnorm(100,0,.05)) > str(dat) > gam(dat$y~s(dat$x)) Error in eval(expr, envir, enclos) : objec

Re: [R] cross-correlation table with subscript or superscript to indicate significant differences

On May 5, 2011, at 12:40 PM, yoav baranan wrote: Here is an example for my earlier question. Say you have a 3x3 correlation matrix: corrs <- matrix(c(0.25,0.32,0.66,0.14,0.24,0.34,0.44,0.34,0.11), nrow=3, ncol=3, dimnames = list(c('varA','varB', 'varC'), c('varA','varB', 'varC'))) And ano

Re: [R] distance matrix

I don't know whether I understood your question, but 1.1, 1.2 , 1.3 all are subsample of 1 , so, rather than comparing 1000 subsample, comparison of 20 pop level makes more sense in my case. thanks for query - Ananta Acharya Graduate Student -- View this message in context: http://r.78969

Re: [R] Using functions/loops for repetitive commands

Thanks a lot, I understand what you say but I'm having problems - maybe with the syntax or the specific command. You are right - I have a dataframe to store the data and want to automate the analysis. i.e. I want do a chisq.test with to know if alcohol intake (Y/N) differs between sexes, then if

Re: [R] Using functions/loops for repetitive commands

On May 5, 2011, at 1:08 PM, dereksloan wrote: Thanks David, I did notice that and I got his code to work using wilcox.test for the continuous variables. The problem is that when I tried to alter the code to do chisq.test on my categorical variables there is something wrong with the syntax a

Re: [R] quantmod's addTA plotting functions

Sorry, I forgot to give you a quick example of 'wrapping in plot' plot(addTA(...)) The objects returned by most of the charting functions in quantmod results from the desire for the functions to be syntactically identical whether called from inside of chartSeries (e.g. TA="addMACD()" ) as they ar

Re: [R] quantmod's addTA plotting functions

Russ, All you have to do is replace addTA(GSPC.EMA.3, on = 1, col = "#ff") with plot(addTA(GSPC.EMA.3, on = 1, col = "#ff")) etc. I can sympathize with the documentation frustration, but I think that much of the documentation in R and in many R packages is actually very good.

Re: [R] perspective plot

The persp function expects z to be a matrix, so you could reshape your data so that z is a matrix (the reshape function or package may help). Or the wireframe function in the lattice package expects data more like what you show, that may be the easiest solution. -- Gregory (Greg) L. Snow Ph.D

Re: [R] Using functions/loops for repetitive commands

Thanks David, I did notice that and I got his code to work using wilcox.test for the continuous variables. The problem is that when I tried to alter the code to do chisq.test on my categorical variables there is something wrong with the syntax and I don't know what. Derek -- View this message i

Re: [R] Boxplot in order

Hi: Try this: karla = data.frame( Groups = factor(rep(c('CPre','SPre','C7','S7','C14','S14','C21','S21'), 11), levels = c('CPre','SPre','C7','S7','C14','S14','C21','S21')), Time = rep(c(0,7,14,21), 11), Resp = rnorm(88) ) boxplot(Resp~Groups, data = karla) Since you didn't h

Re: [R] quantmod's addTA plotting functions

There is a struggle in documentation that revolves around being too brief to be useful and too verbose which then is often ignored. In general, R "proper" is far less verbose than quantmod docs - so if you have trouble with quantmod... from ?addTA Value: ‘addTA’ will invisibly return an S4

Re: [R] Insert values to histogram

Are you really sure that you want to do that? Read the discussion starting with this post: http://tolstoy.newcastle.edu.au/R/e2/help/07/08/22858.html for reasons why you probably don't (yes, the question is about bar plots not histograms, but much of it will still apply). Near the end of the d

Re: [R] Null

with(fcv,shapiro.test(case)) -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FAX:(360) 407-7534 Olymp

Re: [R] Compiling a FORTRAN program under Windows 7

Mikael Anderson wrote: > > Hi, > > I am trying to compile a FORTRAN program to call from R under Windows 7 > but > I am having problem in the compiling step. To demonstrate this is the > program testit.f: > > -- > subroutine TESTIT(x,n,m) > di

Re: [R] quantmod's addTA plotting functions

Thanks. You're right. I didn't see that. I read the ?addTA help page, which (annoyingly) didn't mention that feature, but I didn't read the ?TA page. (That page was mentioned as a see also, but not as a must see.) I don't know what it means to wrap these calls in a plot call. I tried to put the a

Re: [R] Using error in histograms

The logspline package does density estimation in a different way than KDE, but it does allow for interval censored data (I know this value is between a and b, but not where in that range) using the oldlogspline function. This may be what you are looking for. -- Gregory (Greg) L. Snow Ph.D. St

Re: [R] Null

Maybe you were doing something like fcv <- read.csv('fcv.csv') instead of fcv <- read.csv('fcv.csv')[1] (I haven't tested this.) Tom On Thu, May 5, 2011 at 8:48 AM, pcc wrote: > > This is probably a very simple question but I am completely stumped!I am > trying to do shapiro.wilk(x) test on

Re: [R] simulate AR(1) process

?arima.sim - Phil Spector Statistical Computing Facility Department of Statistics UC Berkeley

Re: [R] Compiling a FORTRAN program under Windows 7

You are compiling a subroutine not a program and you compile line should read: gfortran testit.f -c testit.o You then reference that object code testit.o in your final loading stage after compiling other routiens and the main program. -- Clint BowmanINTERNET: cl...@

[R] simulate AR(1) process

Dear R users, May any of you tell me how to simulate data on: y_t = a+b*y_{t-1} + u_t where u_t~N(0,sigma^2), b<1, and for some constant a. Many thanks Tsegaye [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://

Re: [R] cross-correlation table with subscript or superscript to indicate significant differences

On May 5, 2011, at 10:48 AM, yoav baranan wrote: Hi, I wonder whether the following is possible with R, and whether anyone has done that and can share his/her code with me. I have a correlation matrix, and I want to create a correlation table that I can copy to Microsoft Word with a supe

[R] perspective plot

Hello, I`m tryint to plot some variable over x and y coordinate, but can`t figure out hot to do it ( I could do a simple scatterplot3d, but there I can`t change the viewing angle). My dataset looks something like this (dataframe): X Y Q95 21 2628711 1104437 0.7723994 22 2628

Re: [R] Confidence interval for difference in Harrell's c statistics (or equivalently Somers' D statistics)

Dear David, Thank you for your reply. I have come across rcorrp.cens before. However, I'm not sure it does quite what I want it to. It seems to compare whether one predictor is more concordant than another within the same survival function. I want to see whether one predictor is more concordan

[R] Compiling a FORTRAN program under Windows 7

Hi, I am trying to compile a FORTRAN program to call from R under Windows 7 but I am having problem in the compiling step. To demonstrate this is the program testit.f: -- subroutine TESTIT(x,n,m) dimension x(n) do 10 i=1,n 10 x(i)=x(i

[R] Vermunt's LEM in R

Hello- Does anyone know of packages that could emulate what J. Vermunt's LEM does ? What is the closest relative in R ? I use both R and LEM but have trouble transforming my multiway tables in R into a .dat file compatible with LEM. Thanks, David Joubert

Re: [R] Null

There were no problems when I repeated the test with data provided. I simply created a vector <-c(1.600972896,1.534026106,1.633468456,1.69019608,1.686636269,1.713490543,1.460897843,1.604226053,1.547774705,1.575187845,1.50242712,1.489958479,1.555094449,1.56937391,1.46686762,1.583198774,1.5943925

[R] cross-correlation table with subscript or superscript to indicate significant differences

Hi, I wonder whether the following is possible with R, and whether anyone has done that and can share his/her code with me. I have a correlation matrix, and I want to create a correlation table that I can copy to Microsoft Word with a superscript above each correlation, indicating significant d

[R] Boxplot in order

Hi, I need construct box plot graph, but I want keep Groups order karla = data.frame( Groups = factor(rep(c('CPre','SPre','C7','S7','C14','S14','C21','S21'), 11)), Time = rep(c(0,7,14,21), 11), Resp = valor ) boxplot(Resp~Groups, order=T) doesn't work. How do this? -

[R] Tone in mailing lists (was " issue with "strange" characters (readHTMLTable)")

I did read How To Ask Questions The Smart Way and I don't have a problem with calling me stupid and/or a winer for posting this... What's the purpose of an R Mailing list (especially r-help), after all? IMHO it should be about R users being

Re: [R] distance matrix

Em 5/5/2011 07:34, antu escreveu: Hello all, I am wondering if there is anyway to create distance matrix for replicated data for example, I have a data like sample pop id var1var2var3var4 1.1 1 a 1 1 0 1 1.2 1 a 0 0

Re: [R] Confidence interval for difference in Harrell's c statistics (or equivalently Somers' D statistics)

On May 5, 2011, at 8:20 AM, Laura Bonnett wrote: Dear All, I am trying to calculate a 95% confidence interval for the difference in two c statistics (or equivalently D statistics). In Stata I gather that this can be done using the lincom command. Is there anything similar in R? Have y

Re: [R] Using functions/loops for repetitive commands

On May 5, 2011, at 10:01 AM, dereksloan wrote: Your code may be untested but it works - also helping me slowly to start understanding how to write functions. Thank you. However I still have difficulty. I also have some categorical variables to analyse by age & hiv status - i.e. my dataset

[R] Conditional distribution plot using Model-based Recursive Partitioning

Hello, I am using the party module to estimate the relationship between the probability of being a student and number of siblings (alive). However, I need to include a number of relevant covariates. My code is below: fm3 <- mob(Student ~ age + alive + sex2 + cwa + cha + cym | Religion+Servant + L

Re: [R] How to fit a random data into Beta distribution?

Beta is not as general as you think. Its support is limited to [0,1], but you are trying to fit data that lies outside of its support. Please read about the beta distribution from a basic stats/prob book. Ravi. From: r-help-boun...@r-project.org [r-he

Re: [R] two-way group mean prediction in survreg with three factors

Oops, I hope not too. Don't know why I had the brackets around B+C. My model is actually A*B+C. And I'm not sure how to obtain the two-way prediction of AB with C marginalized. Thanks. Pang -Original Message- From: Andrew Robinson [mailto:a.robin...@ms.unimelb.edu.au] Sent: Wednesday

[R] RV: R question

which is the maximum large of digits that R has?, because SQL work with 50 digits I think. and I need a software that work with a lot of digits. Thanks. Pamela Santelices Elgueta Estadístico Instituto Nacional de Estadísticas Fono: (56-2) 7962491 Paseo Bulnes 209 oficina 82 Santiago [

[R] Confidence interval for difference in Harrell's c statistics (or equivalently Somers' D statistics)

Dear All, I am trying to calculate a 95% confidence interval for the difference in two c statistics (or equivalently D statistics). In Stata I gather that this can be done using the lincom command. Is there anything similar in R? As you can see below I have two datasets (that are actually two i

Re: [R] Using functions/loops for repetitive commands

Hi Derek, You can accomplish your loop jobs by following means: (a) use for loop (b) use while loop (c) use lapply, tapply, or sapply. (i feel "lapply is the elegant way ) ---For Loop- "for" loops are pretty simple to use and is almost similar to any other

[R] functions pandit and treebase in the package apTreeshape

Hello. I'm trying to use the functions pandit and treebase. They are in the package apTreeshape. Once I've loaded the package, R responses: - no function pandit/treebase. Somebody knows why or what is the reason? Thanks, Arnau. Arn

Re: [R] Remove all whitespaces

A more elegant way would be: myString<-"1 2 3 4 5" myString<-paste(unlist(strsplit(myString," ")),collapse="") The output will be "12345" Regards, Som Shekhar __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] How to fit a random data into Beta distribution?

Hi, @Steven: Since Beta distribution is a generic distribution by which i mean that by varying the parameter of alpha and beta we can fit any distribution. So to check this i generated a random data from Normal distribution like x.norm<-rnorm(n=100,mean=10,sd=10); Now i want to estimate the par

[R] Null

This is probably a very simple question but I am completely stumped!I am trying to do shapiro.wilk(x) test on a relatively small dataset(75) and each time my variable and keeps coming out as 'NULL', and > shapiro.test(fcv) Error in complete.cases(x) : no input has determined the number of cases

[R] distance matrix

Hello all, I am wondering if there is anyway to create distance matrix for replicated data for example, I have a data like sample pop id var1var2var3var4 1.1 1 a 1 1 0 1 1.2 1 a 0 0 1 0 1.3 1 a

Re: [R] heatmap.3

I am intesrested by your heatmap function (allowing matrix in ColSideColors option). Can you give your complete code of your function? thanks -- View this message in context: http://r.789695.n4.nabble.com/heatmap-3-tp1566584p3498156.html Sent from the R help mailing list archive at Nabble.com. _

[R] lm weight adj r-sq

Dear all, First of all apologies, I'm pretty newbie, and secondly I know this is not a question for the forum but I'd be very grateful if you'd help me. I had problems with normality assumptions, and I tried with weights in the linear regression this way: modelw <-lm(log(v1)~v2+log(v2)+log(v3

[R] Insert values to histogram

I'm trying to add the exact value on top of each column of an Histogram, i have been trying with the text function but it doesn't work. The problem is that the program it self decides the exact value to give to each column, and ther is not like in a bar-plot that I know exactly which values are bee

[R] Using error in histograms

Hello! I am trying to produce a histogram of measurement data (orientation of archaeological structures) that are a subject to measurement error. The normal histogram just computes frequencies, but does not take into account that a particular value is spread over a range of values (in my case the

Re: [R] Convert a presence/ absence matrix to a list of presence only

Thank you, it perfectly worked, except I had to modify some codes on the rep("", maxLen - length(sp)) because it gave me some error, thanks -- View this message in context: http://r.789695.n4.nabble.com/Convert-a-presence-absence-matrix-to-a-list-of-presence-only-tp3468479p3498116.html Sent fro

[R] Draw a nomogram after glm

Hi all R users I did a logistic regression with my binary variable Y (0/1) and 2 explanatory variables. Now I try to draw my nomogram with predictive value. I visited the help of R but I have problem to understand well the example. When I use glm fonction, I have a problem, thus I use lrm. My code

  1   2   >