Displaying 20 results from an estimated 500 matches similar to: "Formatting a Table"
2008 Mar 25
1
Subset of matrix
Dear R users
I have a big matrix like
6021 1188 790 290 1174 1015 1990 6613 6288
100714
6021 1 0.658 0.688 0.474 0.262 0.163 0.137 0.32
0.252 0.206
1188 0.658 1 0.917 0.245 0.331 0.122 0.148 0.194
0.168 0.171
790 0.688 0.917 1 0.243 0.31 0.122 0.15 0.19
0.171 0.174
290 0.474
2008 Jan 30
4
Loop with variable index
I have a list of 20 values. The first time through a loop I want to find the
mean and stnd.dev. of the first two values; the second time through the loop
I want to find the mean and stnd. dev. of the first 3 values, etc. until
the last time through the loop I want to find the mean and stnd. dev. of all
20 values, so I end up with 19 means and stnd. deviations.
How would I construct such a loop?
2018 May 15
0
Systemfit
... and the mailing list is picky about attachments... whatever you attached did not conform to the stringent requirements mentioned in the Posting Guide. Pasting the code right into the email is usually safest, though you DO have to post using plain text (as the Posting Guide indicates) or your code may get mangled by the automatic html format removal.
On May 15, 2018 7:04:31 AM PDT, Bert Gunter
2007 Aug 10
7
Help wit matrices
Hello all,
I am working with a 1000x1000 matrix, and I would like to return a
1000x1000 matrix that tells me which value in the matrix is greater
than a theshold value (1 or 0 indicator).
i have tried
mat2<-as.matrix(as.numeric(mat1>0.25))
but that returns a 1:100000 matrix.
I have also tried for loops, but they are grossly inefficient.
THanks for all your help in advance.
Lanre
2008 Feb 26
4
numeric format
Hi!
I'm an R newbie and this should be a trivial problem, but I can't make it
work and cannot find what I'm doing wrong in the literature.
I entered the the command:
table<-data.frame(x, scientific=F, digits=4)
table
This prints a column of x with 16 useless decimal places after the decimal
point. Also, it prints an unwanted index number (1-20) in the left column.
How do I get rid
2018 May 15
1
Systemfit
Unless there is good reason not to, always cc the list -- there are lots of
smarter folks than I on it who can help.
I may or may not have time to look at this. Hopefully someone else will.
-- Bert
Bert Gunter
"The trouble with having an open mind is that people keep coming along and
sticking things into it."
-- Opus (aka Berkeley Breathed in his "Bloom County" comic strip
2009 Feb 23
1
why results from regression tree (rpart) are totally inconsistent with ordinary regression
Hi,
In my analysis of impacts of insecticide-treated bednets on malaria, I
look at the relationship between malaria incidence and mosquito
behaviors. The condensed data set is copied here. Ordinary regression
(lm) shows that Incidence was negatively related to Mortality. This
makes sense because the latter reflected the strength of killing
mosquitoes by insecticide-treated nets. Since the
2018 May 15
2
Systemfit
OK, Let's try this again! Here is the reproducible script; it is long because I had to copy the panel dataset here. My question is related to systemfit; I don't know how to get the result for the entire panel.
#Reproducible script
Empdata<- read.csv("/Users/ngwinuiazenui/Documents/UPLOADemp.csv")
View(Empdata)
install.packages("systemfit")
2006 Mar 08
1
RES: survival
Dear Thomas,
The head of my dataset
> head(wsuv)
parcel sp time censo treatment
species
1 S8 Poecilanthe effusa ( Hub. ) Ducke. 1 1 1 1
2 S8 Poecilanthe effusa ( Hub. ) Ducke. 1 1 1 1
3 S8 Poecilanthe effusa ( Hub. ) Ducke. 1 1 1 1
4 S8 Poecilanthe effusa ( Hub. ) Ducke. 1 1 1
2003 Jan 20
1
make check for R-1.6.2 on IBM AIX
Dear all,
The 'make check' step fails for the pacakge mva on IBM AIX.
The tail of the Rout log file looks like:
> for(factors in 2:4) print(update(Harman23.FA, factors = factors))
Call:
factanal(factors = factors, covmat = Harman23.cor)
Uniquenesses:
height arm.span forearm lower.leg weight
0.170 0.107 0.166
2008 Jan 28
1
R loops
I'm a new user and am having trouble with loops.
In the following, I'm trying to add the results of "test" and the loops are
not working.
I've simplified the loop. What am I doing wrong?
Thanks!
> test<-numeric(20)
> tot<-numeric(20)
> for(i in 1:20){test[i]<-1}
> for (i in 1:20){tot[i]<-(test[i]+tot[i])}
> tot
[1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
2008 Feb 02
1
Plotting 3 vectors on one graph.
I'm an R newbie and am trying to plot 3 vectors, say a,b,c. I have
downloaded 3 R manuals and searched your forum. There are plenty of X vs Y
examples, but cannot find how to plot 3, or more vectors one one graph. I'm
sure I overlooked something.
Thanks for any help.
CHV
--
View this message in context: http://www.nabble.com/Plotting-3-vectors-on-one-graph.-tp15236552p15236552.html
2009 May 19
1
Wilcoxon nonparametric p-values
When I use wilcox.test, I get vastly different p-values than the problems
from Statistics textbooks.
For example:
The following problem comes from "Applied Statistics and Probability for
Engineers", 2nd Edition, by D. C. Montgomery. Page736, problem 14.7. The
problem is to compare the sample data with a population median of 8.5. The
book answer is p = 0.25, wilcox.test answer is p =
2011 Mar 12
3
betareg help
Dear R users,
I'm trying to do betareg on my dataset.
Dependent variable is not normally distributed and is proportion (of condom
use (0,1)).
But I'm having problems:
gyl<-betareg(cond ~ alcoh + drug, data=results)
Error in optim(par = start, fn = loglikfun, gr = gradfun, method = method, :
initial value in 'vmmin' is not finite
Why is R returning me error in optim()?
What
2009 Jan 14
1
Can't find files after install package (Windows)
I am a total newbie at R but experienced with computers. If this is not the
right forum for this question, please let me know one that is. I searched in
the R manuals in the Help section, the Nabble help and Rseek but no luck.
I am trying to install the package gcl on my recently installed R 2.8.1
I searched on cran for the windows version and found a zip file for the
package
In RGui, I ran
2018 May 16
0
Systemfit
Sadly you failed to set your email program to send plain text and the data is corrupted at my end.
I also think you need to reduce the size of the data set... the intent here is to increase your understanding, not debug your particular analysis.
I will say that I am having a very challenging time understanding what you are trying to accomplish though. What are the equations that you think need
2009 Jul 09
2
How to Populate List
Hi,
I'm new to R and would like to know, how one can populate the list with array data.
I'm reading a tab separated table in R. The data in the table looks something like this.
#Table Data
Comp A B C
Extracellular 103 268 535759
Nucleus 45603 47783 442744
#R code
myData <- read.table("table.data",
header=T,
2008 Aug 21
1
summary.lme and anova question
Dear all,
When analyzing data from a climate change experiment using linear mixed-effects models, I recently
came across a situation where:
- the summary(model) showed a significant difference between the levels of a two-level factor,
- while the anova(model) showed no significance for that factor (see below).
My question now is: Is the anova.lme() approach correct for that model? And why does
2011 Apr 05
6
simple save question
Hi,
When I run the survfit function, I want to get the restricted mean
value and the standard error also. I found out using the "print"
function to do so, as shown below,
print(km.fit,print.rmean=TRUE)
Call: survfit(formula = Surv(diff, status) ~ 1, type = "kaplan-meier")
records n.max n.start events *rmean *se(rmean) median
200.000
2008 Aug 08
2
aggregate
Dear All-
I have a dataset that is comprised of the following:
doy yr mon day hr hgt1 hgt2 hgt3 co21 co22 co23 sig1 sig2 sig3 dif flag
244.02083 2005 09 01 00 2.6 9.5 17.8 375.665 373.737 373.227 3.698 1.107
0.963 -0.509 PRE
244.0625 2005 09 01 01 2.6 9.5 17.8 393.66 384.773 379.466 15.336 11.033
5.76 -5.307 PRE
244.10417 2005 09 01 02 2.6 9.5 17.8 411.162 397.866 387.755 6.835 5.61
6.728