similar to: Incorrect order

Hello, I am a newbie in R and I have a problem. I want this line: newMyData <-MyData[!(is.na(col1)) | !(is.na(col2)) | !(is.na(col3)) | !( is.na(col4)) | !(is.na(col5)),] write in for loop. The number of col[number] may be different from 3 to 8 Is there any ideas? Thank you Laura [[alternative HTML version deleted]]

Dear Laura, It looks like you want to remove all rows for which each column is NA. You can to that with the code below. na.matrix <- is.na(MyData) all.na.row <- apply(na.matrix, 1, all) MyData[!all.na.row, ] Best regards, ir. Thierry Onkelinx Statisticus / Statistician Vlaamse Overheid / Government of Flanders INSTITUUT VOOR NATUUR- EN BOSONDERZOEK / RESEARCH INSTITUTE FOR NATURE AND

Hi, I have used this command : resamples<-lapply(1:1000,function(i) sample(lambs,replace=F)) resamples2<-lapply(resamples,Cusum) to get a list of 1000 samples of my data. The function Cumsum is defined as follows: Cusum<-function(x){ SUM<-cumsum(x)-(1:length(x))*mean(x) min<-min(cumsum(x)-(1:length(x))*mean(x)) max<-max(cumsum(x)-(1:length(x))*mean(x)) diff<-max-min

Hello, my data.frame is sort of a collection of process values, i.e. huge run-chart. It consists of a time-stamp in the first column (date as string), factors in the following columns (used for subset-filtering), and some process-data columns. Hereafter, two examples are listed, showing the problems that occour during print: At first the example, that works fine:

Dear R specialists, I have a multivariate statistics question that I want to submit to the R community (which conveys a very good statistical knowledge). I need to perform an optimal scaling based on a discrete variable and an ordinal variable. The discrete variable, Area, defines a geographical area. The ordinal variable, EducationLevel, describes the education level of individuals (the

Hi Alejo, According to my knowledge the two plots are different because in the first one a point belongs to a group depending on its group in the data whereas in the second plot a point belongs to the group predicted by the linear discriminant analysis. I hope somebody will correct me if I am wrong. Alain Alejo C.S. wrote: > Hi Alain, this is the code: > > > library(MASS) >

I'm trying to assign NAs to values that satisfy certain conditions (more complex than shown below) and it gives the right result, but breaks the loop having done the first one viz: new<-c(rep(5,4),6) for (i in 1:6) {new[new[i]>5.5][i]<-NA} gives the correct result, though an error message appears which causes a break if it's in a loop. If I can get rid of the error message and

Hi, Just a quick query - if I'm creating a function to produce a number of histograms per page of output (one per column from a matrix), how can I pass the column name of the matrix into the title (or indeed to form part of the x-axis label)? TIA, Laura Laura Quinn Institute of Atmospheric Science School of Earth and Environment University of Leeds Leeds LS2 9JT tel: +44 113 343 1596 fax:

Ok, I have reinstalled R-1.9.0 and this appears to have fixed the problems I was having with png(). However, I have a further question regarding png() Is it possible to pass a par() argument to the png() command? I am wanting to produce 4 plots per object, which I normally acheive on an X window by par(mfrow=c(1,4)). I have tried calling a new plot and setting par in this way but this has no

Hello, I requested help a couple of weeks ago creating a dipole field in R but receieved no responses. Eventually I opted to create a 3d sinusoidal plot and concatenate this with its inverse as a means for a "next best" situation. It seems that this isn't sufficient for my needs and I'm really after creating a continuous 3d gaussian mesh with a "positive" and

Dear R users, I am in trouble using the function solve(matrix) in a loop. When it happens that I have a singular matrix, I get an error message saying so and the loop stops. My question is if there in any way to keep on going in the loop. Thank you in advance, Laura *---------------------------------------------------------* Laura Antolini, PhD Department of

Dear R People: I used the "Surv" function to produce the following object: >a <- Surv(1:4,2:5,c(0,1,1,0)) a [1] (1,2+] (2,3 ] (3,4 ] (4,5+] >str(a) Error in "[.Surv"(object, 1:ile) : subscript out of bounds > Why does str(a) give an error, please? Or did I do something wrong? Thanks in advance. R Version 1.9.1 Windows Sincerely, Laura Holt mailto: lauraholt_983

??????? Original Message ??????? On Thursday, April 11, 2019 12:55 AM, John Fawcett via dovecot <dovecot at dovecot.org> wrote: > On 11/04/2019 00:51, Laura Smith via dovecot wrote: > > > ??????? Original Message ??????? > > On Wednesday, April 10, 2019 11:48 PM, John Fawcett via dovecot dovecot at dovecot.org wrote: > > > > > On 11/04/2019 00:18, Laura

Dear list, I just tried to use the function cluster.stat in the package fpc. I just have a couple of questions about the syntax: cluster.stats(d,clustering,alt.clustering=NULL, silhouette=TRUE,G2=FALSE,G3=FALSE) 1) the distance object (d) is an object obtained by the function dist() on my own original matrix? 2) clustering is the clusters vector as result of one of the many clustering methods?

Hello, Is there a simple way of stacking/merging two dataframes in R? I want to stack them piece-wise, not simply add one whole dataframe to the bottom of the other. I want to create as follows: x.frame: aX1 bX1 cX1 ... zX1 aX2 bX2 cX2 ... zX2 ... ... ... ... ... aX99 bX99 cX99 ... zX99 y.frame: aY1 bY1 cY1 ... zY1 aY2 bY2 cY2 ... zY2 ... ... ... ... ... aY99 bY99 cY99 ...

Hello there, I'm trying to plot vectors with p<0.1 in a NMDS ordination plot using p.max. Below the scripts I'm using. I guess I'm missing something! could you please give me a hand? species<-metaMDS(species_matrix)ef<-envfit(species,environmentaldata_file,permu=999,na.rm=TRUE)efplot(species, dis="sites")plot(ef,p.max=0.1) Error in plot.envfit(ef, p.max = 0.1) :

Hi dear list, If anybody could help me, it would be great! I have two files: File 1 is a list (one column and around 100000 rows) File 2 is a list with all the names from file one and a few more (one column and more than 100000 rows) What I want is to add a column in file 2 that says which name appeared in file 1 and which doesn't (yes and no would work as a code) It's very important to

Hi, Is it possible to recreate "smoothed" data sets in R, by performing a PCA and then reconstructing a data set from say the first 2/3 EOFs? I've had a look in the help pages and don't seem to find anything relevant. Thanks in advance, Laura Laura Quinn Institute of Atmospheric Science School of Earth and Environment University of Leeds Leeds LS2 9JT tel: +44 113 343 1596

Hello, Rather new to RoR, so I''m not sure about the terminology and such :-/ . My question is: how could I know the model which is associated to a controller? For example: - controller class is TestController, which is a subclass of ApplicationController - associated model class is Test How could I write some code in a method of ApplicationController to dermine the current

this is hopeless, ive given up trying to update samba, might as well go back to win server this weekend unless someone can send me a ftp mirrorlist or url to samba.org page with this list On Fri, Dec 13, 2019 at 8:29 PM Rowland penny via samba < samba at lists.samba.org> wrote: > On 13/12/2019 10:24, Laura Steynes wrote: > > i have full speed from the mirror, it is samba.org