[R] Bootstrap

Dear Any, Can someone please inform me, if they have a code to estimate the varaince using bootstrap resampling method under a two stage cluster design. Thanks for all your help and time. Murthy.M.N., PhD, Student, University of Canterbury, New Zealand.

Re: [R] png problem

Hello, Thanks for the fast reply. Paul Murrell wrote: Hi Clément Calenge wrote: Dear R-users, I have a small problem with the function png(), when used with the argument colortype=pseudo.cube. png(toto.png, colortype=pseudo.cube) image(matrix(rnorm(1), 100, 100)) dev.off() R is blocked

[R] t test problem?

Hello, I got two sets of data x=(124738, 128233, 85901, 33806, ...) y=(25292, 21877, 45498, 63973, ) When I did a t test, I got two tail p-value = 0.117, which is not significantly different. If I changed x, y to log scale, and re-do the t test, I got two tail p-value = 0.042, which is

Re: [R] t test problem?

Hi Liu, before applying a t-test (or any test) you should first check if the assumptions of the test are supported by your data, i.e., in a t-test x and y must be normally distributed. I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public

[R] t test problem?

Hi, maybe your data are distributed according a log-normal distribution, so logs are normally distributed. But remerber the significancy of t test can applied only on log transformated data and not on original data. See basic hypothesis for t testing, in alternative use non-parametric methods to

Re: [R] t test problem?

Hi Dimitris, you are describing a more stringent requirement than the t-test actually requires. It's the sampling distribution of the mean that should be normal, and this condition is addressed by the Central Limit Theorem. Whether or not the CLT can be invoked depends on numerous factors,

[R]

Dear Any, Is there a fonction in R to change a string to uppercase ? Thanks for all your help __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Re: [R] ordered probit and cauchit

On Wednesday, Sep 22, 2004, at 03:42 Europe/London, roger koenker wrote: What is the current state of the R-art for ordered probit models, and more esoterically is there any available R strategy for ordered cauchit models, i.e. ordered multinomial alternatives with a cauchy link function. MCMC

Re: [R] Ever see a stata import problem like this?

I've had a similar problem once. What may have caused the problem then was a variate for which value lables had been defined for the highest and lowest values. What complicates things is that the file had been originally converted from SPSS to Stata. A workaround was to set convert.factor=FALSE

RE: [R]

?toupper p.s: By default, generally everything on this list *is* regarding R; hence it would be nice to see a more imaginative subject line. -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of [EMAIL PROTECTED] Sent: Wednesday, September 22, 2004 5:43 PM

RE: [R]

The following function will do it. But be warned it will only work if all the input strings are lower case. Make.To.Upper.Case-function(my.string){ paste(LETTERS[match(strsplit(my.string,)[[1]],letters)],collapse=) } Make.To.Upper.Case(jhjhaskjdakdsj) [1] JHJHASKJDAKDSJ -Original

Re: [R] convert string to uppercase; was: nothing

[EMAIL PROTECTED] wrote: Dear Any, Is there a fonction in R to change a string to uppercase ? Thanks for all your help __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

Re: [R]

help.search(upper) chartr(base)Character Translation and Casefolding ?chartr refers to the function toupper() Best wishes, Arne On Wednesday 22 September 2004 10:43, [EMAIL PROTECTED] wrote: Dear Any, Is there a fonction in R to change a string to uppercase ? Thanks for all

Re: [R] t test problem?

Hi, Many thanks for your helpful comments and suggestions. The attached are the data in both log10 scale and original scale. It would be very grateful if you could suggest which version of test should be used. By the way, how to check whether the variation is additive (natural scale) or

[R] Create bitmap graphic in C

Hello, How to create a R graphic bitmap with C or C++ using R api ? Thanks, samples, ... ? [[alternative HTML version deleted]] __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

Re: [R] t test problem?

Hi, Many thanks for your helpful comments and suggestions. You're welcome. The attached are the data in both log10 scale and original scale. It would be very grateful if you could suggest which version of test should be used. I feel that it would be inappropriate. It depends on the origin

[R] t test problem?

Hi Liu, I'd suggest you to use non-parametric tests (see http://www.cas.lancs.ac.uk/glossary_v1.1/nonparam.html) such as: wilcox.test() in stats package pairwise.wilcox.test() in stats package and see the result tou got (significancy/non significancy) and compare it with t test result;

Re: [R] Bootstrap

did you look at library(boot) ?boot On 9/22/04 9:19 AM, nmi13 [EMAIL PROTECTED] wrote: Dear Any, Can someone please inform me, if they have a code to estimate the varaince using bootstrap resampling method under a two stage cluster design. Thanks for all your help and time.

RE: [R] t test problem?

Hi Kan Lui, I've had a quick look at the data. The logged data seems reasonably nicely distributed (roughly symmetrical + equal variance). Indeed the y variable passed the (very strict) shapiro.test for normality. However the main problem is that I do not get the same results as you for the

RE: [R] is.constant

Christian Hoffmann x - c(1, 2, NA) is.constant(x) [1] TRUE For data such as c(1, 1, 1, NA), I should think the safest answer should be NA, because one really doesn't know whether that last number is 1 or not. Andy My version is is.constant - function(x) { if

[R] asypow.noncent. Thanks!!

Thank y'all for a) pointing out what I was doing wrong b) being so patient with what, in retrospect, was such an obvious blunder by me. I'm afraid I was confused by the accompanying pdf file, which is not entirely consistent with the help files.

[R] Re: R-help Digest, Vol 19, Issue 22

I will be out of the office 9/22/04 - 9/27/04. For immediate help, please call the JPSM main number, 301-314-7911. __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

Re: [R] function to change a string to uppercase

On Wed, Sep 22, 2004 at 10:43:18AM +0200 - a wonderful day - [EMAIL PROTECTED] wrote: Is there a fonction in R to change a string to uppercase ? toupper() hint: basic functionality is in package base ;) which is easy to investigate in html-help - packages in my installation

[R] impenetrable warning

Dear R-help, Can anyone explain the meaning of the warning, Singular precision matrix in level -1, block 1 ? Or how to track down where it comes from? More precisely, using the nlme package, I'm issued with the warning itt2 - lme(lrna~rx.nrti+lbrna, random=~1|patid,

[R] capitalize etc string

Dear Any, Is there a fonction in R to change a string to uppercase ? Thanks for all your help Use the following: capply - function(str, ff) { sapply(lapply(strsplit(str, NULL), ff), paste, collapse=) } cap - function(char) { # change lower letters to upper, others leave unchanged if (any(ind

Re: [R] t test problem?

On 22-Sep-04 kan Liu wrote: Hi, Many thanks for your helpful comments and suggestions. The attached are the data in both log10 scale and original scale. It would be very grateful if you could suggest which version of test should be used. By the way, how to check whether the variation is

[R] is.constant 2

x - c(1, 2, NA) is.constant(x) [1] TRUE For data such as c(1, 1, 1, NA), I should think the safest answer should be NA, because one really doesn't know whether that last number is 1 or not. Andy My version is is.constant - function(x) { if (is.numeric(x) !any(is.na(x)))

[R] Residuals, smoothers and estimates of noise

We have a clinical measurement on patients over time. Each patient has about 5 of these measurements over a period of two years, but the measurement are not necessarily taken at equal space in time. We want to use this data to establish test-retest variability. My first thought was to look at the

[R] loops: pasting indexes in variables names

I cannot figure out how, using R, I can paste indexes or characters to the variable names which are used within loops. I will explain this with a simple example: Immagine I have a huge series of variables, each one taken two times, say x1 x2 y1 y2 z1 z2. Now, immagine that I want to compute

Re: [R] is.constant 2

Christian Hoffmann christian.hoffmann at wsl.ch writes: x - c(1, 2, NA) is.constant(x) [1] TRUE For data such as c(1, 1, 1, NA), I should think the safest answer should be NA, because one really doesn't know whether that last number is 1 or not. Andy My

RE: [R] impenetrable warning

Generally you can set options(warn=2), run the code, then do traceback(). Andy From: Simon.Bond Dear R-help, Can anyone explain the meaning of the warning, Singular precision matrix in level -1, block 1 ? Or how to track down where it comes from? More precisely, using the nlme

Re: [R] loops: pasting indexes in variables names

Hi Umberto, look at ?get, ?assign, ?paste and try the following: x1 - rnorm(10) x2 - rnorm(10) y1 - rnorm(10) y2 - rnorm(10) z1 - rnorm(10) z2 - rnorm(10) ## names. - c(x, y, z) for(i in names.){ res1 - get(paste(i,1,sep=)) res2 - get(paste(i,2,sep=)) assign(paste(d,i,sep=),

Re: [R] loops: pasting indexes in variables names

Umberto Maggiore umberto_maggiore at hotmail.com writes: I cannot figure out how, using R, I can paste indexes or characters to the variable names which are used within loops. I will explain this with a simple example: Immagine I have a huge series of variables, each one taken two times,

Re: [R] is.constant 2

On Wed, 22 Sep 2004 12:16:10 + (UTC), Gabor Grothendieck [EMAIL PROTECTED] wrote : Suggest you use an S3 generic and a separate methods for factor, and in the future, other classes. That's not a bad idea, but is it really worth the trouble? Why not piggyback on the unique() generic, and

[R] dot density maps

Dear All, In the moment i'm using the map and maptools package to read shapefiles and display the maps. I'm looking for the possibility to draw points (randomly positioned or positioned according to a grid) into the polygons instead of filling the polygons with colors. For example: a

Re: [R] loops: pasting indexes in variables names

Umberto Maggiore wrote: I cannot figure out how, using R, I can paste indexes or characters to the variable names which are used within loops. I will explain this with a simple example: Immagine I have a huge series of variables, each one taken two times, say x1 x2 y1 y2 z1 z2. Now, immagine

[R] DCOM server high quality graphics ?

I'm using c# and I can produce some graphics with the DCOM Server on R but the graphics's quality is very bad, fonts problems, anti-aliasing ? How to produce high quality graphics with the DCOM server [[alternative HTML version deleted]] __

[R] Multinomial Response Variable: Estimating the parameters

Hello! In a dyslexia-spelling-study we have a multinomial response variable with four categories. There are 12 subjects. Each subject was tested on 30 words, i.e. there are 30 responses for each subject. We are interested in estimating the parameters pi1, pi2, pi3, pi4, if possible point

[R] R course

Apologies for cross-posting We would like to announce a 3 day course: R programming for beginners. When: 10-12 January 2005 (Monday-Wednesday). Location: The Ythan hotel in Newburgh (UK). Newburgh is a small coastal village 10 miles north of Aberdeen airport. Host: Organised by Highland

Re: [R] RMySQL and Blob

Hi Jonathan, Currently RMySQL doesn't handle blob objects. The mechanics of inserting and extracting blob objects by itself is not too hard, but issues such as how should blobs be made available to R, how to prevent buffers overflows, how to prevent huge blobs from exhausting the available

RE: [R] R 1.9.1 Fails to Start on WinXP SP2

One other suggestion: Run the msconfig System configuration utility to turn off most of the software that loads on your machine at startup, and see if that allows R to start. Then gradually add it back until you find the culprit, if there is one. Bingo! You rock Duncan. I had Rage3D (a video

[R] layout for xyplot

Dear all, I tried to use layout argument in xyplot to get one panel per page. I have a dataframe named 'data' with the following variables: x, y = coords, sub, bloc = 2-level factors, etat = 5-level factor, I did : lset(theme = col.whitebg()) xyplot(y ~ x | bloc*sub , data=data,

Re: [R] R 1.9.1 Fails to Start on WinXP SP2

On Wed, 22 Sep 2004 09:28:03 -0500, Scott Higginbotham [EMAIL PROTECTED] wrote : One other suggestion: Run the msconfig System configuration utility to turn off most of the software that loads on your machine at startup, and see if that allows R to start. Then gradually add it back until you

RE: [R] t test problem?

Hi, thanks for all your suggestions. It is realy helpful. The data was not paired. Sorry for a wrong set number of sample, compared to the sample I used for tests. Best Liu Wayne Jones [EMAIL PROTECTED] wrote: Hi Kan Lui, I've had a quick look at the data. The logged data seems

Re: [R] Multinomial Response Variable: Estimating the parameters

hmm, the design sounds for me similiar to choice-data from choice-based-conjoint. IMHO check library(MNP) ? regrads,christian Am Mittwoch, 22. September 2004 15:04 schrieb Martin Plöderl: Hello! In a dyslexia-spelling-study we have a multinomial response variable with four categories.

[R] glmmPQL correlation structure

Running Mac OS 10.3.5 and R 2.0 Does glmmPQL use the same spatial dependence models as gls in nmle? It does not seem to - I get the following, for example: m2 - glmmPQL(S.Early ~ fertilized*watered, data=geodat, family=poisson, random=~1|col) iteration 1 iteration 2 plot(Variogram(m2,

Re: [R] layout for xyplot

Have you read the posting guide, which says: quote For questions about unexpected behavior or a possible bug provide details about your platform (Windows2000, Linux, OS X) and R version (type version at the R prompt). State the full version number, e.g., `1.8.1', not just `1.8'. State whether

[R] block statistics with POSIX classes

I have a monthly price index series x, the related return series y = diff(log(x)) and a POSIXlt date-time variable dp. I would like to apply annual blocks to compute for example annual block maxima and mean of y. When studying the POSIX classes, in the first stage of the learning curve, I

[R] Sample without replacement

Hello, I have a simple problem (checked the archives and the appropriate help pages, to no avail). I want to creat a vector that is length(2000). The vector is to consist of two strings( std and dev) with 1760 stds and 240 devs. Furthermore, for each element of the vector, i want the

Re: [R] Sample without replacement

Mark G Orr [EMAIL PROTECTED] writes: Hello, I have a simple problem (checked the archives and the appropriate help pages, to no avail). I want to creat a vector that is length(2000). The vector is to consist of two strings( std and dev) with 1760 stds and 240 devs. Furthermore, for each

[R] pairs, panel.functions, xlim and ylim

Hi, I have the following problem. I wanted to get a matrix of scatterplots and I used pairs. I wanted to add the line y=x in each plot and I created a panel function for this scope. I used points and abline in the following way: ## put y=x in each plot panel.lin- function(x, y) {

RE: [R] ordered probit and cauchit

The Political Science Computational Laboratory at Stanford has R code for ordered probit (courtesy of Simon Jackman): http://pscl.stanford.edu/oprobit. David Reinke -Original Message- From: David Firth [mailto:[EMAIL PROTECTED] Sent: Wednesday, September 22, 2004 02:04 To: roger koenker

Re: [R] block statistics with POSIX classes

__ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

RE: [R] RMySQL and Blob

Hi David, The application I have in mind is for images. In my case, size of images is known and they are not big. As an example, a 64*32 image will have 2048 pixels. If they are 8-bit grey-level pixels, the image occupies 2KB memory. I may venture to guess that the unknown size and type of a

[R] Issue with predict() for glm models

__ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Re: [R] glmmPQL correlation structure

Have you tried GLMM in lme4? Doug Bates is the primary architect of both nlme and lme4. Therefore, I would think that a spatial dependence model that works in nlme might also work in GLMM. hope this helps. spencer graves Martin Henry H. Stevens wrote: Running Mac OS 10.3.5 and R

[R] THanks, random/repl. prob. SOLVED.

Thanks for the help, it worked. Mark G. Orr Postdoctoral Research Fellow Dept. of Neuroscience RM 825 Kennedy Center Albert Einstein College of Medicine Bronx, NY 10461 718-430-2610 __ [EMAIL PROTECTED] mailing list

Re: [R] impenetrable warning

In this case, this trick will identify the specify command in the lme code that produce the error. If that does not lead you to an answer, have you tried simplifying your lme call to identify more clearly which part (or combination) seems to generate the problem? hope this

[R] Issue with predict() for glm models

__ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

[R] pairs, panel.functions, xlim and ylim

Valerie, A bit ugly, because you must ignore some warnings, but for me works the code below. Best Jens Oehlschlägel x - rnorm(100, sd=0.2) x - cbind(x=x-0.1, y=x+0.1) pairs(x , panel=function(x,y){ function to be called with xlim=, ylim= parameters points(x, y) abline(0,1)

[R] aparchFit()$fitted.value

Dear R people, I'm not able to have the component residuals, fitted.value from an aparchFit() estimation as explain in the Value of aparchFit Help, package fSeries. Could someone help me? Thanks in advance. Lisa __ [EMAIL PROTECTED] mailing list

[R] Issue with predict() for glm models

Hello everyone, I am having a problem using the predict (or the predict.glm) function in R. Basically, I run the glm model on a training data set and try to obtain predictions for a set of new predictors from a test data set (i.e., not the predictors that were utilized to obtain the glm

Re: [R] Issue with predict() for glm models

Dear Mark and Joe, Actually, the problem here appears to be caused by the use of a matrix on the RHS of the model formula. I'm not sure why this doesn't work (I must be missing something -- perhaps someone else can say what), but Joe can get the output he expects by specifying the columns of his

Re: [R] t-test problem

Kan Liu wrote: Hello, I got two sets of data x=(124738, 128233, 85901, 33806, ...) y=(25292, 21877, 45498, 63973, ) When I did a t test, I got two tail p-value = 0.117, which is not significantly different. If I changed x, y to log scale, and re-do the t test, I got two tail p-value

[R] Facing problems with C code compilation - Please help.

Hello, I started using R a month ago - so I am a novice in this area. I am stuck with a problem and need some help urgently. I am using windows version of R 1.9.1. I am trying to compile C code in it. I have my C code - hello.c is lying in C:\Program Files\R\rw1091 This code is - #include

Re: [R] png problem

Hi Clément Calenge wrote: Hello, Thanks for the fast reply. Paul Murrell wrote: Hi Clément Calenge wrote: Dear R-users, I have a small problem with the function png(), when used with the argument colortype=pseudo.cube. png(toto.png, colortype=pseudo.cube) image(matrix(rnorm(1), 100, 100))

[R] Cox proportional hazards model

Good afternoon, I am currently trying to do some work on survival analysis. - I hope to seek your advice re: 2 questions (1 general and 1 specific) (1) I'm trying to do a stratified Cox analysis and subsequently plot(survfit(object)). It seems to work for some strata, but not for others. I

[R] Re: Your Amazon.com Inquiry

Greetings from Amazon.com. We're sorry. You replied to a notification-only address that cannot accept incoming e-mail. But that's OK--this automated response will direct you to the right place at Amazon.com to answer your question or help you make changes to your order. To change any unshipped

RE: [R] Facing problems with C code compilation - Please help.

Read and follow the instructions in c:\Program Files\R\rw1091\README.packages _very_, _very_ carefully. Stray from it even a bit and you get what you deserve. Andy From: neha chaudhry Hello, I started using R a month ago - so I am a novice in this area. I am stuck with a problem and

[R] multinomial logistic regression

Hi, how can I do multinomial logistic regression in R? I think glm() can only handle binary response variable, and polr() can only handle ordinal response variable. how to do logistic regression with multinomial response variable? Thanks __ Y!

Re: [R] multinomial logistic regression

array chip wrote: I think glm() can only handle binary response variable, That's not true, have you looked at ?glm and ?family ? HTH, Kevin -- Ko-Kang Kevin Wang PhD Student Centre for Mathematics and its Applications Building 27, Room 1004 Mathematical Sciences Institute (MSI) Australian

Re: RE: [R] Facing problems with C code compilation - Please help.

Hi Guys, Thanks a ton for all the help. My code is finally compiling and running. -Neha - Original Message - From: Liaw, Andy [EMAIL PROTECTED] Date: Wednesday, September 22, 2004 5:12 pm Subject: RE: [R] Facing problems with C code compilation - Please help. Read and follow the

RE: [R] layout for xyplot

thanks a lot ! My version is 1.9.1. : _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major1 minor9.0 year 2004 month04 day 12 language R and I'm working on a Windows2000 platform. Besides I tried layout = c(1,1,4), but it